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AUTHOR'S  NOTE 


CHAPTER  II,  "Tooting  of  Hypotheses , "  ia  part'  of  the  material 
being  prepared  in  connection  with  a  contemplated  handbook  or  monograph 
on  Statistical  Techniques  in  Life  Testing:  It  is  in  the  nature  of  a 
preliminary  report  on  ono  aapect  of  the  over-all  tinder  taking. 

Earlier  reports  were: 

Technical  Report  Ho*  1,  "Statistical  Developments  in  Life 

Testing,"  June  1,  195?. 

Technical  Report  No*  R,  "Tho  Exponential  Distribution  and  Its 

,  Role  in  Life  Testing,"  May  1,  1958  o 

"An  On lino  of  Throe  Chapters  of  a  Handbook  on  Statistical 
Methods  in  Life  Testing,"  June  5,  1958.. 

Further  material  dealing  with  other  aspects  of  life  testing 
is  in  preparation. 

Comments  and  suggestions  are  invited. 


Benjamin  Epstein 


TABLE  07  CONTENTS 
CHAPTER  2 
TESTING  OF  HYPOTHESES 

Pag* 


INTRODUCTORY  REMARKS*  2.X 

SECTION  I 

(1)  A  lif*  test  procedure,  bac«d  on  th*  first  r  failures 
with  th*  property  that  th*  probability  of  rejecting  a 

lot  with  mean  life  9«9 is  equal  to  oC  »  2.2 

(2)  Cost  considerations  in  choosing  number  of  items  to  be 
placed  on  test  for  a  preassigned  number  of  failures  r 

to  be  observed.  2*5 

(3)  Life  tests  based  only  on  the  rth  failure  for  a  sample 

of  size  a  .  2*7 

(a)  Exact  procedure  2.7 

(b)  Approximate  procedure  2.8 

(4)  Numerical  examples.  2.10 


SECTION  II 

(1)  A  life  test  based  on  the  first  r  out  of  n  failures, 
having  the  property  that  its  O.C.  curve  is  such  that 


L(90)>  l-<*  and  L(9X)<0  .  2.13 

(2)  Numerical  examples.  2.13 

SECTION  III 

(1)  A  life  test,  truncated  in  time,  having  the  property  that 

its  O.C.  curve  is  such  that  L(9  )>  1-oC.  ,  L(9.)<^  .  2.15 

o  x  * 

(2)  Truncated  life  tests  when  one  is  interested  in  the 

probability  of  failing  before  time  To  (Remarks  2  and  3).  2.19 

(3)  Numerical  examples.  2.22 

SECTION  IV 

(1)  Sequential  life  tests.  2.24 

(2)  Graphical  procedures  for  sequential  life  tests.  2.25 

(3)  Numerical  examples.  2.27 

APPENDICES 


(2A)  Proofs  of  fundamental  results  in  life  testing  baaed  on 
the  exponential  distribution.  Proofs  of  results  in 
sections  I  and  II. 


(2B)  Alternative  estimates  for  9  . 

Bibliography  for  sections  I  and  II 


2.33 

2.41 

2.43 


TABLE  OF  CONTENTS  (contd) 

Pag* 

(2C)  Detailed  dlaeuaslon  of  tnme&ted  replacement  and 

nonreplaeement  testa.  2.44 

(2D)  Tests  of  the  fora  9  >  C  considered  as  truncated 

rtn 

tests.  2.44 

(2E)  Illustration  of  material  given  in  section  III  and 

Appendix  (2D).  2.49 

Bibliography  for  Section  III.  2.52 

(2F)  Proofs  of  results  in  section  IV.  2.53 

(20)  Truncated  sequential  life  tests.  2*54 

(1)  Truncation  on  the  nuaber  of  items  failed.  2.54 

(2)  Truncation  on  the  total  observed  life.  2.60 

(3)  Truncation  on  both  the  nuaber  of  items  failed 

and  total  observed  life.  2.67 

(2H)  Probability  of  termination  of  the  sequential  life  test 

at  p reassigned  values  of  rQ  and  VQ  .  2.70 

(21)  Upper  and  lower  bounds  for  L(9)  and  EQ(r)  .»  2. % 

(2J)  Exact  calculations  of  L(9)  and  E0(r)  <•  2.77 

(2K)  An  approximate  formula  for  Eft  (t)  in  the  nonroplacemoat 

case.  2.80 


TABLES 


Values  of  (2r)/2r 


Values  of  0  accepted  with  probability  p  when  9*1 
is  accepted  with  probabilities  .99*  .95,  <90,  .75  and 
.50  . 


2.82 


2.84 


3(a).  Values  of  E(X„  )/9  whore  is  the  rth  smallest 

r»n  r.n 

value  in  a  random  sample  of  size  n  drawn  from  a 

1  -x/9 

distribution  whose  probability  density  function  is  g  e 

for  r«l(l)n  ,  n=ia)20(5)30(10)100,  2.89 

3(b)  Values  of  E(X _ )/©  where  X_  _  is  the  rth  smallest 

PfH  rto. 

value  in  a  random  sample  of  size  n  drawn  from  a 
distribution  whose  probability  density  function  is  ^ 


for  r*l  (1)10(5)50(10)50(25)100  and  n-kr  ,  with 
k«l  (1)10(10)20  . 


2.98 


TABLE  OF  CONTENTS  (contd.) 


Pa^e 

3(c)  Ratio  of  Expected  Waiting  Tino  to  oboorvo  the  rth 

failuro  in  samples  of  size  n  and  r  respectively.  2*99 

4.  Values  of  X?  w-  (2r)E(X  J/2r  for  oC  =  .01, 

A*  T  |  a 

•05i  «10f  and  ®50  *  2*100 

5.  Values  of  r  and  Xi-oe  (2r)/2r  such  that  tho  acceptance 

region  0  >9  "T?  ^.(2r)/2r  is  ouch  that  L(9  )=l-cC 

r  ,n  oAi-<\  o 

and  .  2.105 

6.  Values  of  n  ,  the  sample  size  needed  in  truncated 

replacement  procedures.  2.107 

7.  Values  of  n  ,  the  oample  oize  needed  in  truncated 

nonreplacoment  procedures.  2.111 

8.  Values  of  r  and  of  X?  _,(2r  )/2  ouch  chat  the 

o  l-oc  o 

test  based  on  using  a  sampling  plan  with  sample  size 

equal  to  CXi-oc.^ro^^o^  with  rejection  number 

r  will  have  an  O.C.  curve  such  that  L(p  )=1-<?C  and 
o  o 

</3  .  2.113 

9.  Values  of  hQ,hj,  and  ®  for  various  values  of 

oC  j  ft  5  and  .  2.114 

10.  Approximate  values  of  Eg(r)  for  sequential  tests  for 

9 

various  values  of  k«  o/0^  ,  OC  and  .  2.115 

FIGURES 

A  2 

1.  Operating  Characteristic  Curves  of  ®r  a  >  Ai-oc,^r^^r 

Values  of  l-o^=  *99«  *95 1  «90»  .75  and  .50.  2*116 

"  2.  Graphical  solution  to  example  1  on  sequential  life  tests  2.136 

3.  Graphical  solution  to  example  2  on  sequential  life  tests.  2.137 

4(a)  Comparison  of  Efl(r)  curves  for  sequential  and  truncated 

plans.  2*138 

4(b)  Comparison  of  Eg(t)  curves  for  sequential  and  truncated 
plans. 


2.139 


TABLES  Of  CONTENTS  (contd.) 


FIGPRES  contde 

5*.  The  truncation  on  V(t)  induced  by  truncating  the  life 

testa  at  rar 

o 

6c  The  truncation  on  r  Induced  by  truncating  the  life 

test  at  7( t)fV  , 

o 

7c  The  truncated  life  test  induced  by  truncating  at  r«rQ 
and  V  ( t )  scV 

o 


2®l40 

2,141 


2  d42 


2.1 


Chapter  II 

Testing  of  Hypotheses 

Introductory  remarks:  In  the  following  we  shall  assume  that  the  under¬ 
lying  p.d.f.  of  the  lif a-time  X  is  described  by 

(1)  f(s?©)  *  ~  e'*/®  ,  x  >  0  ,  0  >  0  . 

Our  object  is  to  test  various  hypotheses  regarding  the  mean  life 
0  on  the  basis  of  censored,  truncated,  or  sequential  procedures.  Various 
tables  and  graphs  are  given  and  in  the  appendix,  appropriate  reference^ 


and  proofs  appear. 


Problem  dive  a  censored  life  teat  prooadura  which  will  haw*  the  proparty 
that  tha  probability  of  rejecting  a  lot  with  maan  Ufa  G  ■  ia  equal 
to  ^  ,  Furthermore  giro  tha  operating  oharacteristic  (O.C„)  curve  for 

thia  procedure;  i.e.,  plot  L  (0),  the  probability  of  accepting  a  lot 
having  naan  life  9.  against  0, 

Solution :  A  censored  life  teat  involves  terminating  the  test  after  a 

preaesigned  number,  r,  of  failures  occur.  More  precisely,  n  items 

are  drawn  at  random  from  a  distribution  whose  p.d«fu  is  given  by  (1)  and 

placed  on  life  test.  Observations  become  available  in  order;  i.e., 

x..  _  <  x,  <  . . .  <  x„  <  « . .  <  x„  „  ,  where  by  x.  _  is  meant  the  time 
i)Q  c« a  *•  “  r^n  *•  **  ntn  1  it n 

when  the  i'tn  failure  occurs.  Experimentation  is  terminated  as  soon  as 

the  r'tn  failure  occurs* 

In  the  non-replacement  case  (where  failed  items  are  not  replaced),  it 
can  be  shown  that  an  estimate  of  Q  which  is  "best"  in  the  sense  that  it 
is  maximum  likelihord,  unbiased,  minimum  variance,  efficient  and  sufficient 
ia  given  by 


In  the  replacement  ease  (where  one  immediately  replaces  a  failed  item 
by  a  new  one)  the  approjrlate  "best"  estimate  is  given  by 


>*  nx  /r, 
rvn  fvq 


where  by  x  is  meant  the  total  tine  (measured  from  the  beginning  of  the 
r  ,n 

life  test)  to  observe  the  r ' th  failure  and  where  the  sample  sine  n  is 


maintained  throughout  the  life  teste 

The  following  results  are  very  useful: 

(i)  The  p.d.f.  of  3  _  in  either  the  replacement  or  non-replacement 

r^o 


case  is  given  by 


2.3 


<*>  fr( y)  «  (r/o)V~l  o'"'770  »  y  *  o 

*  Oc  elsewhere 

and  further  the  random  variable  W  ■  2r$  /0  is  distributed  ae  chi-square 

r9n 

a 

with  2r  degrees  of  freedom  (X£(2r))» 

(ii)  The  expected  waiting  time  for  the  r'th  failure  is  given  by 

(5)  £(X  „)  *  0  ZZ  V 

*  »n  j»l 

in  the  non-roplacement  case  and  by 


(6)  EU  )  ■  rS/n 

in  the  replacement  case. 

From  (4)  we  can  now  write  down  a  te3t  procedure  having  the  required 
property  that  the  probability  of  rejecting  a  lot  with  mean  life  C  -  0Q 
is  equal  to  <=<  (such  a  procedure  can  also  be  eaid  to  have  size,  type  I 
error,  or  producer’s  risk  )„  The  region  of  acceptance  is  given  by 


(7)  *r,B  >  0  .  0o  (2r)/2rt 

where  we  define  the  constant  (2k)  by  the  equation 

(8)  Pr  (X2<2k>  >^(2k))  »  #  . 

The  O.Co  curve  associated  with  this  procedure  ie  givon  by 


(9)  UQ)  -  Prob($r  n  >  ®0Xi _*(2r  )/2rl«) 

=  Prob(  %Z(2v)  >  X^^(2r)), 

It  is  convenient  to  cbooee  units  in  such  a  way  that  0Q=  1..  If 
this  is  done,  (7)  and  (9)  become, respectively. 


®r.n  >  X  L*  (2r)/2r 

L(0)  =  Pr  (X2(2 r)  >Xi^  (2r)/0)o 

we  date  giro  the  value®  of  jC^^,(2r)/2r  for  r  e  1(1)10(5) 
for  *  o01„  *05*  olO,  025;>  and  o50«  The  untabulated  value® 
corresponding  to  r  <  100  can  be  found  from  tables  of  chi-square*  For 
value®  of  r  >  100  0  the  normal  approximation  to  chi-square,  or  the  more 
refined  approximations  duo  to  Fisher  or  to  'Vllson  and  Hilferty  can  be  usod. 

In  Tables  2  (a)*(b)>(o)J(d)>(e)  we  give  the  values  of  0  which  are 
accepted  with  probability  p  x  <,99,  „95s  o90t  ,?5,  <>50c  „25,  o10,.  <.05*  and 
oOl  if  a  lot  with  0  =  1  is  accepted  with  probability  *99.  »95t  <*90„  <>75  and 
*50 ^respectively*  0.  C„  curves  baaed  on  Tables  2  (a),  (b)(  (c),  (d)s  and 
(e)  are  given  for  the  r  value®  noted  above  in  Figures  1  (a)  through  1  (e)„ 
All  of  the  teat  procedure®  given  by  (7)  ha^e  the  property  that 
L(0q)  =  1-<*C  independently  of  n.  Changes  in  n  affect  E(Xr  n),  the 
expected  waiting  time  to  observe  the  r'th  value  in  a  sample  of  eize  n* 

The  appropriate  choice  of  n  for  a  given  r  and  hence  fixod  type  I  error, 
oc  ,  depends  on  economic  consideration®  and  involves  balancing  the  cost 
of  increasing  n,  with  the  gains  due  to  decreasing  the  expected  waiting 
time  of  the  experiment.  To  facilitate  the  making  of  such  judgments  we 
rive  values  of  E(X  )  for  r  *  l(l)n  and  n  *  1(1)20(5)30(10)100  in 

°  7|U 

Table  3  (a)  and  for  r  »  1(1)10(5)30(10)50(25)100  and  a  *  kr,  with 
k  =  1(1)10(10)20  in  Table  3  (b). 

It  is  particularly  interesting  to  compare  tbe  following  two  procedures: 

(1)  a  test  based  on  8  where  r  items  are  placed  on  tost  and  where 

r  ,r 

one  waits  for  all  r  items  to  fail  and  (2)  a  test  based  on  ®r^n  whers 
n  items  are  placed  on  test  and  where  one  waits  only  for  the  first  r 


(?*) 

and 

(9S) 

In  Table  1 
30(10)50(25)100 
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failures  to  occur »  From  (5)  the  expected  waiting  tinea  in  the  non-replace- 

*ent  case  a re  given  by 

/  r  r 

E(X  *  0  1/j  and  E(X  )  «*  0  /  l/n-j^l  „  respectively* 

*  r«n  j«l 

Thus  the  ratio  c-  =  E(X  )/E(X  )  is  a  measure  of  the  expected 

r.  n  r 5,  n  r  r  r 

saving  in  time  due  to  using  the  second  procedure  rather  than  the  first 

procedure*.  Values  of  ^  -  can  bo  computed  readily  from  Tables  3  (a) 

rvn 

and  3  (b).  A  brief  table  is  givon  in  Table  3  (c)„  It  follows  from  (6) 

that  the  expected  waiting  times  in  the  replacement  case  are  given  by 

E(X  )  =  0  and  B(X  )  *  rJ/n*  In  this  caso  the  ratio  is  a  imply 

*  'u  *  I  <i  U  X*  f  n 

s-C  r  »  S<x,.  )/E(X  )  tr.  -i-  „ 

rtn  r,n  r9r  n 

Numerical  Example;  Let  ue  compare  the  average  length  of  time  needed  to 
observe  (a)  the  failure  of  the  first  2  out  of  k  items  under  test  with 
the  average  length  of  time  required  to  observe,  (b)  the  failure  of  2  out 
of  2  items*.  The  answer  is  c-c  p  ^  *  E(X2  ^)/E(X2  S)  “  -y>  /  \  m  7/l8  *  .3809  ° 


Hence  it  will  take  on  the  average  only  7/l8  as  long  to  observe  the  first 

2  out  of  4  items  as  2  out  of  2*  If  life  is  exponential  so  know  that  tests 

based  on  either  (a)  or  (b)  have  the  same  0.  C„  curves;  however,  the  time 

required  for  (a)  is  on  the  average  substantially  shorter  than  for  (b) 

[average  tine  for  (a)  is  about  40T.  of  the  average  time  required  for  (b)],, 

Remark?  It  was  noted  above  that  for  given  type  I  error  Cor  producer’s 

riak)e<  and  etoppiag  numb or  r.  all  test  procedures  (7)  have  the  same 

0o  C n  curvo  independent  of  n»  >Ve  wish  to  give  a  method  for  choosing 

a  best  procedure  from  this  class  of  procedures*  It  is  clear  that  for 

fixed  «?'-  and  r,  increasing  a  will  on  the  one  hand  cut  the  expected 

waiting  time,  hut  will,  on  the  other  hand,  increase  the  cost  due  to  placing 

of 

•aore  items  on  test*  More  precisely,  if  is  the  cost^waiting  per  hour 


2.6 


and  c9  is  the  cost  of  placing  an  item  on  test,  then  the  total  expected 
cost  associated  with  a  plan  based  on  (7),  assuming  that  the  mean  life  Is 

close  to  ©o  ,  is  given  by  c^  0  0(|  +  5 =T  +*“+  H^TT)  +  c2»  ln  the 
non -replacement  case  and  c^  Q  £  +  c^n  in  the  replacement  case.  It 
is  clear  that  there  exists  an  n  which  makes  this  quantity  a  minimum  and 
that  this  minimum  depends  only  on  the  ratio  Cg/e^  Table  3  (a)  is  useful 
in  reaching  a  conclusion  as  to  the  best  n. 


To  illustrate  the  point  in  the  non-replacement  situation  consider  a 


case  where  r  ■  10,  0Q  •  1000  hrs,  c^  ■  $1  per  hour, 

tested.  Then  we  have  the  following  table: 

Cost  of  items 

n  Expected  cost  due  to  waiting  tested 

C2  =  $100  p 

Total  cost 

10 

2929 

1000 

3929 

11 

2020 

1100 

3120 

12 

1603 

1200 

2803 

13 

1346 

1300 

2646 

14 

1168 

1400 

2568 

15 

1035 

1500 

2535 

16 

931 

1600 

2531 

17 

847 

1700 

2547 

The  minimum  is  attained  for  n  *  16. 

It  is  easily  verified  that  for  the  given  values  of  r,  Qq,  c and 
c^,  the  optimum  n  is  10  in  the  replacement  case. 

We  have  Just  given  a  numerical  example.  It  is  of  Interest  to  give  a 
general  method  for  finding  the  optimum  n.  In  the  replacement  case  the 


optimum  n  is  the  Integer  nearest  to 


\  ci  9° r 

\  c2 


+  1 
* 


or  either 


m  or  m  +  1 


“  \ 


eo  r 


+  1 
1 


•  2m+l 


for  some  integer  m  ) 
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j  In  the  non-replacement  case,  we  choose  the  smallest  n  ,  such  that 
E(Xr  Q)  -  E(Xp>n+1)  ^  c  '  In  the  abov®  example,  c2/ c^  Q0  *  .1 

and  E(X10>n)  -  B(X10>n+1)  >.l  for  n  *  10,  11,  12,  13,  14,  15,  and 


l(Xio  q)  -  e(X1q  n+1)<.l  for  n>l6.  Therefore,  the  appropriate  n 


16. 


The  procedures  descrihed  hy  (7)  depend  on  knowing  the  first  r  values 
x1  <.  x^  <  . . .  <  x^.  It  is  interesting  that  if  the  underlying  distribution 
is  exponential,  then  one  cam  find  a  truncated  procedure  having  almost  pre¬ 


cisely  the  same  O.C.  curve  as  the  5  procedure.  The  only  requirement 

r,n 

that  this  be  so  is  that  n  be  moderately  larger  than  r.  The  advantage 
of  such  procedures  is  that  they  depend  only  on  an  extremely  simple  observa¬ 
tion,  the  time  of  failure  of  the  r'th  item.  Procedures  of  the  form, 

accept  if  i  „  >  T  and  reject  otherwise,  are  very  simple  to  interpret. 
r,n 

In  words,  the  rule  of  action  is  to  stop  experimentation  at  min  -  [x  ;T] 

r,n 

with  acceptance  of  the  hypothesis  if  min  [x  ;T]  »  T  (since  in  that 

r,n 

case  the  r'th  failure  occurs  after  time  T)  and  with  rejection  of  the 

hypothesis  if  min(x  ;T)  ■  x  (since  in  that  case  the  r'th  failure 

r>n  r,n 

occurs  prior  to  time  T) .  Experimentation  is  actually  truncated  at  T 
and  at  x^q  ,  respectively,  in  the  two  situations. 

To  derive  the  truncated  procedures  one  proceeds  as  follows: 

In  the  replacement  case  (7)  becomes 

(10)  >®,  X  C2r)/8n 

and  (7^)  becomes  (when  is  normalized  as  l)  , 

(101)  Xj.#n>  X  (2r)/2n. 

In  the  non-replacement  case,  the  exact  procedure  is  obtained  by  first 
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'1 

I 

In  the  non-replacement  case,  we  choose  tho  smallest  n  ,  such  that 

c_ 

B(X  )  -  E(X  ,)  <  -—4-  .  In  the  above  numerical  example,  c?/c.Q  *  .1 
v  r,n  r,n+l  ci  ®Q  d  i  o 

and  E(X. „  )  -  E(Xin  .)  >  .1  for  n  «  10,  11,  12,  13.  1^»  15,  ®nd 

E(X,  )  -  E(X,  „  ,)  <  „1  for  n  >  16  *  Therefore,  the  appropriate  n  •*  l6  . 

10, n  lo,n+j  ~ 

The  procedures  described  by  (7)  depend  on  knowing  the  first  r  values 
at,  <  x  <  ...  <  x  «  It  is  interesting  that  if  the  underlying  distribution 
is  exponential,  then  one  can  find  a  truncated  procedure  having  almost  pre¬ 
cisely  the  saute  OcC.  curve  as  the  Or  Q  procedure.  The  only  requirement 
that  this  bo  30  is  that  n  be  moderately  larger  than  r  .  Tho  advantage 
of  such  procedures  is  that  they  depend  only  on  an  extremely  simple  observa¬ 
tion,  the  time  of  failure  of  the  r’th  item.  Procedures  of  the  form, 

accept  if  x  >  T  and  reject  otherwise,  are  very  simple  to  interpret. 

. 

In  words,  the  rule  of  action  is  to  stop  experimentation  at  min 

with  acceptance  of  the  hypothesis  if  min  |jc^  n»T3  *  T  (since  in  that 

case  the  r’th  failure  occurs  after  time  T)  and  with  rejection  of  the 

hypothesis  if  min(xr  a;T)  *  xr  n  (®ince  ^Xl  ^kat  caBe  the  r’th  failure 

occurs  prior  to  time  T).  Experimentation  is  actually  truncated  at  T 

and  at  x  ,  respectively,  in  the  two  situations. 
r,n 

To  derive  the  truncated  procedures  one  proceeds  as  follows: 

In  the  replacement  case  (7)  becomes 

(10)  xr  B  >  90 

and  (71)  becomes  (when  0Q  is  normalized  as  1)  , 

(101)  *rfB  >  Xt*  (2r)/2n  • 

la  the  non-replacement  case,  the  exaot  procedure  is  obtained  by  first 
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writing  down  the  distribution  of  tha  r'th 

size  a.  In  the  case  where  ©  *  lt  the  p. 

o 


(11) 


n! 


*r,n(3t)  3  <r-i7rtn-7rr  Cl~*  J 


Therefore,  Pr(X  >  T)  becomes 
*  •}** 


smallest  value  in  a  sample  of 


d. f *  of  x  „  is  given  by 
r?n 


■(n-r+l)x 


x  >  0, 


(12)  Pr(X  >  T)  » 

rvn 


fV  (x)dx 

1  al  ri  o^n*r+1^ 

(±-T?T(Z7YT  [1"e  3  a  * 


Lotting  ti  =  l-o 


-x 


a©  get 


(13)  Pr(X  >  T)  =  n(n~*)  J  u3vl(l-u)n'’'r  du  , 

r,n  vr-l 

(l-o  ; 

If  we  want  to  solve  for  1'  in  the  equation  '^(X  >  T)  =  l-0*-  v 

wa  can  do  this  readily  from  tables  of  tha  incomplete  Beta  distribution* 
.in  alternative  procedure  is  to  use  tables  of  the  cumulative  binomial 
distribution.  If  ono  definos  B(k;n,p)  as 


k 

(14)  B(k;u;p)  =  2 _  b(V;n,p),  whe.rs 

V-o 

b(V;n,p)  »  (y)  p^(l-p)n~^,  than  one  can  show  that  (13)  becomes 

(15)  l-3(r-l;n,p)  *■  n(^^)  f  ur“1(l-u)n  rdu, 

m  ® 

where  p  =  1-e  .  Tc  solve  the  pt oblem  at  hand,  the  procedure  ic  to  compute 

p  so  that 

(16)  l-B(r-l;n,p)  =  /  _  _  (£)  pk(l-p)n  k  Tt°<- 

k*r 

and  than  sot  T  =  los(^^j)'' 

An  alternative  approximate  procedure  in  tho  non-replacement  case 
which  has  been  shown  to  be  extremely  close  to  the  exact  procedure  is 


obtained  as  follows: 


2.? 


In  (7)  or  (7*)  replace  the  unbiased  estimator  Q  by  /Br  n  xr  n 

where  n  *  1/E(X  q) 0  It  haa  been  observed  that  if  n  is  sufficiently 

larger  than  r,  then  the  O.C.  curve  resulting  from  using  the  acceptance 
region  #  >  C  virtually  coincides  with  the  O.C.  curve  associated 

with  n  >  C.  Consequently  a  test  procedure  having  virtually  the  same 

9 

O.C.  curve  as  (7)  or  (7  )  is  given  by 

(17)  *r,„  >  T  .  0oXU<2r)/2r^r,a  -  s.XSU(2r>E«r,n>'2r 

or 


<17'>  zi\n  >  XjU<2r)/2r4,a  * 


in  the  normalized  case  0  =  1. 

o 


In  Table  4  we  give  values  of  Y?  (2r)E(X  )/2r  for  the  values 

rvn 

of  r  and  «*■  covered  in  Table  1  and  for  various  values  of  n.  If  r 

and  n  are  sufficiently  large,  then  ^  log(--~-)  „  This  approxi- 

^rtn 

station  is  useful  for  computing  some  values  of  T  when  r  and  n  lie 
outside  Table  4„ 

Remark:  An  approximate  formula  for  p*  satisfying  T  (°)p .  k(l-p  -  )n”^» 

k«r  K  * 

is  found  by  solving  for  p  in  the  equation: 

oC 

log(._^  *■■)  »  Y?  (2r)E(X  )/2r.  Solving  for  p  one  gets: 

,  - 

!•••  • 


p  ° 

OC 

Example :  r  *  10  5  n  *  20.  Find  p^  such  that 


k»10 


<?>»« 


20»k 


for  «  .01,  .05,  <10,  .25,  -50-  In  column  2  we  give  the  approximate 
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vadue  of  ,aa  obtained  by  the  formula,  and  in  column  3  wo  give  the 

value  of  p  computed  by  interpolation  from  the  Harvard  tables  of  the 
oc 

Cumulative  Binomial  Distribution. 


©c 

(approximate) 

p  (exact) 
ec 

.01 

.241 

.239 

.05 

.305 

.302 

.10 

.340 

.338 

o25 

>404 

.402 

.50 

.476 

.475 

Numerical  Examples: 

1.  Find  a  censored  life  test  which  will  accept  a  lot  having  a  mean  life 
of  1000  hours  with  probability  .90.  The  experiment  is  to  be  stopped 
after  one  has  observed  the  first  5  failures. 

Solution:  In  terms  of  the  notation  that  we  have  used  0q  =  1000, =  .10, 

and  r  =  5«  In  the  non-replacement  case  =  [x^+X2+Xj+x^+(n-4)x,.]/5s1 

and  in  the  replacement  case  =  nx_/5<>  The  region  of  acceptance  is 

P»n  > 

given  by  ^  >  9Q^^(2r)/2r  *  (1000)  (10)/10 

e  436.5. 

In  words,  place  n  items  on  test.  Wait  until  the  first  5  failures  occur. 
Compute  .  Accept  the  hypothesis  that  9  *  1000  if  „  >  486. 5« 

0  7,n 

reject  otherwise. 

2o  If  the  above  procedure  is  used,  find  the  probability  of  accepting  a 
lot  having  mean  life  9  *  500  ,  9  ■  250. 

Solution:  The  result  caul  be  obtained  readily  from  the  0,C.  curve  in 
Figure  1.  Analytically  we  use  formula  (9).  From  this  we  get 


I 


1 
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L(JOO)  «  Pr  (X2<10)  >  2X29  (io>)  ■  Pr(X2<10)  >  9.73oj 

*  .4? 

L(250)  »  Pr^ 2 (10)  >  4^29(10)j  «  Pr  j^2(10)  >  19.46^ 

»  .037. 

3.  If  the  above  procedure  is  used,  for  what  value  of  0  la  L(0),  the 
probability  of  acceptance  =  .50?  *  .10?  =  .05? 

Solution:  Using  table  2  (c),  L(0)  *  „50j  for  0  «  521;  L(0)  «.  .10,  for 

0  a  304;  and  L(9)  *  .05.  for  0  a  266. 

4.  Find  the  expected  waiting  tine  in  the  non-replacement  case  for  the 
following  choices  of  n:  5,  10,  20,  assuming  0  «  Qo  =  1000. 

Solution:  From  formula  5  or  more  easily  fron  Table  3  (a)  or  3  (b)  we  get 


E(X5^5>  «=  2283,  ®<*5#10>  a  ^5o6,  and  ECXj^) 


279.5  . 


5„  Same  as  4  in  the  replacement  case. 

Solution:  E(X_.  _)  ■  rO./a 

■  r,n  o 

and  E(X?  $)  «=  1000,  E(X5  1Q>  *  500  and  E(X^  2Q)  *  250. 

6.  Find  a  truncated  procedure  based  on  X^  in  the  non-replacement 
case,  with  Type  I  error  oC*  .10  . 

Exact  solution:  Compute  p  so  that 


£  (“)  pNi-p)10-11 

k-5  * 


.10  . 


Using  binomial  tables  we  find  p  »  .267.  Thus 

T  *  log  ■  log  ■  .3106. 

Thus  the  test  procedure  is:  Accept  if  x,.  >  311  and  reject  otherwise. 


2ol2 

According  to  (1?)  an  excellent  approximation  is  given  by  lotting 

T  *  Of?  (2r)E(X_  )/2r  *  (486.5)  (.6456)  •  314  . 

o/vi^  rtn 

See  also  table  4  (c)  with  raj  and  n  a  2r  =  10. 

The  two  values:  311  (exact)  and  314  (approximate)  are  very  close. 
7.  Find  a  truncated  procedure  based  on  x^  in  the  replacement  case. 
Solutions  According  to  (10)  the  appropriate  truncation  procedure  is  to 
accept  if 

x_  >  (1000) (4.865 )/20  a  243.2 
and  reject  otherwise. 
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Se  ction  2 

A  censored  toot  having  tho  property  that  its  O.C.  curve  is  such  that 

L(Q  }  ■  1  -  oC  and  U0. )  <  8  , 

O  X  “  ' 

One  frequently  wishes  to  design  a  life  test  which  requires  that  the 
O.Co  curve  meet  the  following  proscribed  conditions: 


(i) 

If 

©  -  ©o, 

then  L(©  ) 
0 

«  1  -  OC 

(ii) 

If 

®  3  9lt 

than  LCGj^) 

% 

VI 

where  9  >  ©,  • 

0  X 

Put  into  words  wo  have  a  situation  where  lots  having  mean  life 
©  >  ©o  are  considered  desirable!  lots  having  mean  life  9  <  ©^  are 
considered  undesirable.  The  interval  (©^,  Oq)  is  essentially  a  zone 
of  intlifferenco.  The  oC  and  /5  can  be  thought  of  as  producer’s  and 
consumer’s  risks  or  as  errora  of  the  first  and  second  kind,  respectively. 

The  problem  amounts  to  choosing  r  in  (7)  in  such  a  way  that  not 
only  is  L(©o)  *  1  -  dC  ,  but  also  that  L(©^)  <  j8  c  From  (7)  and  (9) 
it  is  clear  that  those  two  conditions  are  met  if  r  is  such  that 

(18)  xf  .*(2 r>2  X|  <2r>  or  <  X*  <2r>/ X*  (2r)  . 

1  0 

More  precisely  we  want  the  smallest  r  meeting  this  condition; 
i,e.t  we  want  that  integer  r  which  is  such  that  the  associated  O.C.  curve 
passes  most  nearly  through  the  points  l  ©o  ,  L(©q)  •  1  -  oC  ]  and 
[  ©^  ,  1.(0^)  a  fi  3  „  It  is  readily  verified  that  as  r  goes  through  the 

values  1,  2,  3.  the  ratio  X*  _Q(. (2r)/  XJf  ( 2r )  Increases  monotonically 

to  unity.  Consequently  we  can  always  find  an  integer  r  such  that 

Q 

Xi  (2r)/Xa2  <2r)  >^i->  Xi  (2r-2)/ <*r-2>  • 


(19) 
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f 


{ 


Ttyis  is  the  value  of  r  which  we  want  to  uea.  Ueing  this  value  of 
r„  the  region  of  acceptance 


(20>  V» >  °o  Xf-.(2r)/2r 

is  such  that  its  associated  0UC.  curve  hae  L(0o>  =  l-«<*c.  and  1,(0^)  <  „ 

In  Table  5  we  give  the  appropriate  values  of  r  and  of 

for  the  16  number  pairs  («<  ,  >3)  which  can  be  made  with  the  numbers  (.01, 

~  0 

o050  .-10,  „25)  and  the  values  k  »  2«  3s  5i  and  10. 

U1 

Remark •  It  will  be  noted  that  the  valuos  of  r  required  for  fc  = 
are  quite  large*  It  is  our  feeling  that,  generally  speaking,  it  io  rare 
that  one  would  want  to  work  with  values  of  k  <  — ,  In  case  this  .is  so, 
however 9  we  should  like  to  indicate  how  vre  would  find  the  required  r 
and  ^  ^  (2r)/2rf. 

Since  r  must  be  large  we  use  the  approximation  that  ^2(2r)  is 
distributed  approximately  normally  with  mean  2r  and  standard  deviation 
2'/  r  .  Thus  to  require  that  the  O.C.  curve  bo  such  that 

C 

L(Qo)  =  1-wc  and  K^)  <  with  =  k  <  -f-  , 

we  choose  the  smallest  integer  r  such  that 


XilflL 


k0  This  means,  using  the  normal 


approximation,  finding  r  such  that 


2r  +  2C^j  V  r 
Hr  +  2C  .  V~r~ 

oc 


4 


or  95  «  2.326,  1.645,  1,282,  .674  for  or  jB  »  .01,  .05u 

,10,  „ 25, respectively. 


Men** 
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r 


Once  we  have  found  r, 

2r  •»  2(^^  \Tt~ 
2r 


th«n  ^^(<(2r)/2.r  is  given  by 

,  C<*- 
=  1-  - -  * 


Numerical  Example a 

lo  Design  a  censored  life  test  which  will  meet  the  following  conditions: 

When  9  *  1500,  L(©  )  a  .95  and  when  0,  *  500,  L(0. )  <  ,05. 

Solution:  In  the  problem  k  =  0  /8,  *  3,.  oC  .05,  and  £  *  .05. 

17  1  1  ■ 1  oi  ~ 

Using  Table  5»  wo  eeo  that  r  =  10  and  the  region  of  acceptance  is  given 

by 

V»  >  0  ■  °o  X?95  <20)/2°  *  8l5‘ 

For  this  procedure  1,(0^)  -  «9?  exactly  and  L(C^)  a  ,033  . 

^  2,  Find  the  appropriate  v  for  the  case  where  k  *  0o/®x  -  1*1  and 

«*.  *  /5  =  ,05. 

Solution:  Using  the  remark,  wo  know  that 

r  .  [(2i)a.6wjj  2 

(34. 55)2  «  1194. 


£ 


Truncated  life  teats  having  the  property  that  the  associated 
O.C.  curve  is  such  that  L(®0)  £  l*^  and  L(Q^)  - 

In  tho  previous  section  wo  considered  censox*ed  life  tests,  i.'e.c 
teats  in  which  life  testing  stops  after  a  prescribed  number  of  failures, 
r,  have  occured.  While  such  tests  do  in  general  have  the  desirable 
effect  of  shortening  experiment  time,  there  is  nevertheless  the  ever 
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present  feature  that  one  does  not  knqv  precisely  when  the  experiment 
will  end,  since  this  depends  on  the  random  time  Xr^Q.  As  a  matter 
of  fact,  It  Is  frequently  necessary  because  of  practical  considerations 
to  terminate  a  life  test  by  a  preassigned  time  TQ  ,  a  requirement 
which  censored  tests  do  not  meet. 

If  we  wish  to  truncate  en  experiment  by  a  preassigned  time  TQ 
we  are  led  to  truncated  life  tests  in  which  it  is  decided  In  advance 
that  the  life  test  will  be  terminated  at  mln(Xr^n;  TQ)  where  Xj. 

is  the  time  when  the  rQ’th  failure  occurs  and  TQ  is  the  truncation 
time  beyond  which  the  life  test  will  not  be  allowed  to  run.  (Both  rQ 
and  Tq  are  preassigned. )  If  the  life  test  is  terminated  at  Xj.  n 

(l.e.,  rQ  failures  occur  before  time  TQ),  then  the  action  taken  will 
be  to  reject.  If  the  experiment  is  terminated  at  time  T0  (l.e. ,  the 
rQ'th  failure  occurs  after  TQ),  then  the  action  in  terms  of  hypothesis 
testing  is  acceptance.  It  can  be  shown  that  three  functions  characterize 
the  test  procedures  in  either  the  replacement  or  non-replacement  case. 
These  are: 


(i)  E0(r),  the  expected  number  of  items  falling  before  reaching 

a  decision, 

(ii)  Eq  (T) ,  the  expected  waiting  time  to  reach  a  decision,  and 

(iii)  L(  o)  ,  the  probability  of  accepting,  if  the  true  value  of  the 
mean  life  is  0  . 

In  the  non-replacement  case 


(21) 

where 


Eq (r)  *  np0B(ro  -  2j  n  -  1,  pQ )  +  rQ[l  -  B(r0-  1}  n,p0  )] 


P0  =  1 


-Ve 


-  e 


S-MaHt*!  IT'S** 
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The  probability  distribution  of  r  is  given  by 

Pr(r-k|©)  »  b(k;n,pe),  k  »  0,1,2,..  .,r  -l 


(22!) 

Further,  one  haa 


Pr(r«ro\0)  »  1  -  B(ro-l;n,pe)  . 

r 

o 

vt>  -  Y1  %  <\,d) 


where  Sg( p)  ^an  be  found  froa  '(5),  end 

V1 

(2k)  1(0|  a  ^rsJt|®)  38  ®(VljD/I>S^  * 

fe 0~ 

In  the  replacement  case  the  probability  distribution  of  r  is  given  by 
(25)  Prir^kjs)  n  pt'kjh^),  k  «  0,1,2,  ...  ,r  ~l 

V1 

(25’)  ?ri.r-«roJ0)  «  1  -  p(kj>.d)  ~  1  -  jt'r^-lj.^) 

feo 

In  (25)  and  (35’),  \  *■  ,  pCk;^)  ■  e  ®*£/k* 


■=» 


y~~  p(k;kQ)  . 


Further,  one  lias 


Vr)  •’  V^V2i  V  +  rJ^k(r0-i;^)] 

E  (t)  -*■  9E  (r)/n 


L(9)  a  a(ro-l;^) 
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Ws  have  just  given  formula#  for  the  0,C.  curve,  the  expected  waiting 

time,  and  expected  number  of  items  failed  in  the  course  of  reaching  a 

decision  for  any  preassigned  a,  T  ,  r  .  The  problem  io  to  find  the 

o  o 

appropriate  truncated  test  (i.e.,  to  find  r  and  n)  when  the  trun- 

o 

cation  time  riQ  is  preassigned  and  the  O.C.  curve  is  required  (for 
preassigned  type  I  error  oc  and  typo  II  error  /3  )  to  bo  ouch  that 
L(0o)  >  l-«<  and  L(O^)  <  jg  .  It  can  be  shown  that  for  both  tho 
replacement  and  non-replacement  cases  the  appropriate  is  precisely 

the  one  used  in  the  censored  test  (20)  and  tabulated  in  Table  5-  In 
the  replacement  case,  the  appropriate  value  of  n  one  should  choose 
is  given  by 

(29)  a  «  [°o^  1-*, ( 2ro J /'2To] 
where  [x]  aeansj  the  greatest  integer  <  x. 

In  the  non-roplacoment  situation  a  good  approximate  value  of  av 
in  case  0q/To  i8  substantially  mors  than  one  (say  >3),  Io  given 

by 

r  t  /c  1 

(30)  a.  =  jr^/d-e  °  )J 

where  C  =  ©  (2r)/2r  and  where  r  io  the  same  as  in  the 

oAi”ot>  o  '  o  o 

replacement  casc0 

In  Table  6  (7)  we  give  the  appropriate  volute  of  n  to  use  in 
the  replacement  (non-replacement)  case  when  <X  *-•  .01,  .05,  .10,  .25; 
jd>  *  .01,  .05,  .10,  .25;  ®0/®i  ■  ”2~»  2t  3,  5*  10?  and 

®q/To  »  5,  5,  10,  20.  The  values  of  n  have  been  checked  by  computing 
L(9q)  and  LtO^)?  the  0„CU  curve  of  the  truncated  test  does  come  very 
close  to  meeting  the  requirements  L(0q)  >  1-^i  and  L(©^)  < 
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car. 


it  at  this  point  and  ffcouLd  bo  ■»<]«» 
the  teat  iasert  on 

n  ic 


Several  remarks  appear  to 
Ramark  li  In  the  replacement  ease  the  O.C. 

»in(X  |T  )  ,  where  the  values  of  r  are  given  in  Table  5  an^S 
ro,n  00 

given  by  (29),  is  such  that  L(0q)  >  1-x  ,  but  in  some  cases  it  my  happen 

that  L(01)  nay  be  slightly  >  yS  •  Thia  c&n  be  avoidod  in  either  of  two  ways, 

Ono  way  is  to  give  the  experimenter  the  freedom  to  une,  instead  of  i-he  truncation 

time  ,  the  slightly  larger  truncation  time  »  0qX^_  x(2rQ)/Jn  ,  The  teat 

based  on  min(X  |T>)  will  have  L(Q  )  »  1-<jC  and  L(« ')  <  A  .  The  other. 
o,n  °  °  ~  ' 

way  is  to  use  (n+1)  items  throughout  the  test  and  i*>  aac* ,  instead  of  T  ,  the 

o 

slightly  smaller  truncation  time  T£  ■  ®0^  (2r0)/2(n*l)  ..  The  beet  based  cn 

nd.n(Xr  jT^)  will  have  L(G  )  ■  1-ot  and  L(0,)  <  /3  °  In  moat  cases  it 

will  be  a  matter  of  indifference  which  procedure  one  adopts,, 


Rem*rk  2:  A  good  approximate  solution  for  finding  a  truncated  non-replacement 
test  procedure  was  given  by  (30).  An  alternative,  more  direct  (and  also  more 
lengthy)  procedure  for  finding  a  truncated  non-replacement  teat  meeting  the 
conditions  prescribed  is  to  note  that  such  a  test  is  equivalent  to  a  binomial 
situation  in  which  we  test  pQ  ■  1-e  ^  0  against  »  l-e  ^  ^  and  want 
the  O.C,  curve  to  be  such  that  L(p0)  >  1-ct  and  Wp^)  <  y3  .  Stated  in  the 
language  of  sampling  inspection,  we  are  seeking  a  sample  size  n  and  a  rejection 
number  rQ  such  that  the  resulting  O.C.  curve  has  the  property  that  L’p)  »  l.-©c 
for  lots  with  p  <  p0  and  L(p)  <  y3  for  lota  with  p  >  p^  .  The  detailed 
calculations  necessary  to  determine  n  and  r^  can  be  carried  out  in  any  give- 
situation  by  using  the  Binomial  Tables  or  Tables  of  the  Incomplete  Beta  Function. 

A*8***  1*  Ife  la  appropriate  to  mention  that  truncated  test  procedures  of  the 
kind  considered  in  this  section  are  good  nlee  of  action  in  cases  where  tho 
underlying  life  distribution  is  not  necessarily  exponential.  More  precisely, 
wo  mean  the  following:  Suppose  that  on  acceptable  lot  of  electron  tubes 
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is  :m  fntr  vhich  the  probability  of  failing  betort  n0  ’  ’  ^  ?0  «** 

t&at  -an  unacceptable  lot  la  ooa  for  which  tha  probability  of  failure  before 

uojhw  tins  1'  la  >  p,(p,  >  p  )  and  suppose  we  want  to  O.C.  curwe  to  be  a-  a 

that  L(pq)  >  1-a  and  L(pj)  <  0  .  It  la  clear  that  the  conmsnts  eade  in 

remark  2  are  relevant  here  and  that  the  teat  procedure  involves  finding  a 

jaapl*  site  n  and  rejection  number  rQ  such  that  we  will  accept  the  hypo  hesia 

that  p  *  po  if  the  nuaber  of  defectlvea  (failures  before  Tq)  in  the  samrle 

<  (r  -  1)  and  reject  the  hypothesis  that  p  «  pQ  (accept  p  •  p^)  if  the 

number  of  defectives  in  the  saaple  >  rQ  .  This  test  procedure  clearly  is 

truncated  and  ha*  the  property  that  L(p  )  >  1-a  for  any  distribution  F  (x , 

fi’  °  ~  ° 

which  Is  such  that  dF  (x)  <  pQ  and  L(p^)  <  P  for  any  distribution  F^.x) 

which  i«s  such  that  IF^(x)  >  p^  .  If,  in  particular,  Fq(x)  *  l-e~x^9o  , 

with  Sy  m  T^/log  ^ —  and  F^(x)  »  l-e_x^9l  ,  with  0^  *  TQ/log  y-i — 

the  test  procedure  Just  described  has  the  property  that  L(0q)  >  1-a  and 

1,(9^}  <  P  .  Recalling  tljat  the  rule  of  action  can  be  written  as  accept  if 

■in(Xr  niTQ)  *  PQ  and  reject  if  mln(Xr  n;TQ)  •  *r  #n  »  »  bave  precisely 

the  truncated  procedure  which  one  gets  in  the  exponential  case  when  testing 

0Q  against  with  L(0q)  >  1-a  and  1.(0^)  <  P  .  But  from  the  preceding 

argument  the  test  procedure  is  distribution  free  in  the  sense  that  it  la  the 

appropriate  one  to  use  when  we  wish  to  distinguish  between  two  distributions 

F  (x)  and  F, (x)  with 

O  A 


j  °»0(«)  t  p„ 

o 


and 


o 

«Fj(x)  >  P 


1 


ror  all  Sucu  caaeti  i.(F  )  >  l-°  and  L(F^)  <  P  .  In  this  connection  it 
is  useful  to  point  out  how  Table  8  can  be  used  to  design  Single  saaqple  plans 
which  are  such  that  when  p  -  pQ  ,  the  probability  of  acceptance  equals  1-a 
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ud  when  p  ■  p1  ,  the  probability  of  rejection  is  greater  than  or  equal  to 

1  -  P  •  It  is  assuaed  that  0  <  pQ  <  p  <  .10  and  that  the  producers'  risk 

Of  and  consumers !  risk  P  are  allowed  to  assuae  the  three  values  .01,  .05,  .10, 

thus  giving  rise  to  nine  possible  pairs  (a  ,  p)  .  In  entering  the  table,  the 

ratio  p./p  should  be  identified  with  0  /o.  .  When  this  is  done  the 
ao  Ox  2 

rejection  nuaber  is  given  by  rQ  and  the  sample  site  is  given  by  £)[  (2rQ  >/<| 

where  CxD  wans  the  greatest  integer  <  x  .  In  this  table,  rQ  is  the  upper 
number  and  Xl-O'^po^2  **  the  lover  number. 

The  Justification  for  the  procedure  stems  from  the  fact  that  the  probability 
of  draving  a  defective  item  from  a  lot  having  fraction  defective  p  is,  for 
small  p  ,  essentially  the  ssoe  as  the  probability  that  a  failure  occurs  when 
one  observes  a  Poisson  process  having  failure  rate  p  (per  unit  time)  or  mean 
lift:  1/p  ,  for  a  unit  of  observation  time.  Thus,  draving  a  sample  of  site  N 
from  a  lot  having  fraction  defective  p  can  be  thought  of  bb  draving  one  Item 
at  random  from  a  lot  vhose  items  follow  a  life  distribution  vhich  is  exponential 
vith  mean  life  0  *  l/p  ,  placing  this  item  under  test;  replacing  it,  vhen  it 
fails,  by  a  nev  item  drawn  form  the  lot,  and  terminating  experimentation  either 
as  soon  as  rQ  failures  occur  or  vhen  Tq  *  5  units  of  time  have  elapsed, 
whichever  comes  first.  This  means  placing  n  ■  1,  0q  m  l/pQ  in  formula  (29) 
and  getting  the  sample  site  N  -  T,  r.[f  1-a(2r0)/2poj  . 


( 


numerical  Examples 


X.  Find  a  truncated  replacement  plan,  for  which  Tq  <*  500  hours,  which  will 
accept  a  lot  with  mean  life  «  10,000  hours  at  least  95  P^r  cent  of  the  time 
and  reject  a  lot  with  mean  life  *  2000  hours  at  least  95  per  cent  of  the  time. 
Compute  L(0),Ed(T)  ,  and  Efi(r)  at  0  =.  10,000  and  0  »  2000,  respectively. 
Solution:  In  this  case  ©Q  *  10,000,  ©^  m  2,000,  0L  »  »  .05. 

Since  © o  5  ,  it  follows  from  Table  5  that  rQ  *  5  .  From  Bible  6, 
we  find  that  corresponding  to  ©q/®^  *  5  >  ©q/T^  *  20,  &  »  .05,  n  «  39* 

Thus  the  following  truncated  replacement  plan  meets  the  requirements:  Start 
the  life  test  with  n  a  39  items.  As  soon  as  one  itam  falls,  replace  it  by 
a  new  item.  Accept  the  lot  if  mtn(X_  500)  ™  500  and  reject  the  lot  if 

J»  27 

min(X^  500)  =  If  the  lot  is  rejected,  experimentation  is  stopped  at 

>  the  tioe  of  occurrence  of  the  fifth  failure. 

For  9  ^  10,000,  Xq  -  nT^©  *  (39)(5«>)/lO,0OO  *  1.95»  Using  Molina’s 
Tables,  one  finds  from  (28)  that  L(©)  *  .95©*  Substituting  in  (26)  and  {27), 
respectively,  gives  EQ(r)  »  1.93  and  Bg( T)  «  495.  For  ©  .»  2,000, 

Kq  *  nl^/9  *  (39)(500)/2,000  »  9  75.  For  this  value  of  ©  ,  !(©)  *  .034, 

EG(r)  *  4.95,  and  E^(T)  *  254. 

2.  Same  as  1  except  that  we  want  a  non-replacement  procedure. 

Solution:  r  »  5  .  According  to  Table  7,  the  sample  size  is  n  »  42.  For 
'  —  0  -T/0 

0  m  10,000,  T./8  «  .05,  and  pQ  3  1-e  ■  .049.  Using  the  Binomial  Tables 

one  finds  from  (£4)  that  L(©)  »  .94 6.  Substituting  In  (21 )  and  (23)  gives 
E0(r)  «  2.02  and  2g(T)  »  494  .  For  ©  »  2,000,  T J*  «  .25.  For  this 
value  of  ©  ,  L(©)  -  .031,  Eg(r)  «  4.91  »  and  Eg(T)  «•  248. 

3.  Consider  the  truncated  replacement  plan  meeting  the  conditions  of  Problem  1. 
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For  what  values  of  9  Is  L(e)  »  .5 T  What  are  S?(r)  and  Eg(T)  for  thi* 
value  of  0  ? 

Solution:  To  find  the  «  such  that  L(0 )  »  .5  aeans  finding  such  that 
7R4j?^)  »  .5  •  Using  the  Molina  tables,  ve  see  that  tola  aeans  *  4.67  • 
Since  n  rtsjB  ,  with  n  »  39  ,  TQ  ■  500  ,  one  finds  0  ■  4100  .  Frost  (26) 
and  (87)  ve  find  Bg(r)  -  3=97  end  Bg(T)  »  424  . 

4=  Consider  the  truncated  non-re  placement  plan  neeting  the  conditlona  of 
Problem  2  .  For  vhat  values  of  9  is  L(0 )  ■  =5  T 

Solution ;  Ibis  aeans  finding  p  such  that  B(4j42,p  )  »  .5  „  Using  the 
Binomial  tables  this  aeans  pg  ■  .1104  .  Since  pg  m  1-e  ,  the  appropriate 

9  v  4874  .  In  this  problem  Eg(r)  and  Eg(T)  will  be  approximately  the  same 
as  in  the  replacement  case.  The  exact  calculations  can  be  aide  from  (81)  ,  (28), 
and  (23)  and  are  left  to  the  reader. 

5.  Find  a  life  test  having  the  following  properties:  It  will  accept  at 
leaet  95$  of  the  lots  for  which  the  probability  of  failing  before  tarn  time 
Tq  le  <  .01  and  will  reject  at  leaat  90$  of  the  lots  for  which  the  proba¬ 
bility  of  failing  before  T  is  >  .05  . 

o  — 

Solution:  In  line  with  what  was  said  in  Remark  3  thla  means  finding  tbs 

appropriate  sampling  plan  for  the  special  case  where  pQ  »  .01  ,  p^  .  .05  , 
o  -  .05  ,  0  -  .10  .  in  this  case  Pj/po  -  5,rQ  -  4,  and  from  Table  8,  the 
ea^le  else  R  «  (1.37)/,01  ■  137.  Thus  the  life  test  la  as  follows:  Place 
137  items  on  teat.  If  4  or  more  failures  occur  before  time  Tq  ,  reject,  if 
3  or  fever  failure*  occur  before  time  TQ  ,  accept.  It  should  he  noted  that  we 
are  not  making  any  assumption  about  the  underlying  distribution  of  life. 


( 


Section  4 


Sequential  Lift  Tests 

It  can  be  shown  tbat  sequential  life  teats'  are  superior  to  either  censored 
or  truncated  life  teats.  It  is  shown  in  a  paper  by  Spate in  and  Sobel  that  the 
sequential  probability  ratio  test  of  A. Wald  can  be  applied  to  life  testing. 

The  Interesting  point  now  is  tbat  decisions  can  be  made  continuously  in  tine. 

At  each  moment  t  ,  one  can  decide  either  to  accept,  to  reject,  or  to  continue 
the  life  test.  If  ve  are,  as  before,  testing  B Q:0  *  0Q  against  ^  :  0  « 
m  e1  (0q  >  0^)  with  Type  I  error  «  c  and  Type  II  error  »  B  ,  then  the  decision 
as  time  unfolds  depends  on 


where  A  and  B  can  for  all  practical  purposes  be  taken  a? 

(38)  A  *  (l-P)/a  and  B  *  P/(l-a) 

In  (31)  r  is  the  number  of  failures  observed  by  time  t  and  V(t)  is 
a  statistic  which  equals  the  total  number  of  houra  lived  by  all  items,  failed 
and  unfailed,  up  to  tl'.ne  t  .  In  the  replacement  case 

(33)  V(t j  -  nt  , 


while  in  the  ncn-replacemect  case 


(3^ 


«  (n-i+l)  (x  -  x  .)  ❖  (nr)  (t~x  ) 

i-1  1  r 

«  £  x.  ♦  (n-r)t  „ 


In  W  \ 


denotes  the  time  of  the  i*th  failure. 
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R emark:  The  decision  to  continue  exper imentat ion  is  made  as  long  as  the  inequality 
(51)  holds.  As  soon  as  (31)  is  violated,  one  accepts  1^(1. e;0  s  0o)  it  the 

function  of  t  in  (31)  ls<B  ,  and  one  rejects  IQ  (accepts  %)  if  the 

function  of  t  in  (31)  is  >A. 

Remark:  It  should  be  noted  that  in  the  non-replacement  case  a  special  problem 
arises  if  all  n  items  fall  vlthout  reaching  a  decision.  This  eventuality  can 
be  taken  care  of  in  various  ways  and  will  be  discussed  in  the  Appendix. 

If  we  wish  to  graph  the  life  test  data  continuously  in  time,  it  is  convenient 
to  write  (31)  as 

(35)  -h]_  +  rs  <  V(t)  <  hQ  +  rs 

where  hg  ,  ,  and  s  are  positive  constants  given  by 


A  good  way  to  describe  hQ  ,  ,  and  s  ,  is  as  follows: 

hQ  is  the  intercept  on  the  total  life  axis  of  the  accept  line;  -h^  is  the 
intercept  on  the  total  life  axis  of  the  reject  line;  and  s  is  the  common  slope 
of  the  two  straight  lines. 

The  O.C.  curve,  i.e. ,  the  probability  of  accepting  HQ  when  ©  is  the 
true  parameter  value,  is  given  approximately  by  the  parametric  equations 


by  letting  the  parameter  h  run  through  all  real  values.  The  values  of  L(  ©) 
at  the  five  points  ©■  O,01,s,©o,eo  enables  one  to  sketch  the  entire 
curve.  These  values  are, respectively,  0  ,  /3  ,  log  A/(log  A  -  log  B)  ,  1  -a  , 
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w 


ud  1.  lot#  that  In  rlew  of  (36),  L(s)  •  log  A/(log  A  -  log  B)  -  h1/(h.  ♦hj),  An 
approximate  formula  for  Z^(r)  ,  the  expected  number  of  items  to  reach  a  decision, 
when  6  is  the  mean  life  is  given  by 


K0)  KttB.fr-  Uofl  log 

10*<8A>  -  «<?,  -  r  > 


*  hrL<e)<Vbi>  n , 

“  -  - - -  9  f  « 


(38) 


E9(r) 


h_  h. 


—  log  A  log  B  o  1 

-  — x—  ,  0  «  a  • 
flog  ejej  a 


If  we  let  k  ■  9q/0^  ,  the  approximate  values  of  E0(r)  become  particularly 

simple  when  0  »  0,  ,  e  ,  or  9 
j.  o 

They  are 


(39) 


E0^(r)  [P  log  B  +  (l-P)  log  aJ  j  [log  k-(k-l)/kj 
E0(r) ~  -  log  A  log  B/(log  k)2 

Eg  ~  [(l-°)  log  B  ♦  a  log  a]  j  [log  k-(k-l)]  . 


It  can  further  be  shown  that  Ea(V(t))  the  expected  amount  of  total  life 
observed  in  reaching  a  decision  is  connected  with  g  (r)  by  the  identity. 

U 

(40)  E0(V(t))  -  9E0(r) 


in  either  the  replacement  or  non-replacement  case.  Since,  in  the  replacement 
case  V(t)  ■  nt  ,  it  follows  that  Eg  (t)  ,  the  expected  waiting  time  to  reach 
a  decision,  is  related  to  Kg(r)  by  the  formulae 


(41) 


Eg(t)  .  0Eg(r)/n 
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Id  the  non-replacement  case, 


(W)  E9(t.).gPr(r-k|e)E9(Xk>n) 


where  Q)  Is  given  by  (5)  .  A  good  approximation  for  2g(t)  is  given  by 

(43)  E„(t)~  9  log^jy)  . 


Practical  Applications 


It  will  be  convenient  to  -normalise  the  preceding  situation  in  6uch  a  way 

that  9  m  1  .  If  this  is  done  it  is  convenient  to  calculate  once  and  for  all 
o 
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the  values  of  hQ  ,  h^  ,  and  s  for  the  cases  where  0,  “  3  '  2  *  3  >  aod 

for  a  >•  .01  ,  .05,  and  P  *  .01,  .05.  In  Table  9  ,  we  give  the  values  of 

h  ,  h,  ,  and  c  for  each  of  these  cases, 
o  1 

In  the  event  that  9q  is  not  equal  to  one  and  that  k  *  0q/b^  *  |  ,  2,  3, 
one  can  readily  find  the  appropriate  equation  for  V(t)  by  multiplying  hQ  ,  h^  , 
and  s  by  0  .  In  Table  10,  we  give  approximate  values  of  E  fr)  for  the 

values  of  a  ,  p  ,  and  k  «  given  above. 


Numerical  Examples 

1.  Find  a  sequential  replacement  procedure  which  will  accept  a  let  with  mean 
life  SQ  *  1500  hours,  95jJ  of  the  time  and  will  reject  a  lot  with  mean  life 
9^  »  500  hours,  95$  of  the  time.  The  constant  number  of  items  under  test 
is  n  m  SO  .  In  this  case,  9q  *  1500  ,  0 ,  *.  500  ,  a  *  0  «  .05  . 

Solution;  Substituting  in  formula  (31)  we  get 


„-t/37.5 


<  19  • 
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In  this  ease  (35)  become* 


-  no  +  klx  <  t  <  HO  +  4lr  , 

vaere  t  represents  the  length  of  tire  that  the  life  test  ha*  been  in  progreoa 
and  r  denotes  the  number  of  failures  obtained  up  to  time  t  .  Itoe  experiment 
is  continued  as  long  as  the  inequality  holds  and  is  stopped  as  soon  as  the  in¬ 
equality  does  not  hold.  If,  at  the  time  of  stopping,  t  is  less  than  the  left» 
hand  member  of  the  inequality,  we  reject  0Q  ■  1500  (accept  0^  •  500)  j  if,  at 
the  time  of  stopping,  t  is  greater  than  the  right-hand  member  of  the  inequality , 
we  accept  0Q  *  1500. 

2.  Compute  Eg(r)  and  EQ(t)  for  0  *  0  ,  0^(»  5°0)  ,  e(*»  8®3)  »  0O{**  1500)  , 
and  on 

Solution:  From  Table  10  we  get  E^(r)  »  3  *  Ea  (*■)  *  6.14  ,  E_(r)  <*  7.18  , 

. . .  o  o 

Eq  (r)  *  2.94  ,  and  E^Cr)  *  0  . 

o 

In  the  replacement  case  E0(t)  is  found  most  easily  for  all  values  of  e(j  00) 
by  using  (41)  ,  Ea(t)  ■  0Efl(r)/n  .  This  gives  E  (t )  »  0  ,  E  (t]  -  155  , 

u  U  O 

E  (t)  «*  295  ,  and  E  (t)  •  220  .  For  0  ■  oo  ,  the  expected  waiting  time  is 
8  0 

,  ^  ^  -t/37.5  1 

given  by  tffi  ,  where  e  m 

This  gives  t  m  E  (t)  -  HO  . 

0  on  cc* 

Remark ;  More  generally,  in  terms  of  B  ,  n  ,  0q  ,  and  k  we  find 

t^  ■  0Q  log  B/n(k-l)  . 

This  means  that  if  no  items  fail  by  t^  ,  w  stop  experimentation  at  tw  with 
acceptance  of  HQ  . 

3.  Assume  that  we  are  testing  the  hypothesis  in  problem  1.  A  sample  of  slse 
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20  is  placed  on  tent.  Items  which  fail  aw  replaced  at  once  by  new  items  drawn 
from  the  same  lot.  The  experiment  is  started  at  time  t  »  0  and  the  first  five 
failures  occur  at  «  20.1  hours,  xg  *  100.5  nours,  x^  *  121.7  hours, 
x ^  m  167*4  hours,  and  «  179*2  hours,  all  times  being  measured  from  t  «  0 

(a)  Verify  that  no  decision  has  been  reached  by  tic*,  x,.  - 

(b)  Verify  that  if  the  sixth  failure  has  not  yet  occurred  at  315  hours, 
measured  from  t  «  0  ,  we  can  stop  experimentation  at  that  time  with  acceptance 
of  Hq  ,  namely  that  9q  *  1500  . 

Solution :  VFe  remarked  in  the  poiutiors  to  '(l)  that  (35}  becomes 

-  no  +  Ur  <  t  <  no  +  4lr  . 

This  region  is  drawn  in  Figure  (2).  Toe  life  test  data  are  plotted  by  moving 
vertically  so  long  as  we  are  waiting  for  the  next  failure  to  occur.,  and  moving 
horizontally  by  one  unit  (in  r)  at  each  failure  time.  In  Figure  (2)  the 
path  crosses  into  the  region  of  acceptance  when  r  «  5  ,  at  time  t  «  llo  •> 

♦  (fcl)  (5)  *  315  .  Siace  the  sixth  failure  has  not  yet  occurred  we  can  stop 
life  testing  t  «  315  hours,  with  acceptance  of  Hq  . 

Remark;  As  a  matter  of  fact, we  happen  to  know  in  this  example  that  the  sixth 
failure  occurs  at  x ^  *  346.7  hours.  Thus,  as  indicated  in  Figure  (2)  ,  we 
saved  346.7  -  315  «  3d *7  hours  by  virtue  of  the  "act  that  life  test  data 
were  becoming  available  continuously  in  time. 

4.  The  first  six  failure  times  in  a  sample  of  20  (with  replacement)  ere 

x^  *  19*3  *  *2  m  "*5*®  ,*  *  49.9  ,  *  96*7  ,  *  115.2  ,  x^  •  127.7  • 

Verify  that  if  the  hypotheses  being  tested  ere  those  in  Problem  1,  then  fiQ 
is  rejected  at  time  x^  «  127*7  hours. 

Solution;  x^  ,  xg  ,  ,  and  x^  all  fall  within  the  region  bounded  by 
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the  two  straight  Unas.  However* when  r  •  6  ,  -  110  +  4lr  »  136.  Since 
Xg  «  I27.7  <  136  ,  Hq  is  rejected  at  iiime  Xg  ■  127-7  hours.  A  graphical 
solution  is  given  in  Figure  (3)  . 

Remark;  While  the  acceptance  of  Hq  (9q  ■  1500)  in  problem  3  is  made  between 
failure  times  and  Xg  ,  rejection  of  HQ  in  Problem  4  is  made  at  the 
failure  time  Xg  ,  with  an  excess  over  the  boundary.  ‘Phis  illustrates  the  point 
that  acceptance  of  Hq  is  always  made  between  failure  times,  whereas  rejection 
of  Ho  is  always  made  at  a  failure  time. 

5"  Find  a  truncated  (nonsequential)  replacement  procedure  for  testing  the 
hypothesis  in  Problem  1,  using  a  constant  sample  site  n  *  20. 

Solution ;  From  cur  earlier  results  dealing  with  truncated  replacement  procedures 
it  can  readily  be  verified  that  the  truncated  replacement  procedure  meeting  the 
requirements  is 

(i)  Tf  min  [ f  407- 5J  *  407 • 5  ,  truncate  the  experiment  at  407. 5 
with  acceptance  of  Hq  . 

(ii)  If  min  [x? Q  .  407-53  *  xjo  *  tru®cate  the  experiment  at  x1<?  with 
acceptance  of  Hj 

The  O.C.  curves  of  this  test  procedure  and  of  the  one  in  Problem  1  are 
essentially  the  same. 

6.  Compute  E0(r  j  and  E0(t)  for  the  plan  in  Problem  5  for  9  *  0,9^1,9^00 . 

Solution;  Using  the  formulae  given  in  the  section  on  truncated  replacement 
procedures  and  recalling  that  E  (t)  -  (r)/n  ,  one  gets  E.(r)  »  10  ,  9.93  , 

u  u  u 

fi>75  ,  5-39  *  0  and 

l0(t)  .  0  ,  248  ,  36O  ,  404.5  ,  407.5  for 

8  «  0  ,  0^500)  ,  0(823)  >  9o(1500)  f  and  ,  respectively. 


■  ! 

*  » 

1 
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7.  Compare  Eg(r)  and  E0(t)  for  the  test  procedure*  obtained  as  eolutlons 
to  Problem*  l  and  5  » 


Solution :  Using  the  solutions  to  Probleas  2  and  6  ,  one  has  the  following 

comparisons; 

B„(r> 


Truncated  with  replacement  rule 

Em 

0  «*  500- 

0  «  823 

9  m  1500 

0  m  oo 

10 

9.93 

8.75 

5.39 

0 

Sequential  rule 

3 

6.14 

7.18 

2.94 

0 

S9(t) 


'Truncated  with  replacement  rule 

9  m  0 

0  *  500 

0  M  823 

9  »  1500 

0  *  00 

0 

248 

360 

404.5 

407.5 

Sequential  rule 

0 

155 

295 

220 

no 

These  tables  give  e  fairly  good  idea  of  the  savings  associated  with  adopting  a 
continuous  sequential  rather  than  a  truncated  plan  and  are  typical  of  what  may  be 
expected  to  happen..  A  graphical  comparison  of  the  two  procedures  is  given  in 
Figure  4  «, 

8.  Find  t  in  Froblem  1  if  a  «  0  ■  .01  . 
oo 

Solutions  t„  -  -  e_  log  B/n(k-l)  «  230  .  T*1®  u  *bout  twic®  the  value  of 

CO  O 

t  when  a  «  p  «  .05  . 
oo 

9»  Find  L(s)  for  all  nine  conbirtt'.or.s  of  a.  .01  ,  .05  ,  .10  and  0  »  .01,  .05  , 

.10 

Solution;  Since  L(e)  -  \K*Q  +  it  followe  from  Table  9  that  L(s)  ie 
given  by  the  valuee  in  the  following  Table. 
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L(a) 


.01 

.05 

.10 

.01 

.500 

.604 

.662 

•05 

*396 

.500 

.562 

.10 

.338 

.438 

.500 

ftsmarh  ■  Since  LtO)  *  0  ,  1,(6 J  -  1-ct  ,  it©  )  «.  ft  ,  or.a  L(«*  )  ■■ 
readily  drnvf  tbs  C.C.  curve©  lor  ail  cf  the  case©  treated 

3.0  Fiii  a  sequential  test  fo>  the  case  when  «  -  , 0  t-  .»"5  , 
and  p  r  3.00  . 

Solution  i  From  TaKie  5,-  oae  1‘indr  that  b_  »  h,  •*  1.4722  (sine* 
SQ  is  corsialized  ac  1')  .  Therefore  the  region  (35)  i  :  given  by 

;-00( -1.4722  *■  «5493r)  <  V(t)  <  300(1.4722  •*  .5493*“'  . 


I 


d 


we  can 


300  , 


After  simplify!!**  this  becnsee 


-  4*2  4  l65r  <  V(t)  <  442  *  l65r 


The  life  tent  it?  cor.tinved  r,o  long  as  V(fc’;  satisfies  V>th  inequalities. 
As  soon  as  the  iueqvul Ifcice  are  vinluted,  one  accepts  H  (i.e-*  8  *  300)  if 
V(t)  >  442  a  .165 n  and  reject**  Hq  (accepts  84  (i.e..  9j  *  300))  if  V(t)  < 
*-  44 2  +  l65r. 


Host  of  the  results  in  sections  1  and  2  of  Chapter  II  are  proved 
in  the  following  reference: 

Be  Epstein  and  M„  Sobol,  "Life  Testing".  Journal  of  the  American  Statisti¬ 
cal  Association  4£,  486-502,  1953c 

In  Appendix  1  of  that  paper  it  is  shown  that  3  as  given  by  (2) 

T  fU 

is  a  "best"  estimate  of  0  in  the  non-replacement  case  end  the  p.d.f. 

(4)  is  derived.  The  expected  waiting  time  formula  (5)  for  E(X  )  is 

r  ,n 

derived  in  Appendix  2  of  the  reference.  The  "best"  test  based  on  the 
first  r  out  of  n  failures  having  the  prescribed  properties  that  its 
O.Co  curve  is  such  that  L(Oo)  3  1~°C  and  L(Q^)  <  ^  is  obtained 
directly  from  the  Neynan-Pearson  lemma  in  Appendix  3  of  the  reference. 

There  is  little  point  in  writing  down  detailed  proofs  when  they  are 
readily  available  to  the  interested  reader  in  the  reference  just  cited. 
There  is,  however,  good  reason  to  give  some  supplementary  material  which 
is  very  helpful  in  understanding  the  various  results.  This  we  shall  do 
in  what  follows. 

In  life  teoting  problems  where  one  makes  the  assumption  that  the 
underlying  distribution  is  exponential,  the  following  results  play  a 
fundamental  role: 

(1)  Given  a  Poisson  process  for  which  the  rate  at  which  events 
occur  per  unit  time  is  X  »  Let  the  random  function  X(t)  (  X(0)  is 
assumed  equal  to  zero}  be  the  number  of  events  occuriug  in  (0,t)  „ 

Then 

(i)  Pr(X(t)  -  k)  -  s~At(Xt)k/kl*  k  -  0,1,2,. op  o 

Mors  generally,  if  >  t^ ,  then 
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P 


-Mt-.fc,)  . 

(ii)  PrMtgJ-Ktj)  »  k)  «  »  21  W.(t2-t1)]Vkl, 

k  ®  0  { 1  ^  2  (i  c  y . 

(ill)  Let  the  random  variable  T  be  the  waiting  time  until  the 
first  count,  or,  more  generally, the  waiting  time  botaeon  successive  counts i 
then  the  p.d.fo  f(t)  of  T  is  given  by 

f(t)  »  *0“**,  t  >  0 
a  0,  elsewhere 

and  the  c,d.f.  F(t)  is  given  by 

F(t)  a  Oct  <  0 

_  t  >  Or. 

In  Feller's  book  on  Probability  Theory  the  fundamental  postulates 
for  the  Poisson  process  are  given  as  follows:  whatever  the  number  of 
changes  during  (Ovt)  the  (conditional)  probability  that  during  (t„t+h) 
a  change  occurs  is  Xhfo(h).  and  the  probability  that  more  than  one 
change  occurs  is  of  smaller  magnitude  than  hu 

Proof:  (i)  and  (ii)  ar®  direct  consequences  of  the  definitions  of  a 
Poisson  proceed.  To  prove  (iii)  let  T  be  the  random  variable  repre¬ 
senting  the  waiting  ti.ae  until  the  first  count  occurs  (measuring  time 
from  the  origin  cf  tiw?  t  =  0)  or  the  waiting  time  between  two  succes¬ 
sive  counts  (wtaero  ae  vould  now  measure  time  from  the  moment  when  the 
last  count  was  recorded  and  would  wait  for  the  nest  count),  then 

(2Aol)  Pr(T  >  t)  =  Pr  (0  counts  occur  in  an  interval  of  length  t) 

=  Prob  (X(t)  a  0)  a  e*Atc 

Therefore,, 

(2A,2)  F(t)  =  Pr(T  <  t)  a  l-Pr(T  >  t)  *  t~>  0 


«  0,  elsewhere 


t  >  0 
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and 

(2A.3)  f(t)  *  F'(t)  a 

a  08  else whore, 

(2)  A  Poisson  process  baa  tho  following  interesting  feature  by  definition0 

If  0  <  <  t.,  <  c  • «  <  t^  <  <.  #  o ,  then  the  random  variables 

^X(t.)c  Xtt^)  -  X(t^){,,„B  X( -  X(t^  ^)j  &ra  Mutually  independent*, 

(3)  Consider  the  random  variable  T  distributed  with  p»d»fu 

f(t)  b  X*  ^  *  t  Oc  Then  Pr(T  >  t+tjl  >  t)  =  e  ^  ,  Put  into  words: 

If  one  observes  a  Poisson  process  for  a  length  of  time  t  and  no  events 

occur,,  then  the  probability  of  no  overts  occurring  in  an  additional  amount 

^  v> 

of  time  Z  is  given  by  e  ,  This  it  r  special  case  of  (2)  where  one 
considers  only  the  two  intervals  (06t)  and  (tvt+ lr) 

(4)  Putting  (3)  into  life  testing  language  we  haver  given  that  an  item 
has  lived  for  a  length  of  time  t»  then  the  conditional  probability  of 
surviving  an  addition .0.  T  time  units  is  given  by  e~  „  But  this  is, 

of  course t  tho  probability  of  cn  iter  surviving  X  units  ob  initio.  Thus 
if  the  underlying  distribution  is  exponential,  items  that  have  survived 
up  to  any  given  time  arc  "as  gcod  e.3  new"  and  "have  not  aged"0  The  proof 
is  very  easy: 

(2A„4)  Pr(T  >  tt-X  |t  >  t)  *  Pr(T  >  t+*t)/Pr(T  >  t). 


This  was  the  result  that  we  wanfcid  to  prove, 

(5)  If  a  ifceno  each  having  the  p.d,f,  of  life  f(t)  =  are 

placed  on  test  nimultanoously  *?':  some  tirw  t  =  0f  than  the  first  failure 
to  occur,  t^v  is  exponentially  distributed  with  rate  n^  <,  Two  short 
proofs  follow. 


Proof  1;  Recalling  the  connection  between  Poisson  processes  and  the 


exponential  p.d.f. ,  we  can  Imagine  that  we  superimpose  n  Poisson 
processes  each  having  failure  rate  X  .  The  result  Is  a  Poisson  process 
with  failure  rate  n  Hence  from  (l)(lll)  the  p.d.f.  of  the  waiting 
time  for  the  first  failure  Is  given  by  f^(t)  •  nXe”nXt,  t  >0  and  the 
c.d.f.  Is  F1(t)  ■  1-e  t  >  0. 

Proof  2:  An  alternative  proof  is  to  recall  that  the  first  future  is 
the  smallest  In  a  sample  of  size  n  drawn  from  an  exponential  p.d.f. 

The  c.d.f.  of  the  smallest  value  Is  given  by 

F^t)  ■  l-(l-F(t))n  ■  l-e-11^*,  t  >  0.  This  is  the  result  obtained  before. 

We  now  use  these  results  to  obtain  formulae  (2)  through  (6)  In 
Section  1.  First  we  note  that  given  a  Poisson  process  with  rate  X , 
then  the  associated  waiting  time  random  variable  T  has  expectation 


(2A.5)  E(T)  *  f"  Xte“Xtdt=^. 

In  the  life  testing  context,  where  T  Is  thought  of  as  life,  E(T)  ■ 
mean  life  =  0  ■  ^  .  Thus  items  exponentially  distributed  with  mean  life 
0  can  be  thought  of  as  waiting  times  between  successive  occurrences  of 
a  Poisson  process  with  rate  A  :  i  . 

We  now  obtain  formula  (2)  (non-replacement  case)  using  the  ideas  of 
Poisson  processes.  Placing  n  items  on  test  at  time  t  =  0,  with  each 
Item  having  an  exponential  p.d.f.  with  mean  life  0,  Is  equivalent  to 
considering  the  superposition  of  n  Poisson  processes,  each  having 
rate  X  ■  ^  .  The  process  obtained  by  superposition  is  still  Poisson 

with  parameter  A  .  n  .  The  first  failure  observed  at  time  x,  _ 

n  0  £  1»n 
is  exponentially  distributed  with  mean  life  ^  and  so  nx-^  is 

exponentially  distributed  with  mean  life  0 .  Consider  now  what  happens 
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after  time 


*  At  tine  x. 

.fU  X) 


on*  has  (a-1)  items  loft  each  with 


mean  life  ©  (this  is  a  consequence  of  3  and  4)u  Thus  one  is  now  deal¬ 
ing  with  a  superposition  of  (n-1)  Poisson  processes  each  having  rate 

^ ,  and  hence  the  superposition  is  a  Poisson  process  with  rate  (^j~)  « 

© 

Therefore,  x ^  a«x^  Q  is  exponentially  distributed  with  mean  life 

and  (n-l)(x_  -x,  )  is  exponentially  distributed  with  mean  life  9U 

c„n  i9n 

Also,  X,  and  (x-  -x,  )  are  mutually  independent „  Continuing  in 

itn  ^en  X;,n 


the  same  way 


r 


^l.n5  Ji2,n“*l,n  xi(n“xi-l,n***#*  *r,n"*r-l,i 


are  mutually  independent  and  drawn  from  exponential  distributions  having 


mean  lives 


©  J5_ 

n*  n-l‘ 


Q  _ 0 _ 

J>n-i®  n-iv-l 


\ 


respectively.  Jfore 


simply  the  r  random  variables 


y.  c  t  n-i+lXx^-x^^)*  i  -  lr>2,,  „z| 


1  ~y/9 

where  x  cl  are  mutually  independent  with  common  p.,  d„  f*  ~e  y  >  0C 


TheroforQj 


(2Aq6) 


i=l 


T</-  =  £_ 

i=l 


(n-i+l}(x  >-x, 

X  X* 


}/r 


(n-r)x 


r,n 


is  unbiasodo  The  other  properties  (such  as  maximum  likehood,  unbiasedness, 
minimum  variance^  efficiency,  and  sufficiency)  are  proved  in  the  reference 
cited  above.  The  p„d.:?»  (h)  of  follows  directly  from  the  fact  that 
the  cun  of  independent  random  variables  each  of  which  is  exponential 
follows  a  Type  III  distribution.  More  precisely  consider  the  random 

r 

variable  U  «  ?■  y, »  This  can  be  considered  as  the  waiting  time  for 

i^T  1 

the  rth  event  in  a  Poisson  process  with  parameter  X  =  «  The  pid.fc 


2.38 


of  U  is  found  by  using  the  fundamental  postulates  for  a  Poisson  process,. 
Thus 

(2A.?>  Pr(t  <  U  <  t+4*it)  *  Pr(r-1  events  occur  in  (Oft)  and  1 

event  occurs  in  t„t4^t) 

*  PrCr~l  events  occur  in  (Oet)j0  PrCl  event  occurs  in  (t,t+it)] 


r-l  -  | 

6 


£1 1 
© 


TZ&T 

Therefore, the  podcf.  cf  U  is  given  by 


(2Ao8) 


h(t)  *  Tr^Ti  o  (o)  e*t/C»  t  >  °* 


But.  $  *  U/r  and  using  simple  transformations  the  p0duf^  of  d 

r,n  r.,n 

becomes  (4)e  To  prove  (5)->  we  note  that 


(2A<>9)  (X_  -X,  )  *  o»c+  (X  ~X  ,  ) 

/g&  Xg&  2  git  1  t*n  ^ 

Y  I  Y 

=  Zi  +  JA  +  +  g- 

u  +  n-1  +  •••  +  n-r+1 

1  -x/O 

But  the  are  each  distributed  with  the  p„4.f„  ^  e  o  There  i  ore*, 

.  r 

C2A.10)  E(X  >  «  9  (*  +  r*T  +  •••♦  TT77)B  0  HI  1/a-J+l 
rtn  n  n~l  n-r+1  — 

and  thus  (5)  holds 0 

Incidentally «  since  T^’s  ore  also  mutually  independent  it  follows 

that 

,  v  P 

(2A.11)  ParU  )  =  0  I _ l/<n-J*l)  * 

r«n  7a 

also 

(2A.12)  covCX  s  X.  )  =  var  X„  =  02  jfl  l/(n-j+l)2,  if  s  >  r 
r3n  efn  r,n  ^ 


r.n 


c  var  X  *  ©' 

8 


l/(n»j+l)  *  if  r  >  s. 


2a  39 


For  example,  suppose  s  >  r0  Then 
(2A.X3)  X 


.  »  X  +  (X  -  X  ), 

s(n  r(n  s,n  rtn 


Btaoe 

(2A.14) 


coy  (X  ;  X  )  *  COY  (X  ;  X  +  (X 

r,n*  a,n  v  r,n*  r,n  s,a 

“  C0T  (Xr  n*  Xr.  J  +  00T  <Xr  n}  Xe  n 

rfn  r6n  r?n  ■  *& 


x_  _» 

r,n 


Xr 
r  t>n 


Noting  that  X  (the  wait lag  tine  for  the  r 
r  ,n 


th 


Xa,n  “  Xr,n 


th 


failure)  and 
th 


and  b 


failure)  are 


(the  waiting  time  between  the  r 
independent,  we  get  that 

coy  U  :  X.  )  -  Var  (X„  )  if  a  >  r  „ 

r  |D  5||&  P 

Similarly 

coy  (X  :  X.  )  =  Var  (X  )  if  r  >  So 

r,n  s,n  ecn  *• 

Thus  (2Ad2)  i 3  proYedo 

Op  to  thie  point  we  have  doalt  exclusively  with  the  non-replacement 
case.  If  itena  are  replaced  as  they  fail, then  it  is  clear  that  placing 
n  items  on  teat  and  replacing  failed  items  at  once  by  new  items  ie 
equivalent  to  observing  a  Poisson  process  with  rate  A  b  n/©„  If 
x,_<x__<...<x  are  the  first  r  failure  times  (time  being 

Xfii  *  C|Q  •*  •  r | n 

measured  from  the  beginning  of  the  experiment),  then 

/x,  :x_  -x,  .  ...,x _ -x„  ,  l  ere  mutually  independent  and  idea¬ 
ls  A,n  *,n  i,n  r,n  r-i,nj 

cally  distributed  with  common  p.d.f.  The  random  variables 

^s^  *  a  (Xj-x^),  1  *  1,2, ...,rj  ,  where  xq  *  0,  are  mutually 
independent  with  common  p.dof.  x  >  0*  Therefore, 


(2A.15) 


£  ■  III  *j«/r  *  J * 

fcl  1  r'* 


is  unbiasedo  Other  optimum  properties  are  easy  to  show  and  as  before 
the  p.d.f.  of  $  is  given  by  (4). 


To  proTo  (6)  note  that 


2.40 

fn  “  *  *°*  *  (Xr,n  “  X*-l8aJ 

O 

•d  with  p.cUf.  ^  #-x/®o 

£2 

n  ° 

he  Zi,a  it  also  follows  that 

war  X  ■  if  a  >  r 

jTflU  ' 

A 

war  X  «  •5§—  *  if  r  >  s. 

o.n  Z 
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We  hart  seen  that  in  Ufa  testa  where  items  that  fail  are  not 

replaced, then  the  statistic  «  where 

r  9  a 


(2B«1) 


♦«.. +x_  ,  +  (n-r+l)x, 


$  — itiiL — _ £~)>  a.P. 


r  .n 


r,n  r 

la  a  "best"  estimate*  It  can  further  be  shown  that  the  "best"  test  for 

0*0  against  alternatives  ©  <  ©  is  given  by  an  acceptance  region 
o  o 

a  >  Co  If  the  Type  I  error  is  controlled  at  *><  t  then  C  *  Oq)£  **  r.(2r)/2r<- 

Inspection  of  this  statistic  reveals  that  if  n  >  r,  xr  n  is  weighted  more 

heavily  than  the  earlier  observations  x,  „x_  r..,x  ,  ...»  This  would 

x«a  r^AnD 

lead  one  to  suspect  thc.t  estimates  based  only  on  x  (i„8»(  the  r 

a  jU 

failure  time  only)  may  be  highly  efficient  when  compared  with  estimates 

based  on  $  and  further  that  rules  of  action  based  on  x  _  have  0oC„ 
r,n  r,n 

curves  very  close  to  those  based  on  £ _ .  This  question  has  been  stud-> 

r,n 

led  In  detail  in  the  estimation  case  in  a  report  by  B»  Epstein  entitled 
"Estimates  of  mean  life  based  on  the  smallest  value  in  a  sample  of 

sins  n  drawn  from  an  exponential  distribution"*  Wayne  University  Tech¬ 
nical  Report  No*  2,  July,  1952*  It  is  shown  in  this  report  that  a  highly 

efficient  estimate  of  C  is  given  by  /?  x  where 
.  A-/r»n  r,n 


&.n  *  */£  A.«*r.n  “  “ 

It  can  be  verified  readily  that 

(2B.2)  »«  0  *r,„>  ■  V/ 


unbiased  estimate  of  ©c 


where 

(2Bo3) 


E  ,  =  7^’  l/(n-J+D2 

r*a  th: 


/  j-1 


Efficiency  of  the  estimator 


n*r„n 


1/n-j+ir 

relative  to 


r.n 


is  given  by 
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(2B„4)  „  -  Var  “0  /tar  (  >4  x  )  «  — i— 

r#n  r fff  V/r,n  r,n  rKr  n 

la  the  report  to  which  we  just  referred,  tables  are  given  for  ,K 

r,n  r,n 

•®*  ®r,n  *or  n  *  1(1)20(5)50(10)100  and  r  «  l(l)n.  An  inspection  of 

these  tables  rereale  that  E>  „9  for  ^  <  t.  Furthermore  it  ia  shown 

r,n  —  n  5 

that 


s 

r.n 

>  .9990 

if 

r 

—  a 

n 

,1 

w 

>  .9960 

if 

»  a 

.2 

ft 

>  .9893 

if 

»  3 

°3 

•1 

>  .9784 

if 

•»  e 

.4 

n 

>  ,9608 

if 

"  3 

.5 

it 

>  «9329 

if 

'*  a 

.6 

it 

>  .8874 

if 

T»  _ 

.7 

it 

>  .8094 

if 

"  3 

08 

ti 

>  ,6548 

if 

"  . 

.9 

It  has  similarly  been  observed  that  the  O.C,  curvo  resulting  from  ucing 
the  aooeptanee  region  a*r  tt  >  C  virtually  coincides  with  the  0„C, 
curve  assoeiated  with  'd  >  C, 
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A  detailed  discussion  of  truncated  replacement  and  non-replacement 
testa  ie  given  in  the  following  paper: 

Bo  Epetein(  "Truncated  Life  Tests  in  the  Exponential  Case",  Annals  of 
Mathematical  Statistics  2£,  555-564,  1954. 

Section  2  of  thin  paper,  pp„  555-558p gives  proofs  of  formulae  (21) 
through  (28)  inclusive <, 


Appendix  2  D 


Tests  of  the  form  $  >  C  considered  ae  truncated  tests. 

*  «p _ _ _ _ _ 


The  following  material  follows  very  closely  section  3  of  the  paper 
cited  in  Appendix  2  C„ 

We  have  seen  that  when  testing  :  0  »  ©q  against  any  simple 
alternative  ©  *  0,(0.  <  C  )„  the  "beat"  region  of  acceptance  for  H 
(in  the  sense  of  Neyman  and  Pearson),  based  on  the  first  r  out  of  n 
ordered  observations  from  an  exponential  distribution, is  of  the  form 

r 


$  „  >  where  $ 
r,n  r,n 


C  ) _  x,  „  +  (n-r)x„  „)/r  in  the 

&  i*n  r’n 


non-replacement  case  and  §  „  =  nx _ /r  in  the  replacement  case. 

" on  r^n 

One  could  interpret  the  decision  rule  >  C  to  mean  that  we 

wait  until  time  _  (the  time  when  the  r**1  failtire  occurs),  then 

rtn 

compute  and  make  the  appropriate  decision,  however,  in  the  event 

r0n 

that  we  are  able  to  observe  the  life  toet  continuously t  this  clearly 

wastes  information.  Indeed,  we  assert  that,  if  continuous  observation 

is  t alien  into  account,  we  can  frequently  shorten  the  waiting  time  to 

reach  s  decision  and  reduce  the  number  of  items  failed.  To  see  this  we 

note  that  „  >  C  becomes  D  x.  „  ♦  (n-r)  x_  1  >rC  in  the 
r,n  jSl  r,a 


B.k5 


non-replacement  cut  and  ax_  _  >  rC  in  the  replacement  oa se. 

*  y  ll 

Bat  tV "  x.  _♦  (n-r)  x _ ]  la  the  total  observed  life  up  to  time 


fcr 


r.n 


*r  a  in  the  non-replacement  case  (note  that  >  x^  Q 


r 

1*1  i,fJ 


is  the  total 


life  of  the  r  items  which  failed  and  (n-r)x..  _  is  the  amount  of 

r,n 

time  lived  by  the  (n-r)  items  which  did  not  fail)  and  nx„  is 

rtn 

the  total  observed  life  up  to  time  x  in  the  replacement  case 

(note  that  in  the  replacement  case  n  items  are  constantly  on  test  for 

a  length  of  time  x  )o  Thus  accepting  H  when  $  „  >  C  ie  equiv- 

ren  o  r,|B 

lent  to  accepting  H  if  the  total  life  obsorved  up  to  time  x  is 

o  P 

greater  than  rC.  Suppose  now  that  at  some  moment  t  there  are 
exactly  k  failures,  0  <  k  <  r-1,  and  that  the  observed  total  life 


7(t)  given  by  V(t)  «  > _  x,  +  (n-k)t  in  the  non-replacement  case 

i*l  i»n 

and  by  V(t)  =  nt  in  the  replacement  case  exceeds  rC  (Note  that  in 


k 

the  non-replacement  ce.se  > _  x.  is  the  amount  contributed  to  V(t) 

i«l  * 

by  the  k  Items  which  failed  by  time  t  and  (n-k)t  is  the  amount 
contributed  by  the  (n-k)  items  which  have  not  failed.*  In  particular 0 

5  r 

if  ft  *  X_  .  '  then  V(x.  )  =  )  x .  ♦  (n-r)x„  o  The  formula  for 

r^n1  r,n  lvn  r,n 

7(t)  in  the  replacement  case  is  obvious *,)  Since  V(t)  is  monotonically 

increasing  in  t,  we  know  that  V(xr  )  >  V(t)  >  rC,  and  therefore  we 

should  stop  the  life  test  at  time  t  and  accept  Eqo  More  generally  a 

deolsion  rule  having  preciaely  the  same  O.Cc  curve  as  >  C  but 

rtn 

requiring  on  the  average  fewer  failures  and  a  shorter  decision  time  is 
hasad  on  terminating  at  total  observed  life  *  oin  (V(x_  )„  rC)  (where 
both  r  and  C  are  preassigned).  If  the  experiment  is  terminated  at 


2.k6 


total  lifo  V<xp^a)(i.e.t  If  the  total  Ufa  required  to  observe  r 

failures  is  <  rC)f  then  the  aotion  in  terms  of  hypothesis  testing  is 

the  rejection  of  the  null  hypothesis.  If  life  testing  is  terminated 

with  total  life  ■  rC  (i.e.,  if  V(x  ).  the  total  lifo  required  to 

r,n 

observe  r  failures*  exceeds  rC) ,  then  the  action  taken  is  to  accept 

the  null  hypothesis,  (Note  that  in  the  replacement  case  (V(x  )t  rC) 

r  ,n 

becomes  min  (x  j  where  x  is  the  time  of  the  failure 

r  *n  n  rfn 

and  is  a  truncation  time.) 

Described  in  more  detail  the  decision  rule  is  as  follows: 

(a)  Continue  life  testing  so  long  as  V(t)  <  rC  and  0  <  k  <  r-l. 

(b)  Stop  experimentation  at  time  t  with  acceptance  of  Hq  as 
soon  as  V(t)  >  rC  and  0  <  k  <  r-l. 

(c)  Stop  experimentation  at  time  xp  q  with  rejection  of  Hq 

if  V(t)  <  rC  for  all  t  <  x_  <,(Note  that  acceptance  of  H  takes 

*  o 

place  bo  tween  failure  times  t  and  always  before  time  x  ,,) 

r6r* 

We  now  proceed  to  find  some  useful  properties  of  the  truncated  rule 
based  on  V(t).  To  find  these  properties*  we  remark  that  (defining  x 

oBn 

as  zero) 


(2D.1) 


r 


x  +  (n-r)  x  (n-i+lKx.  ~x,  .  ) 

ian  r,n  i  ,n  x-lsn 


in  the  non-replacemen b  case  and 


(2D, 2)  nx 


L  „  *  / _  n(x.  -x,  ,  ) 

r»n  1«1  i,L  i"1,n 


in  the  replacement  case. 

Introducing  (as  ?as  done  in  Appendix  2A)  new  random  variables  defined  by 


(2D. 3) 


yl  “ 


“I.. 


and  y{  «  (n-i+lXx^-x^^h 


2,3,, 


in  the  non-replacement  case  and 


(2D04) 
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xu^  ^  and  y^  >  n(x, 


..  _**x,  t  _)u  i  *  2,3i««tjr 
i«n  i-l,n 


la  the  replacement  cast. 


(2D<>5) 


V(tu)  >  rC  caa  be  written -as 
r£- 


i«l 


7-  >  rC, 


We  saw  in  Appendix  2A  that  the  y^  are  mutually  independent  xrndom 
▼nriables,  each  distributed  with  common  p.,  d,f«  ~  x  >  O*  0  >  0o 

If  we  interpret  y^  &s  the  time  interval  between  the  (i-l)Bt  and 


,th 


event  in  a  Poisson  proeesa  having  mean  occurrence  rata  A  *  1/0 v 


it  is  clear  that  j  y.  >  rC9  if  and  only  if  k»  the  number  of  events 

i*l  1 

in  a  time  interval  of  length  rC4  is  0  <  k  <  r~l„  If  the  number  of 

r 

events  in  such  an  interval  ia  >  r  „  then  V y  <  rC„  Thus  the  probe- 

i«l  “ 

bility  of  reaching  a  decision  requiring  exactly  •  f  «  k  failures  is 


(2Dn6) 


Pr(  k  8)  *  p(k  k  «  0,1  *2,. ,r-l 


Pr(  f  »  r|o)  «  1  p(k,V^)  »  1  -/^(r-l;/^), 

where  /U^  •  rC/©0  The  expected  number  of  observations  to  reach  a  decision 
is  given  by 

(2D„7>  E(j(f)  »  JL^  XPrif  .  k|0)  »  ^7(r-2i^)  ♦  rCl-y^ (r-li^]. 

Eg(Vfc)),  the  expected  total  life  in  reaching  a  decision, 

is  given  by 

(2Do8)  EQ(?<t))  -  0Eq(  J1)  e 


Bq(T) t  the  expected  waiting  time  to  reach  a  decision,  is  given  by 


(2D. 9) 


V«  '  Sr  Pr<p*  kl»)  Wa' 


where  (Xk>a)  :  © 


in  the  non-replacement  case  and 


W  :k0  /n  In  the  replacement  case.  In  the  replacement  case 


(2D. 10) 


e0(t)  »  y“-  Pr(p  =  k|e)^ 

**  E_(V(t)) 

:  £  s(p)  =  - 


Finally  L(©  ),  the  probability  of  accepting^  ■Qo  when  Q  is  true, 
is  given  by  L  (0 )  *  *(r-lj  p^).  Note  that  in  the  replacement  case  (25)_ 
through  (28)  coincide  with  what  we  have  Just  done  if  we  set  Tq  ■  rC/n 
and  r  ■  r  •  If  this  is  done  then  ■  nT0/ Q  -  tC/q  ~  Pq  * 

o  v  K 


Remark:  In  the  above  we  considered  a  test  based  on  &  r  n  as  a  ^'ruIlca^e<^ 
test.  This  involved  consideration  of  total  life.  The  essence  of  what 
was  said  is  a  special  case  of  the  following:  Suppose  that  the  experimenter 
wishes  to  expend  no  more  than  total  life  V*  in  experimentation  and  that 
he  employs  the  following  rule  of  action:  Reject  if  rQ  failures  occur 
before  total  life  V*  has  been  used  up;  accept  if  fewer  than  rQ  failures 
occur  by  the  time  one  has  observed  a  total  life  of  V*.  In  the  event  that 
one  rejects,  experimentation  stops  at  V(  ) ,  the  total  life  observed 
up  to  and  including  T  ,  the  r  th  failure  time.  In  the  event  that  one 
accepts,  the  toted  life  observed  will  be  V*.  It  follows  directly  from 
the  properties  of  Poisson  processes  that  the  probability  of  reaching  a 
decision  requiring  exactly  P  ■  k  failures  is 
(2D.ll)  Pr (p  *  kjo)  s  p(k;  Hq),  k  *  0,l,2,...,ro-l 


r0-l 

Pr(p  =  r0  |  0)  ■  1  -  y~~  P (kJ  M9)  *  1  " 


2.4 9 


wher*  ^  •  V*/0.  'fhc  expected  number  of  observations  to  reach  a  decision 
is  given  bj 

(2D.12)  E0(/)  »  ♦  r  Cl  -T(r-l;^)Jo 

S^(7(t))f  tha  expected  total  life  in  reaching  a  decision  is  givon  by 

(2D.13)  E0(V(t))  *  0Ee(/) 

and  L(0),  the  probability  of  accepting  0  *  0Q  when  0  io  true, is 
given  by  p  ml 

(2D.14)  L(0)  »  x?~~  p(k}/^)  «  T(ro-lj>^)«. 

k«o 

The  considerations  involving  “6  are  a  special  case  of  what  we  have 

r,n 

just  done),  with  V*  =  rC  and  rQ  =  rQ 


As  an  illustration  of  the  theory  presented  in  Section  3  and  Appendix 
2  D  are  consider  three  test  procedures  which  have  virtually  the  same 
operating  characteristic  curve*  Specifically  it  is  assumed  that  we  wish 
to  test  H0:®0  ®  1500  hours  against  =  500  hours  with  *  .05 

i.e.(  we  want  L(Qo)  ~  l-<  *  .95  and  L(0^)  a  JS  *  .05  (actually  wa 
have  to  be  satisfied  with  L(O^)  <  .05).  The  three  procedures  are: 

(a)  20  items  are  taken  at  random  from  the  lot  and  placed  on  life 

test.  Items  which  fail  are  not  replaced.  At  each  moment  t,  compute 
the  total  life 

V(t)  a  y _  x.  ..  ♦  (n-k)  t.  where  k  is  the  number  of 

W  l»n 

failures  which  have  occurred  before  time  t  and  n  *  20  (if  i  »  0t 
define  total  life  aB  nt).  If  V(t)  exceeds  8150  for  any  k,  0  <  k  <  9» 
stop  the  experiment  at  tine  t  and  accept  Hq(0  *  1500  hours)*.  Otherwise 


2.50 


the  act ion  taken  i*  to  re ject.  This  teat  i a  equivalent  to  accepting 

<  815 »  (Prom 


Ho  if  *10,20  >  815  <wd  n0  if  & 


10,20 

Table  5  we  see  that  if  ■  3,*t  ■  .05,  then  r  *  10  and 

^2r)/2r  »  <>54263  Therefore  the  aeoeptance  region  is 
*10,20  >  <1500(,5426)  »  815.  ) 

(b)  20  items  are  taken  at  random  from  the  lot  and  placed  on  test, 
railed  iteas  are  not  replaced.  If  min  CX^  2Q,  5403  *  540  (i.e,„  the 
tenth  failure  occurs  after  540  hours )* truncate  the  experiment  at  540 


hours  with  acceptance  of  H 


If  *la  “lO.ZO’5403 


^0  2q  ( i « e .  ^  the 


tenth  failure  occurs  before  540  hours),  truncate  the  experiment  at  X1Q  2Q 
with  the  rejection  of  (Proa  Table  4a,  using  r  «*  10,  n  »  2r  *  20„ 

we  see  that  the  truncation  time  To  *  1500(,363)  =  540c) 

(c)  20  items  art*  taken  at  random  from  the  lot  and  placed  on  test* 


An  item  which  fails  is  replaced  at  once  by  a  now  item  from  the  original 

lot.  The  time  "hen  the  i**1  failure  occurs  is  measured  from  the 

beginning  of  experimentation,  If  min  CX1Q  20«407.53  *  407.5,  truncate 

the  experiment  at  407«5  hours  with  the  acceptance  of  E  >  If  oin 

o 

^S.0,20****^*^  “  X10,20*  truncat*  experiment  at  X1Q  2Q  with  the 

rejection  of  H  «  (In  the  replacement  case  the  truncation  time  T  is 
0  o 

given  by  «0Xf.(K.(2r)/2n,  This  gives  407*5  for  the  value  in  this 
problem.) 

In  the  table  bolcw  we  give  L(Q),  EQ(r)t  and  EQ(T)  for  the 
tests  A,  B,  and  C  for  selected  values  of  Q„ 
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Detailed  proof*  of  the  result*  on  sequential  Ilf*  tost*  sketched  in  flection 
k  of  Chapter  ZZ  are  given  in  the  following  reference}  B.  Spate  in  and  M.  86bel, 
"Sequential  Life  Tests  la  the  Ixponentlsl  Case, "  Annals  of  M&theaatlcal  Statistics 
86,  88*93*  1955*  In  flection  8  of  this  reference  one  will  find  a  derivation  of 
formulae  (31)  through  (39)  Inclusive.  In  flection  3  of  the  reference  the  basic 
Identity  *g(V(t))  «  e  *Q(r)  ,  relating  the  expected  sonant  of  total  life  observed 
la  reaching  a  decision  and  the  expected  number  of  failures,  is  derived.  This  formula 
holds  In  general,  whether  or  not  items  on  test  are  replaced. 

Appendix  80 
1.  Introduction 

It  la  Interesting  to  ask  the  question:  Bow  will  truncation  of  the  sequential 
life  test  affect  the  type  I  error  O  and  type  II  error  0  Y  We  know,  froai  con¬ 
siderations  analogous  to  those  of  A.  Wald,  that  the  sequential  life  test  procedure 
based  on  using  (31)  will  eventually  terminate.  But  this  any  be  inordinately  ex¬ 
pensive  la  terms  of  either  the  tine  involved  in  the  life  test,  or  In  terms  of  the 
number  of  items  failed,  or  both,  there  are  nany  situations  where  it  is  desirable 
and  even  necessary  that  we  place  a  definite  upper  Halt  on  either  the  number 
of  items  failed  or  on  the  total  length  of  the  life  teat  (or.  If  necesaary  on  both). 

In  what  follows  we  study  how  much  one  changes  the  type  I  and  type  II  errors,  if 

vferiow.fr 

one  truncates  the  sequential  life  test  in  gngssmne  ways. 

Remark:  Prom  this  point  on  wi  follow  closely  considerations  in  Wald's  book, 

pp.  61*65* 

\ 
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2.  Truncation  on  the  nunber  of  items  failed. 

Suppose  firnt  that  w*  Mt  a  definite  Unit.  rQ  ,  on  the  number  of  items 

failed.  We  can  achieve  this  by  truncating  the  sequential  life  teat  at 

r  m  rQ  ,  i.e.,  by  giving  a  new  rule  for  the  acceptance  or  rejection  of  H0:  0  ■  9q 

when  rQ  failures  have  occurred  if  the  sequential  life  test  did  not  lead 

to  a  decision  for  r  <  rQ  .  A  simple  and  reasonable  truncation  rule  after  the 

r^th  failure  is  the  following:  If  the  sequential  probability  ratio  test  given 

by  (31)  does  not  lead  to  a  decision  for  r  <  r  ,  accept  H  :  0  •  9  after  the 

—  0  00 

r  °th  failure  has  occurred  if 

o 

1  1 

(26.1)  log  B  <  r  log  -2  .  (JL  _  i)  V(x  )  <  0 

c  ox  0o  rQ 

and  reject  Hq  (accept  :  0  •»  8^)  after  the  rQ'th  failure  has  occurred  if 

(20.2)  0  <  r  log  -2  -  (i.  -  i.)  V(x  )  <  log  A  . 

o  0X  91  eQ  rQ 

Truncating  the  sequential  life  test  after  rQ  items  fail  will  change  the 

Type  I  and  T^rpe  II  errors.  They  will  no  longer  be  a  and  3  ,  the  Type  I  and 

Type  II  errors, respectively, in  the  untruncated  sequential  case.  The  effect  of 

the  truncation  on  the  a  and  0  depends,  of  course,  on  tq  . 

The  larger  one  makes  rQ  ,  the  smaller  are  the  truncation  effects  on  a 

and  0  .  let  us  denote  the  resulting  Type  X  and  Type  XX  errors  as  a(ro)  and 

0(r  )  ,  respectively,  if  the  sequential  life  test  is  truncatsd  at  r  ■  r 
o  o 

failures  at  the  latest.  We  now  derive  upper  bounds  for  a(rQ)  and  0(rQ) 

To  obtain  an  upper  bound  for  a(rQ)  vs  have  to  consider  the  cases  in 
which  th*  truncated  life  test  leads  to  the  rejection  of  Bq:  9  ■  ,  while 

the  nan-truncated  sequential  life  test  lends  to  the  acceptance  of  H0  . 


2.55 


Supposing  that  Hq  hold*  (i.e.,  that  6  •  dQ)  ,  let  p@(rQ)  bt  the  probability 

that  the  sample  random  function  associated  with  a  Ufa  teat  is  such  that  the 

truncated  life  test  leads  to  rejection  of  H  ,  while  the  n on-truncated  life  test 

o 

leads  to  the  acceptancs  of  HQ  .  Clearly  we  see  that 
(20.3)  a(r0)  <  a  +  P0(ro)  . 

The  reason  for  the  inequality  rather  than  the  equality  Is  that  there  my  be 
sample  raadoa  functions  associated  with  life  tests  for  which  the  truncated  life 
test  leads  to  acceptance  of  HQ  ,  while  the  non-truncated  life  test  leads  to  the 
rejection  of  Hq  .  To  obtain  an  upper  bound  for  a(rQ)  >  «  need  mrely  derive 
an  upper  bound  for  o  (r  )  Assuming  that  B  is  true,  p  (r  )  is  the 

probability  that  the  random  function  associated  with  a  life  test  is  such  that 
the  following  three  conditions  hold  simultaneously: 

(2Q.b) 

(i)  log  B  <  r  log  s—  -  (jr-  -  £-)  V(t)  <  log  A  ,  for  r  -  1,2,..., r  -1 

®1  ®1  9o  ° 

and  for  all  t  <  x  ; 

0  9 

(ii)  0  <  r  log  ~  -  (i  -  ~)  V(x  )  <  log  A  j 
°  90  ro 

and 

(hi)  When  the  sequential  life  test  is  continued  beyong  the  rQ’th  failure, 
it  terminates  with  the  acceptance  of  Bg  . 

Assuming  that  HQ  is  true,  let  Pa(r0)  k*  the  probability  that  condition 
(11)  holds,  l.e.. 


<W-5)  »„<*„)  -  »  [o  «  r0  lo*  -  i)  »(^)  <  lo«  *j»  -  •„  J 
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8iae«  the  probability  that  (ii)  ia  fulfilled  cannot  ba  — Mtr  than  tha  pro- 

babiiity  that  conditions  (i)  ,  (ii)  and  (iii)  art  fulfilled  sinultansoualy, 
we  hare 


(a0*6)  Po<ro>  *  Po<ro> 

and  therefore, 

<20-7)  a(ro>  <a  '  p0(r0) 

°  *  Potro>  *•  “  u»*r  bound  for  o(ro)  .  lib  shou  further  on  that 
po(rQ)  can  ba  confuted  easily. 

To  obtain  an  upper  bound  for  P(rQ)  ,  let  ua  assune  that  9  -  ia 
true  and  let  then  be  the  probability  that  the  truncated  life  teat  leads 

to  the  acceptance  of  Hq  while  the  non-truncated  life  teat  leads  to  the  re¬ 
jection  of  Ho  .In  other  words,  ^(ro)  la  the  probability  (assuming  that 
“  dl  1*  trv*)  that  the  sample  random  function  associated  with  a  life  teat 
ia  such  that  the  following  three  conditions  bold  alaultaneously: 


(i)  log  B  <  r  log  !2  -  -  i.)  v(t)  <  log 


for  r.  1,2 


,  e  e  e. 


r  -  1 


and  for  all  t  <  xy  j 

(»•*)  °  , 

<U)  log  1  <  r„  log  -2  -  <i.  -  *.)  vU  )  <  o  , 

*  1  O  o 

(iii)  If  the  sequential  life  teat  is  continued  beyond  the  rQ'th  failure, 
it  terminates  with  the  acceptance  of  Bjj  e  •  9^  . 

Clearly 

(80 


<Xre)  <  »  ♦  ft(r0)  . 
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8 lac*  it  Is  difficult  to  detend.no  ( rQ )  ,  v*  give  a  single  upper  bound  first. 
A— uuing  that  B^  la  true,  let  ^(r0)  be  the  probability  that  condition  (11) 
holds,  i.e.. 


(20.10) 


Pr  {  log  B  <  rQ  log  g~ 


^-)  v(xr  )  <  o|e 


Then  p,  (rQ)  >  ^(r^  and  hence 

(20.11)  P(rQ)  <  8  +  p^o*  . 


We  now  show  how  to  coapute  pQ^r9)  and  ^(r^  .  To  coapute  pQ(rQ)  ,  we 
recall  that  If  tq  is  preassigned,  then,  under  the  hypothesis  that  Hq:  9  •  0Q 

A 

Is  true.  2V(Xy  )/bq  Is  distributed  as  7Cr(2rQ)  .  Consequently. 


(2C.12) 


•  Pr 


log  k  -  log  A) 
T-rj - 


< 


2rol°8k-) 
■k'-r-  J 


where  k  -  9  Jb^ 

In  a  slnllar  way  one  can  conpute  f  ^0)  If  rQ  Is  preassigned,  then, 
under  the  hypothesis  that  B^:  9  ■  01  Is  true.  2V(xr  )/b1  Is  distributed  as 
7^(2rQ)  .  Consequently 


(20.13) 


k  log 


k 


2k(rQ  log  k  -  log  B 

- - m — 


J 


Thus  we  can  suanarlse  our  results  as  follows: 


a(ro)  <«(.♦  Pr 


log  k  -  log  A) 
“TH - 


<Tt?<*0> 


2ro  iog  k  •) 

c-rrr-j 


(20.1k) 
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It  is  our  fueling  that  the  upper  bounds  for  a(rQ)  end  d(rQ)  that  we  have 
obtained  are  substantially  above  the  true  values  of  a(rQ)  end  0(rQ)  . 

It  seeaa  appropriate  at  this  point  to  give  a  nunerical  example.  Consider 
the  problSs  o t  testing  HQ:  9  »  9Q  against  H^s  9  »  0^  with  a  ■  fi  •  .05  «®d 
k  •  0^/0^  >  3  .  one  can  readily  verify  from  Table  5,  that  a  non-eequent’a1.  ^Ife 

test  requires  r  ■  10  .  She  hypothesis  that  9  «  6  is  accepted  if 

■  s'  *  • 

V(Xio  n)  >  ®0  X  ,^(20)/2  n  5.4260 0  >  and  rejected  otherwise.  Nov  let  us  see 

what  happens  in  a  sequential  life  teat  for  the  four  truncation  nunbcar  of  f  Allures, 
r0  -  10  ,  15  ,  20  ,  30  . 

If  rQ  m  10  ,  then 

a(io)  <  a  ♦  Pr  {.10  log  3  -  log  19  <  75(20)  <  10  lag  3} 

-  a  ♦  pr  {8.05  <  7?(20)  <  10.99}  &  .05  ♦  .04  »  .09 


0(10)  <  0  ♦  Pr  {JP  lug  3  <7C?(20)  <  3  (10  log  ?  ,  log  19)} 

-  0  ♦  Pr  {32.96  <T?(20)  <  41.79} oe  .05  ♦  .04  .  .09  . 

If  r  -  15  ,  then 
o 

0(13)  <0  ♦  PS  {15  log  3  -  log  19  <*)?( 30)  <  15  log  3} 

’  °  *  *  C13-*  <  <0»)  <  -  .05  ♦  .0*9  .  .07J 


and 


2,i;> 

9(19)  <  3  ♦  Pr  {43  log  3  <7^(30)  <  3(15  log  3  ♦  log  i9)} 

•  3  +  Pr  {49.44  <7(?(30)  <  58*26}  *05  ♦  .015  -  .065  * 

If  rQ  ■  20  ,  than 

0(20)  <  a  ♦  Pr(20  log  3  -  log  19  <'X?(40)  <  20  log  3) 

-  a  ♦  Pr(19.04  <  ?C2(40)  <  21. 98)  £*  .06 

and 

3(20)  <  3  +  Pr  {60  log  3  <^f(4o)  <  3(20  log  3  ♦  log  19)} 

-  3  ♦  Pr  {65.92  <7?(40)  <  74.76}  2*  .055  * 

If  rQ  -  30  ,  than 

a(30)  <  a  ♦  Pr(30  log  3  -  log  19  <7C2(60)  <  30  log  3) 

-  O  ♦  Fr( 30.04  <X2(6o)  <  32.96.  &  *051 

3(30)  <  3  ♦  Pr  { (90  log  3  <*><?;60)  <  3(30  log  3  ♦  log  19)} 

•  3  ♦  Pr(98.88  .<o£(60)  <  IO7.7O)  A*  .051  . 

Thus  vs  sea  in  this  'staple  that  if  va  truncate  ttw  sequential  life  teat 
at  rQ  -  30  ,  l.e.,  at,  3  tint*  10  ,  tbe  r  required  for  the  non- sequential  Ufa 
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teat,  than  o(r0)  la  approximately  equal  to  a  and  0(rQ)  la  approximately 
equal  to  P  .  Tables  are  being  calculated  for  other  value  a  of  a,0  and  & 
and  the  indications  are  that  what  we  observed  in  the  example  holds  nor*  generally. 
That  la,  truncation  of  the  sequential  life  test  at  three  tinea  the  nunber  of 
failures  required  in  the  non- sequential  life  teat  will  have  virtually  no  effect 
on  either  a  or  P  . 

3*  Truncation  on  the  total  observed  life. 

We  have  up  to  this  point  truncated  the  sequential  life  test  by  setting  a 

definite  Unit  on  rQ  ,  the  nunber  of  ltens  failed.  We  now  wish  to  truncate 

the  sequential  life  test  by  placing  a  definite  Unit  VQ  on  the  total  observed 

life.  We  can  achieve  this  by  truncating  the  sequential  life  teat  at  V(t)  »  VQ  , 

l.e..  by  giving  a  new  rule  for  the  acceptance  or  rejection  of  H  :  9  ■  0  when 

o  o 

V(t)  »  Vq  if  the  sequential  life  teat  did  not  lead  to  a  decision  for  V(t)  <  Vq  . 
A  staple  and  reasonable  truncation  rule  at  total  life  VQ  is  the  following: 

If  the  sequential  probability  ratio  test  given  by  (31)  does  not  lead  to  a  final 
decision  for  v(t)  <  ,  accept  HQ:  9  »  8q  et  total  life  VQ  if 

(20.16)  logB<r  log  i  )  Vo  <  0  . 

1  1  0 

and  reject  HQ  (accept  H^:  9  ■  9X)  if 

1  t 

(80.17)  0  <  r  log  -2  -  (i-  -  ~)  V  <  log  A  . 

*1  *\  %  0 

Truncating  the  sequential  life  teat  at  total  life  VQ  will  change  the  Type  Z 
and  Type  II  error*.  They  will  no  longer  be  o  and  0  ,  the  Type  I  and  Type  IZ 
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errors,  respectively,  la  the  untruncated  sequential  case.  The  effect  of  the 

truncation  oo  a  and  8  depends,  of  course,  oo  Vq  . 

The  larger  om  makes  V  ,  the  smaller  art  the  truncation  effects  on  O 

sad  B  .  Let  us  denote  tbs  resulting  Type  I  end  Type  II  errors  ss  o(VQ)  sad 

P(VQ)  ,  respectively,  If  the  sequential  life  test  is  truncated  at  total  life 

V(t)  ■  Vq  at  the  latest.  We  now  derive  upper  bounds  for  a(VQ)  and  3(Vq)  . 

To  obtain  an  upper  bound  for  a(^0)  we  have  to  consider  the  carea  in 

which  the  truncated  sequential  life  test  leads  to  the  rejection  of  B  :  e  ■  9 

o  o 

while  the  non-truncated  sequential  life  test  leads  to  the  acceptance  of  Hq 
Supposing  that  Hq  holds  (l.e.,  that  6  ■  0 )  ,  let  pq(Vq)  be  the  pro¬ 
bability  that  the  sample  random  function  associated  with  a  life  test  la  such  that 
the  truncated  life  test  leads  to  rejection  of  Hq  ,  while  the  non-truncated 
life  test  leads  to  the  acceptance  of  Hq  .  Clearly  we  get 

(20.18)  a<v0Ua  +  fo(?o>  * 

The  reason  for  the  inequality  rather  than  the  equality  Is  that  there  aay  be 
saunple  randan  functions  associated  with  life  tests  for  which  the  truncated 
life  test  letds  to  acceptance  of  BQ  ,  while  the  nan-truncated  life  test  leads 
to  the  rejection  of  HQ  Assuming  that  BQ  is  true,  .  ;(Vq)  is  the 
probability  that  the  random  function  associated  with  a  life  test  is  such  that 
the  following  three  conditions  hold  simultaneously: 

(1)  log  B  <  r  log  -2  -  (f-  -  i)  ?(t)  <  log  A  ,  for  7(t)  <  V 

®1  *1  8o  ° 

,  i 

.  „  (11)  0  <  r  log  s2  -  -  i)  ?  <  log  A  and 

(20.19)  *1  ®1  V  ° 

When  the  sequential  life  test  is  continued  beyond  total  life 
Vq  ,  it  terminates  with  the  acceptance  of  BQ 


(iii) 
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Assuming  that  HQ  la  true,  let  po(VQ)  be  the  probability  that  condition 
(ii)  holds,  i.e., 

(20*20)  p0(V0)  .  fr  ^0  <  r  l0«  ^  -  (i  -  i)  V0  <  los  *  j9  .  6,^  . 

Since  the  probability  that  (ii)  la  fulfilled  cannot  be  smaller  than  the  pro¬ 
bability  that  conditions  (l)  ,  (ii)  ,  and  (ill)  are  fulfilled  simultaneously, 
we  hare 

(“•2D  fo<To>  i  f>„(Yo> 

and, therefore, 

(20.22)  o(VQ)  <  o  ♦  po(yQ)  . 

Thus  o  +  ?  (Vo)  is  an  upper  bound  for  a(VQ)  .  We  show  further  on  th.'s.t 
po(Vo)  can  be  computed  easily. 

To  obtain  an  upper  bound  for  0(VQ)  ,  let  us  assume  that  E^:  0  ■  0j  ,  is 
true  and  let  p1(VQ)  then  be  the  probability  that  the  truncated  life  test  leads 
to  the  acceptance  of  Hq  ,  while  the  non-truncated  life  test  leads  to  the 
rejection  of  BQ  In  other  words,  p^(VQ)  is  the  probability  (assuming  that 
0*9^  is  true)  that  the  saaple  random  function  associated  with  a  life  test  is 
such  that  the  following  three  conditions  hold  simultaneously: 

0_  t  i 

log  B  <  r  log  -  (±-  -  i-)  V(t)  <  log  A  ,  for  V(t)  <  V 
91  el  ®o  ° 

lo*  *  <  r  1-*  ^  -  (^-^)vo<0 


(1) 

(20.23) 

(ii) 
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(20.23) 

.  and 

(ill)  If  the  sequential  life  test  is  continued  beyond  total  life  Vq  , 
it  terminates  vltb  the  acceptance  of  H^:  9  ■  9^ 


(20.24)  Clearly  P(VQ)  <  6  *  h<Vo> 


Since  it  is  difficult  to  determine  p^(VQ)  we  give  a  simple  upper  bound  first. 
Assuming  that  is  true,  let  ft(Vo  )  be  the  probability  that  condition  (li) 
holds,  l.e.. 


(20.25)  P]L(Vo)  -  Pr  {log  B  <  r  log  ^  -  (j-  -  VQ  <  0  |  9  -  ej  . 
Then 

(20.26)  fcCV  >  fX(V0) 

and  hence 

(20.27)  0(Vo)  <  P  +  p1(VQ)  . 


Ve  now  show  how  to  compute  po(VQ)  and  p^(VQ)  .  To  compute  j>Q(V0)  > 
we  recall  that  if  VQ  is  preassigned,  then  under  the  hypothesis  that 
Hq:  9  m  dQ  is  true,  we  are  observing  a  Poisson  process  with  rate  parameter 
\Q  ■  l/0Q  for  a  length  of  time  Vq  .  Consequently, 
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where  k  ■ 


0.  V«  / 
j-  ,  »o  -  (k-1)  y-  J  log 


k  and 


“o- 


(k-1)  ~  ♦  log  A 
Bo 


/ 


log  k 


la  a  similar  way  one  can  compute  p1(VQ)  .  If  VQ  la  p  reassigned,  then  under 
the  hypothesis  that  H^:  0  ■  0^  is  true,  we  are  observing  a  Poisson  process  with 
rate  parameter  ■  l/fl^  for  a  length  of  time  Vq  .  Consequently, 


(20.29)  ^(V^)  .  Pr 


L  log  k 


s-  +  log  B 


(k-1) 


<  r  < 


log  k 


■  zZ  ?(*> 


■1 


<  r  <  n. 


0, 


where 


“l 


-  j(k-l)  y*  +  log  .]/ log  k  and  nx  -  (k-1)  ~  j  log  k  . 


Thus  we  cam  summurlse  our  results  as  follows: 


(20-30) 


a(V0) 


- 1  v 

<  a  ♦  P(n  ~) 


m  <  r  <  n 
o  o 


and 


V  1  v 

(20.31)  P(v  )  <  0  +  2. _ »  s2) 

v  m^  <  r  <  n^  1 

where  m  ,  n  ,  m^.  ,  n^  are  defined  above.  We  are  quite  sure  that  the  upper 
bounds  for  a(VQ)  and  P(VQ)  that  we  have  obtained  are  substantially  above  the 
true  values  of  <*(V0)  and  P(VQ)  . 

WS  now  give  a  numerical  example  to  illustrate  what  we  have  Just  discussed. 

Consider  the  problem  of  testing  HQ:  0  -  0Q  against  H^r  0  ■  0^  ,  with  a  •  p  .  .05 

and  k  m  0  /0,  ■  3  .  Let  us  see  vhat  happens  for  truncation  times  V  such 
01  o 

that  VQ/0o  .3, 6, 9, 12, and  15. 

If  Vo/9q  ■  3  ,  then 


A*  -  •«**  —  * 
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0(0  <  <*  ♦  ^  1  »(*!  3) 

■  TTTa 
o  o 

y 

■o  .  (k-l)  ^  j  log  k  .  6/1.10 


,  where 
9  \ 


-  5.45 


Bence 


Similarly 


n° "  [k  "1)  K  * 108  A] 

y/  log  k  .  -t:i^  -  y$b  -  B-13 

a(VQ)  <  *°5  +  ^ 

}  p(r;  3)  »  *05  ♦  .08  ■  .13 
> 

0<VQ)  <  p  +  ^ 

■— > 

p(r;  9)  ,  where 

r  v 

^“l 

J  j  log  k  .  [6  -  2.*]/l.l0 

»  |(k-l)  g-  ♦  log  B 

.  3.06/I.IO  -  2.78 


and 


Hence 


V  , 

1^  .  (k-l)  ~  /  log  k  -  5.V5 


P(V  )  <  ♦  E  P<**5  9)  •  .16 


If  5—  ■  6  ,  then 

w  _ 


C*(V  )  <  a  ♦  p(r;  6)  m  .05  ♦  .039  -  .009 


6(9.)  <  e  ♦ 


6 ,(ri 


10)  •  .05  ♦  .023  -  .073 
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If  x—  ■  S>  »  then 
9o 


a(ve)  <  a  ♦ 


p(r;  9)  -  .060 


and 


P(VG)  <  P  ♦  >  I  p(r;  27)  ■  .064 


v 

If  s—  ■  12  ,  then 
o 


a(VQ)  <  a  +  p(r;  12)  .  .055 

r. 


and 


P<V0)  <  P  ♦  2_  P(ri  36)  -  .054 


If  «•  15  «  then 


0(VQ)  -  °  P(r*  15)  -  .052 


and 


P(v0)  <  P  ♦  Zi  P<*i  45)  -  .052 
r-25 

fhue  we  see  In  thle  example  that  if  ve  truncate  the  eequential  life  teat 
at  V  »  15  0  .  then  a(v  )  and  P(?M)  are  approximately  equal  to  a  and 

OO  w  O 
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£  t  respectively.  For  tho  non- sequential  life  toot,  truncation  occurs  At 
7  ■  $ ,1*269 Q  •  Thus  in  this  example  truncating  tho  sequential  life  toot  at 
three  tines  the  7  required  for  the  non-aequential  test  has  virtually  no 
effect  on  either  cc  or  •  Tables  are  being  calculated  for  other  values  of 
of ,  and  8 ^8^  and  there  are  indications  that  what  we  observed  in  the 
exanple  holds  nore  generally. 


Truncation  on  the  number  of  itsws  failed  and  total  life . 

Nov  it  nay  happen  that  we  would  like  to  truncate  the  life  test  both  with 

respect  to  the  number  of  failures  r  and  total  life  7  •  We  first  note  that 

o  o 

our  truncated  sequential  life  tests  considered  up  to  now  are  of  this  kind. 

Indeed,  suppose  that  one  truncates  at  r  ■  rQ  .  Then  we  assert  that  this  induces 
a  truncation  on  total  life. 


ro  1o*  r 

T.'-r — r-’Sro  • 


[See  Figue  (5).] 


Thus  if  the  randan  Amotion  representing  the  life  test  is  such  that  one  attains 
total  life  7  (reaohes  BC)  with  fewer  than  r  failures,  than  one  knows  that  if 
one  continues  the  test  until  rQ  failures  occur  then  the  sample  random  function 
snot  cross  either  BD  or  DC  and  in  either  case  we  would  accept  HQ  .  Banco  if 
one  attains  total  life  V0  ■  a  rQ  with  fewer  than  rQ  failures,  one  can  step 
with  acceptance  of  H  ,  Similarly,  if  one  truncates  at  total  life  7  »  7  ,  them 
we  assert  that  this  induces  a  truncation  on  the  nunber  of  failures 


r0  •  Vdj  "  (T^10*  8"  *  Figure  (6).J 
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T tone  if  the  randon  fanotlon  representing  the  lif«  teat  it  such  tbit  the  x^'tb  failure 

occurs  before  total  Ilfs  ▼  raaebaa  CO  )  ,  then  om  know*  that  If  ana 

continues  the  tast  until  total  Ufa  ?_  ,  than  tha  saaple  randan  function  nust 

o 

croaa  either  ■  or  Cl  and  in  aithar  eaaa  va  would  reject  H  •  Banca  If 

o 

the  r  *th  failure  occurs  bafore  total  life  V  one  can  step  with  tha  rejection 
o  o 

of  H  . 
o 

Suppose  now  that  one  preassigns  both  tha  nuaber  of  failures  rQ  and  total 

life  V  and  truncates  tha  sequential  life  tast  at  V  «  V  and  r  ■  r  .  Than 
o  o  o 

frae  the  foregoing  one  can  iapoae  the  following  equivalent truncations: 


(1)  For  r  truncate  at  r*  ■  nin 


-  Vo(^-^)/1°*  y  •  “1,“  [w*l 


(2)  For  V  truncate  at  V 


r  r°  io# 

’*  ■  Bin  |?o  |  -y -  jj^1" 

l  ar,  ~  y 

1  o 


I V  8  ro  ’ 


Clearly  r*  and  V*  will  met  the  condition 


[  See  Figure  (7)«  ] 

The  truncation  rule  is  as  follows:  if  the  sequential  probability  ratio  tast 
given  by  (31)  does  not  lead  to  a  decision  for  ?(t)  <  V#  and  r  <  r*  (i.e.,  if 
neither  AB  nor  FO  are  intersected  or  crossed)  ,  than  if  the  saaple  randosi 
function  associated  with  the  life  test  hits  the  boundary  V  ■  V*  (reaches  K) 
before  reaching  r  «  r*  (CO)  ,  accept  1Q  .  If,  however,  the  sample  randan 
function  hits  the  boundary  r  -  r*  (CO)  before  reaching  ?  -  ?*  (BC)  , 
reject  IQ  (accept  H^) 


Truncating  the  sequential  lift  test  at  failure  ouster  r*  and  at  total 
Ufa  V*  ,  i.e.,  accepting  Hq  ,  If  the  eaaple  random  function  acsoeiated  with 
the  life  test  eeets  AB  or  1C  tefort  crossing  VO  or  eseting  00  ,  and 
rejecting  Hq  if  the  Maple  randoa  function  crosses  VO  or  Mete  CO  before 
nesting  AB  or  K  ,  will  change  the  Type  X  and  Type  IX  err oa.  Let  o(r*  ,  V*) 
and  P(r*  ,  V#)  be  the  Type  I  and  Type  XX  errors,  respectively,  associated 
with  the  truncated  test .  It  is  clear  froa  what  we  have  said  above  that  O.C  . 
curves  associated  with  truncating  at  r  •  r*  ,  V  -  ,  rhere  V*  «  er*  coincide 

with  those  based  cm  truncation  at  r  ■  r*  or  V  ■  V*  .  Consequently  o(r*,V*) 

■  ®(r*)  and  P(r#,  V*)  -  P(V#)  .  Upper  bounds  given  previously  for  a(r*) 
end  P(V*)  are  automatically  upper  bounds  for  a(r*,  V*)  and  0(r*,V*)  . 
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Probability  of  termination  of  the  —  guential  ljf»  test 


at  preassigned  valuta  of  rQ  and  Vfl 


It  ia  interesting  to  aak  tba  question:  what  ia  the  probability  that  the 
sequential  life  teat  will  terminate  with  a  number  of  failuree  leaa  than  or 
equal  to  soam  preaaaigned  number,  rQ  ,  or  after  total  life  leaa  than  or  equal 
to  some  preaaaigned  value  f  Uaing  eooalderationa  analogous  to  those  in 

Maid's  bock  on  Sequential  Analysis,  pp.  58-60,  ve  can  state  the  following 
results  which  give  lower  bounds  for  the  probability  that  the  sequential  procedure 
will  terminate  with  a  number  of  failures  r  <  rQ  for  the  two  values  0  ■  0Q 
and  9  m  Pj, 

Consider  the  question  of  evaluating  the  probability  that  the  sequential 
life  test  terminates  with  a  number  of  failures  <  rQ  •  Then  using  consider¬ 
ations  like  those  in  Wald,  we  can  assert  that 

(2H.1)  Pr(r  <  r  |  9  -  0Q)  >  Pr(rQ  log  ^  -  (^-  -~)  V(xf  )  <  log  B  |  9  -  9Q) 

1  1  o  o 

-  Pr[  7?  (2 rQ)  >  jfj  (ro  log  k  —  log  B)  J  , 
since  2  V(xr  )  j  9Q  is  distributed  as  (2r0)  under  BQ  . 


And  similarly 


(2B.2)  Pr(r  <  r0  I  9  -  #x)  >  Pr  £rQ  log?S-(A.-  i)  f(x^)  >  log  A  |  B  -  dj 

•  Pr('3f  (2r0)  <  $  (rQ  iog  k  -log  A  )) 
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since  2  V(*r  )  /  is  distributed  u  96  (2rQ)  under  . 

In  a  similar  way  w  eaa  evaluate  the  probability  that  the  sequential  life 

test  teralnates  at  total  life  V(t)  <  ?  .We  can  give  the  following  lower 

bounds  for  the  probability  that  the  sequential  procedure  will  terminate  at 

total  life  V(t)  <  V  for  the  two  values  8  «  0  and  0  •  0.  : 

—  o  0  1 


(2H.3)  Pr(v(t)  <  VQ  |  0  -  80)  >  Pr(r  log  ^  -  (i  -  ?o  <  log  B  |  0  -  0Q) 

1  1  o 


^  J  p(r;  ~  )  ,  where  mq  «  [(k-1)  ~  +  log  b]  j  log  k 


—  -  o 


Since  under  the  hypothesis  that  HQ:  0  •  dQ  is  true,  we  are  observing  a 

*  ..  _ 

Poisson  process  with  rate  parameter  \q  ■  l/0Q  for  a  length  of  tins  Vq 


(2H.k)  Pr(V(t)  <  VQ  j  0  -  9J  >  Pr(r  log  -  (i  ~i-)  Vq  >  log  A  |  0  .  0^ 

1  1  o 

p(r;  ~  )  ,  where  ^  «  [(k-1)  ~  ♦  log  a|  j  log  k  , 


since  under  the  hypothesis  that  0-0^  is  true,  we  are  Observing  a 
Poisson  process  with  rate  parameter  -  1/0 ^  for  a  length  of  tine  Vq  . 

We  now  give  a  numerical  example  to  illustrate  what  we  have  Just  discussed. 
Consider  the  problem  of  testing  BQ:  0  ■  0Q  against  H^j  •  «  0^ ,  with  o  .  0  . 

•  .0^  and  k  »  0^0^  •  3  •  Let  us  compute  lower  bounds  for  Fr(r  <  rQ  |  0  ■  0q) 
and  Pr(r  <  rQ  |  0  ■  0^)  for  rQ  ■  10  ,  15  ,  20  ,  25  ,  30  and  lower  bounds 
for  Pr(V(t)  <  ?6  I  0  •  0X)  for  V0foQ  •  3  *  6  ,  9  ,  1 2  ,  and  15  .  We  first 
compute  Pr(  r  <^ro|  0^)  ,  i  ■  1  ,  2  for  this  example: 


If  rQ  m  10  ,  than 

Pr(r  <  10  |  0  •  eQ)  >  Pr(*jC2  (20)  >  10.99  ♦  2.94) 

-  Pr(%2  (20)  >  13-93)  -  .83 
and 

Pr(r  <  10  |  9  .  flA)  >  Prf'XJ  (20)  <  3(10.99  -  2-94)  ) 
.  .  PrtfC2  (20)  <  24.15)  -  .77  . 

If  rQ  ■  15  *  than 

Pr(r  <  15  i  0  -  eQ)  >  P rtf  (30)  >  19^3  )  -  .94 

and 

P r(r  <  15  I  0  -  0j)  >  PrCjS  (30)  <  40^65  )  -  .90  • 

If  rQ  m  20  ,  than 

Pr(r  <  20  |  0  -  9q)  >  Pr(X2  (40)  >  24-92  )  «  .97 

and 

Pr(r  <  20  |  0  m  9X)  >  (40)  <  57-12  )  -  .96  . 

If  rQ  ■  25  ,  than 

Pr(r  <  25  I  0  -  0Q)  >  P>r(%2  (50)  >  30.42  )  .  .99 


Pr(r  <  25  I  9  m  0t)  >  p (50)  <  73.62  )  -  .96  . 
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(  If  ro  •  30  ,  then 

Fr(r  <  30  |  9  -  90)  >  PrtK?  (60)  >  33-91  )  -  .993 

and 

Pr(r  <  30  |  9  -  9X)  >  Pr( 'Xf  (60)  <  90.13  )  -  -998  • 

Blailarly  we  compute  Pr(V(t)  <  VQ  |  9  -  0^  ,  i  «  1.2  ,  for  this  example. 

vo 

If  g~  •  3  1  then 


Bence 


and 


Pr(V(t)  <  V  |  9  -  9  )  > 


p(ri  3)  »  »*»•*• 


o  rr  _ 

—  -  o 


»o  .  kk-1)  f-  ♦  log  B  j  /log  k  .. 

*■  o 


J  ^lo g  Hr  js  -  2-9*j  j 


1.1  -  2.8  . 


Pr(V(t)  <  V  |  9  -  9  )  >  ^  p(r*  3)  -  *2 


fs6 


Pr(V(t)  <  V0  i  ®  •  ®x)  >  2^  '  wh*rt 

^  .  j(k-l)  ^  ♦  log  aJ  /  log  k  .  £  ♦  2.9* J  y 

Bence  Pr(V(t)  <  ?Q  |  9  -  9X)  >  p(r;  9)  -  -53  . 


1.1  -  8.1  . 


i 


If  a—  ■  6  ,  then 
9 0 
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?n 


and 


Pr(V(t)  £V0I9,90)^  P<*l  6)  . 

* 

<  V„  I  *  -  »!>  >  £  P(r;l8) 

r=l3 


If  ^*9,  than 


Pr(v(t)  <  v0  I  e  -  eo)  >  p(p,  9) 


and 


Pr(V(t)  <  V  |  0 


P(r;  27)  ■ 


If  ^  «  12  ,  than 
®o 


Pr(V(t)  <  ?0  I  0  -  0o)  >  p(r;  12)  - 


and 


Pr(V(t)  <  VQ  |  0  .  91)  >  Sp(r  i  36)  • 


If  sr*  ■  15  #  then 
9 o 


Pr(V(t)  <  V0  |  9  -  0Q)  >  2_  P(r?  15)  - 


>-£ 


and 


.85 


.908  . 


.926 


•  956  . 


•979 


.986  . 


•989 


Mv(t)  <  v0  I  «  -  ©j)  >  X.  Pto  *5)  -  -996  . 
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Tables  for  ths  probability  that  a  sequential  life  test  will  teralnate  by 
a  preassigned  rQ  and  VQ  are  being  coaputed  for  other  values  of  a  t  0  end 
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Appendix  81 

Ppper  and  lover  bound*  tor  L(0)  and  Efl(r) 

Am  formulae  for  L(t)  and  K0(r)  ,  given  by  (37)  and  (38),  respectively, 
are  approximations  to  tbs  actual  L(0)  and  actual  Sg(r)  arising  from  the  use 
of  tbs  sequential  rule  specified  by  tbe  inequalities  (31).  The  question  arises 
as  to  bow  good  these  approximations  are.  A  modification  of  the  results  of  Wald  on 
bounds  for  the  O.C.  and  ASH  curves  In  the  binomial  case  and  of  results  of 
Herbach  on  the  discrete  Ralston  yields  the  following  bounds  on  the  actual  L(0) 
and  S0(r)  s 


where  the  upper  inequality  sipis  hold  for  0  <  s  and  the  lover  Inequality  signs 
hold  for  0  >  8  . 

One  unpleasant  feature  of  the  bounds  given  in  (11)  le  that  they  Involve 
L(0)  ,  which  is  unknown.  However,  this  matters  little  lu  actual  practice  because 
the  limits  on  L(0)  given  by  (1)  ere  quite  close  together  for  the  range  of 
values  of  k  sad  (0,0)  covered  in  Table  10.  Thus,  for  example,  for 

k  -  9jBx  -  3  ,  a  -  0  •  .05  ,  A  -  (1  -  P)/a  -  19  ,  b  -  0/(1  -  0)  -  1/19  , 
ve  get  that  .95  <  L(0O)  <  .963  and  .00  <  L(0j)  <  .000  .  The  upper  and  lover 
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bounds  for  S0(r)  given  by  (il)  an  elm  together  for  6  ■  and 
ooaperatlvely  far  apart  for  6  •  8^  .  Thus  for  tha  ease  k  •  3  and 
o  ■  0  m  .05  ,  the  difference  between  the  typer  and  lower  bounds  Is  <  .06 
for  9  m  6  and  is  about  2.5  for  8  »  8,  . 

The  left  aide  of  (il)  is  the  approximate  fonaula  (38)  for  I0(r)  except 

that  the  L(8)  in  (il)  refers  to  the  exact  value  and  the  L(8)  in  (38)  is 

given  by  the  approximation  (37).  In  view  of  the  preceding  paragraph,  the  values 

of  Kfl  (r)  given  in  Table  10  are  very  close  to  the  correct  values,  while  the 
"o 

values  of  E„  (r)  are  essentially  lower  bounds  for  the  correct  value.  He 
cannot  say  more  unless  we  go  through  more  extensive  calculations  of  the  sort 
to  be  described  in  Appendix  2J. 


Appendix  2J 

Some  exact  calculations  of  L(6)  wad  B0(r) 

Wald  pointed  out  that  in  order  to  have  a  test  of  exactly  strength  (ft ,fi  )  , 

f  f 

the  A  and  B  in  (31)  should  be  replaced  by  A  and  B  ,  where 
A*  <  A  •  (l  -  p)/a  and  B*  >  B  »  8/(1  -  a)  .  In  the  present  case,  with 
information  available  continuously  In  time,  B  -  B  ■  0/(1  -  a)  since  the 
acceptance  of  Hq  Involves  no  excess  over  the  boundary.  However,  acceptance 
of  does.  In  general,  entail  a  positive  excess  over  the  boundary,  and 
all  we  can  say  Initially  about  A*  is  that  it  should  lie  between  A  8j /0O 

and  A  •  Thus  using  A  »  (1  -  0)/a  Instead  of  A  is  an  approximation . 

The  approximate  test  based  on  using  A  and  B  is  suitable  for  all 
practical  purposes,  since  one  consequence  of  the  inequalities  in  Appendix  21 
is  that  the  strength  (a'  ,8')  is  such  that  a'  <  a  ,  ^0'  <0/(1  •  a)  . 

Since  O  and  0  are  generally  snail  (<  .10  say)  a  procedure  based  00  A  and  B 
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provides  essentially  the  mm  protect  ion  against  err  ora  of  tha  first  and 
saoond  kind  as  doaa  tha  taat  baaad  on  using  A*  and  B*  •  However,  tha 
nsa  of  A  rathar  than  A*  la  (3D  will  antall  a  aaall  inorsass  in  £Q(r) 
particularly  for  0  <  a  • 

A a  a  practical  matter,  oua  would  usually  be  content  with  a  taat  baaed 
on  (3D  which  uses  A  and  S  •  As  a  matter  of  fact,  this  is  what  la  dona 
all  tha  tins  by  people  faced  with  a  practical  decision  problem,  for  moat 
sequential  problems,  tha  problem  of  finding  tha  A*  and  B*  which  will 
give  exactly  strength  ( OC .  /3 )  has  not  bean  solved.  One  has  to  rely, 
in  such  oases,  on  tha  results  of  Held  which  indicate  that  tha  errors  in¬ 
volved  in  using  A,  B,  and  approximate  formulae  for  L(0)  and  EQ(r) 
are  "reasonably"  small* 

In  tha  problea  at  hand  we  know,  in  view  of  the  continuous  availa¬ 
bility  of  information,  that  B*  a  B  a  /*/(l  -  <f )  •  Furthermore,  formulae 
are  available  for  computing  A*  and  for  computing  O.C.  and  ASH  curves 
exactly.  The  formulae  for  accomplishing  these  tasks  are  available  in 
papers  by  Burman  and  by  Dvoretsky,  Kiefer,  and  Wolfowits.  While  the 
computational  labor  involved  in  any  special  oase  is  exceedingly  heavy, 
the  results  of  suoh  computations  do  throw  some  light  on  how  exact  O.C.  and 
ABN  curves  compare  with  those  computed  by  using  approximations. 

Formulae  (4.17)  and  (4.23)  in  the  Dvoretsky,  Kiefer,  and  Wolfowits 
paper  (similar  formulae  are  given  in  Barman's  paper,  p.  102)  were  used 
to  oompute 

(i)  the  exact  O.C.  and  Eg(r)  curves  for  the  sequential  rule  (31) 
with  B  a/ff/(l  -  oC)  and  A  a  (1  -  )/<*  •  This  was  done  for  the  case 
k  a  0 a  3  and  <Ca  (3  •  .03,  and 
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(ii)  A*  (whara  A  ®1/®0  £  i'  5  *)  auch  that  tha  daeision  rula 

^/(l  -  <c)  <  <#0/«1)r  axp£  -  (1/»1  -  l/®0)  ▼  (t)3  <  A* 

haa  aa  O.C.  ourra  for  which  L(©0)  ■  1  -  OC  and  1(0^)  a  0  axaetly,  and 
thaa  to  coaputa  SQ(r)  for  tha  (B,A*)  rula.  This  was  dona  for  tha 


oasaa  A  * 

0  -  .05 

and  k 

.  eo/»i  , 

■  3/2,2 

t  and  9, 

and  also 

and  k  •  3 

• 

Tha  rasult  of  (i)  was 

K«0) 

■  .988, 

L(.)  . 

.529, 

!(©!> 

«  .051  . 

Bg  (r) 

0 

-  9.09, 

B,(r)  . 

8.10, 

(r) 

*1 

*  7.00  . 

Computation  (ii)  gara 


OC  «/S 

k 

A* 

(r) 

O 

E^Cr) 

(r) 

W1 

3 

13.25 

2.94 

7.22 

6.21 

.05 

i  2  „ 

15.1 

8.64 

18.0 

13.8 

3/2 

16.6 

27.9 

52.8 

96.8 

.01 

3 

68.9 

5.00 

17.5 

10.5 

Baaring  la  mind  that  tha  eoaputatioaa  wara  carriad  through  only  in 
a  small  nuabar  of  caaas,  ona  can  maka  thraa  obsarrations: 

(a)  For  tha  easa  k  ■  3  and  Am  fi  m  .05,  tha  uea  of  B  *  1/19 

t  • 

and  A  *  19  rasults  in  gatting  oC  •  .092  and  /d  >  .051  as  eoaparad 
with  dC  •  /3  »  ,05  whan  ona  usas  B*  *  B  ■  1/19  and  A*  ■  19*25.  Also, 
Eg(r)  ia  inoraasad  by  .09,  .88,  and  *79  at  ©  »  ©Q,  a,  ©1  raapaotiyaly. 
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(b)  Of  more  interest  is  the  fact  that  tha  axact  values  of  fig(r) 
for  tha  (B,A*)  rula  practically  coinoide  with  tha  approximate  values  of 
Bg(r)  computed  for  tha  (B,A)  rula  by  formulae  (58)  and  (59)  and  given 
in  Table  10. 

(e)  In  tha  range  of  values  of  k  a  0^®^  and  of  oC  and  covered 

by  Tabla  10,  a  good  guass  at  tha  value  of  A*  la  tha  value  A**  lying 

* 

midway  batwaan  A  *nA  A/k  ,  tha  uppar  and  lowar  limits  on  A  .  This 
naans  that  A**  a  (k+1)  A/2k  .  On  tha  basis  of  our  liaitad  calculations 
wa  conjecture  that  in  tha  ranga  of  values  oovarad  in  Tabla  1,  a  sequential 
dacision  rula  based  on  (31)  with  A  replaced  by  A**  will  hare  almost 
exactly  strength  ( <£,  fi)  .  The  values  of  Eft(r)  associated  with  a 
(B,A**)  rula  will  be  given  to  a  close  approximation  by  (38). 


Appendix  2K 

An  approximate  formula  for  £Q(t)  in  tha  nonraplacamant  case 

A  useful  approximation  to  Bg(t)  in  tha  nonreplaooment  ease  is  given 
by  Ee(t)'v'Q  log  (n/tn  -  E^r)])  •  This  approximation  is  obtained  by 
replacing  EgU^)  in  (42)  by  its  approximation  0  log  (n/ILn  -  k]f). 

Thus  (42)  becomes  (43) 

*.(*>  *.  [l»«  (r?i)]~®  1o«  (.-^w)  • 

This  approximation  has  bean  tasted  numerically  by  calculations  on 
trunoated  nonraplacamant  dacision  procedures,  where  tha  exact  valuas  of 
S0(t)  can  be  oosputad  and  compared  with  tha  suggested  approximation.  Tha 

agreement  is  dose. 


r 


2.61 

. 

REFERENCES  FOR  SECTION  4 

Q.  E.  Albert.  "Accurate  sequential  teste  on  the  seen  of  an  exponential 
distribution ,"  Ann.  Math.  Stat. ,  Vol.  27*  pp.  460-470,  1956. 

J.  P.  Bureau,  "Sequential  saapling  formulae  for  a  binomial  population." 

J.  Roy.  Stat.  Soo.,  Ser.  B.  Vol.  8,  pp.  96-105.  1946. 

J.L.  Doob,  Stochastic  Processes,  John  l/iley  and  Sons,  1955* 

A.  Dvoretsky,  J.  Kiefer,  and  J.  Wolfowitn,  "Sequential  decision  probleas 
for  processes  with  continuous  time  parameter.  Testing  hypotheses," 
Ann.  Hath.  Stat.,  Vol.  24,  pp.  254-264,  1955* 

J.  Kiefer  and  J.  Wolfowitn,  "Sequential  tests  of  hypotheses  about  the 

seen  occurrence  time  of  a  continuous  paraaotor  Poisson  prooess," 
Naval  Research  Logistics  Quarterly  5.  PP*  205-214,  1956. 

B.  Epstein  and  11.  Sobel,  "Life  testing,"  J.  Aaer.  Stat.  Assoe.,  Vol.  48, 
pp.  486-502,  1953. 

B.  Epstein,  "Statistical  probleas  in  life  testing,"  Proceedings  of  the 
Seventh  Annual  Convention  of  the  Aaeriean  Society  for  Quality 
Control,  pp.  385-398,  1953. 

B.  Epstein,  "Truncated  life  tests  in  the  exponential  case,"  Ann.  Math. 
Stat.,  Vol.  25,  pp.  553-564,  1954. 

B.  Epstein  and  M.  Sobel,  "Sequential  life  tests  in  the  exponential  ease," 
Ann.  Math.  Stat.,  Vol.  26,  pp.  82-93,  1955. 

L.  H.  Berbaeh,  "Bounds  for  some  functions  used  in  sequentially  testing  the 
seen  of  a  Poisson  distribution,"  Ann.  Hath.  Stat.,  Vol.  19, 
pp.  400-405,  1948 

A.  Wald,  Sequential  Analysis,  John  Wiley  and  Sons,  1947. 


t 


Tmx  i 


Valuta  of  Xl^(2r)/2r 


rNTM  .01 

.05 

.10 

.25 

.56 

1 

.010 

.052 

.106 

.288 

.693 

2 

.074 

.178 

.266 

.481 

0839 

3 

<,145 

.272 

.367 

.576 

.891 

4 

.206 

.342 

.436 

<,634 

.918 

5 

.256 

.394 

.486 

.674 

.934 

6 

.298 

.435 

.525 

.703 

<945 

7 

.333 

.469 

.556 

.726 

.953 

8 

.363 

.498 

.582 

o?44 

.959 

9 

.390 

.322 

n  604 

,,760 

.963 

10 

.413 

.543 

,622 

.773 

.967 

15 

.498 

«6l6 

.687 

.816 

.978 

20 

.554 

0663 

.7 26. 

.84  2 

.983 

25 

o594 

0695 

.,754 

.859 

.987 

30 

.625 

.720 

.774 

.372 

.989 

40 

.669 

.755 

.803 

.839 

,992 

50 

.701 

.779 

.824 

.901 

.993 

75 

.751 

.818 

.855 

.920 

.996 

100 

.782 

.841 

.874 

.931 

.997 

goto;  This  table  la  used  in  the  following  way: 

Aoeept  0  •  0Q  if  >  ®oXi-^2r^/2r  and  reject  otherwise. 

Suppose,  for  exanple,  that  we  want  to  discontinue  a  life  test  after 
r  *  3  failures  have  occurred  and  that  we  want  the  life  test  to  be  such  that 
a  lot  having  oe an  life  1000  hours  is  accepted  with  probability  »90<» 

Using  formula  (7)  and  table  (1)  the  region  of  acceptance  is  given  by 

*3|a  >  (1C00) (..I*?)  -  567. 

Zn  words,  one  places  n  iteas  on  life  test  and  stops  testing  after 

3  iteas  have  failed.  One  then  ooapufees  ,  an  estimate  of  the  mean 

p.n 


2.83 


life  after  3  failures,  using  formula  2  in  the  non-roplacemeat  case  and 
formula  3  1a  the  replacement  case*  One  accepts  the  lot  if  6  >  36?  and 
rejects  otherwise .  Suppone,  for  example,  that  we  place  a  *  10  items 
on  test,  do  not  replace  items  as  they  fail,  and  that  the  first  3  failure 

times  are  50,  125,  250,  In  this  case  &3el0  3  B  Q \21  ,  725. 

Since  725  >  3^7,  *e  accept  the  lot. 


Number  of  failures 


I 


2.84 
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Values  of  0  accepted  with  probability  pf  when  0  ■  1  is  accepted  with 


probability  .99.  Rule  of  action  is  to  accept  if  $  >  (2r)/2r 

- 1- -  - • _ m®  A*99 


AL 

■■Lag- 

_ *i5_ 

*30 

»75 

. . .50 

.  .25 

.10 

.05 

.01 

1  -- 

1.0 

.194 

.095 

.035 

.014 

.007 

,004 

.003 

o002 

2 

1.0 

.418 

.279 

.154 

.088 

.055 

.038 

.031 

.022 

3 

1.0 

«533 

.396 

.252 

.163 

.111 

.082 

.069 

.052 

4 

1.0 

.602 

.472 

.325 

.224 

.161 

>123 

.106 

.082 

5 

1.0 

.649 

.526 

.380 

.274 

.204 

.160 

.140 

.110 

6 

1.0 

0683 

o$66 

.423 

.315 

.241 

>193 

.170 

.136 

7 

1.0 

.709 

.598 

.458 

.349 

.272 

.221 

.197 

ol60 

8 

1.0 

.730 

.624 

.488 

.379 

.300 

.247 

.221 

.182 

9 

1.0 

.747 

.646 

.513 

.405 

.325 

.270 

.243 

.202 

10 

1.0 

.761 

.664 

.535 

,427 

.347 

.291 

.263 

.220 

15 

1.0 

.809 

.726 

.611 

.510 

.430 

.371 

.342 

.294 

20 

1.0 

.836 

.763 

.658 

.563 

.486 

.428 

.398 

.348 

25 

1.0 

0855 

.788 

.692 

.602 

.527 

.470 

.440 

=  390 

30 

1.0 

.868 

.807 

.717 

.632 

.560 

.504 

.474 

.424 

40 

1.0 

.837 

.833 

.753 

v  675 

.608 

.554 

.526 

.477 

50 

1.0 

0899 

.851 

.777 

.705 

.642 

.591 

563 

.516 

75 

1.0 

.918 

.878 

.317 

.754 

.699 

.653 

.627 

.583 

100 

..hi  - 

.930 

.895 

.840 

.785 

.734 

.692 

.669 

.627 

N0te:  Tables  2(a)  through  2(e)  give  O.C.  curves  associated  with  test  pro¬ 
cedures  of  the  form:  Accept  if  >  ^2_^(2r)/2r,  Time  units  are 

chosen  in  such  a  way  that  the  probability  of  accepting  0  »  1  is 
1  •of  . 

E3WP?r9a  at-Mg.  ?*  W 9  2  a,, 

from  table  (1)  we  know  that  the  acceptance  rule  4-  >  1000  (*256)3256 

/in 

will  lead  to  the  acceptance  of  a  lot  with  mean  life  0  3  1000,  with 
probability  .99.  From  table  2(a)  we  can  say  that  a  lot  with  mean  life 
526  will  be  accepted  with  probability  .90,  that  a  lot  with  mean  life 
274  will  be  accepted  with  probability  .50,  and  a  lot  with  mean  Ilfs 
160  will  be  accepted  with  probability  .10. 


i 

\ 

« 

1 

i 

! 

\ 

i 

! 


i 

s 

I 

1 


i 

I 

1 

! 

1 

I 

I 
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TABU  2(b) 


Valuta  of  0  accepted  with  probability  p,  whan  0  ■  1  ia  acceptad  with 
probability  .95.  Sul©  of  action  ia  to  aocapt  if  >  ^^j(2r)/2r. 


.99 

«95 

.90 

.75 

.50 

.25 

.10 

.05 

.01 

1 

5.150 

1.0 

.488 

.179 

.074 

.037 

.02 2 

.017 

.011 

2 

2.394 

1.0 

.668 

.370 

.212 

.132 

.091 

.075- 

.054 

3 

1.875 

1.0 

.742 

.473 

.306 

.209 

.154 

.130 

.097 

4 

1.66o 

1.0 

.783 

.539 

.372 

.26? 

.205- 

.176 

.136 

5 

1.540 

1.0 

.810 

.585 

.42  2 

.314 

.246 

.215 

.170 

6 

1.463 

1.0 

.829 

.619 

.461 

.352 

.282 

.249 

.199 

7 

1.410 

1.0 

.844 

.646 

.493 

.384 

.312 

.277 

.225 

8 

1.370 

1.0 

.855 

.668 

.519 

.411 

.338 

.303 

.249 

9 

1.339 

1.0 

.864 

.687 

.542 

.435 

.361 

.325 

.270 

10 

1.314 

1.0 

.872 

.702 

.561 

.455 

.382 

•  345 

.289 

15 

1.257 

1.0 

.898 

i755 

.630 

.531 

.459 

.422 

.363 

20 

1.196 

1.0 

.912 

.788 

.674 

.581 

.512 

.475 

.416 

25 

1.170 

1.0 

.922 

.810 

.705 

0617 

.550 

.515 

.456 

30 

1.152 

1.0 

.930 

.C26 

.728 

.645 

.581 

.546 

.489 

40 

1.128 

1.0 

.940 

.849 

.761 

.685 

.625 

.595 

.538 

50 

1.112 

1.0 

.946 

.865 

.785 

.714 

.658 

.627 

.574 

75 

1.089 

1.0 

.956 

.889 

.822 

.761 

.711 

.683 

.635 

100 

1.076 

1.0 

.962 

.904 

.844 

.790 

.745 

.719 

0675 

ExaaDla:  If  the  life 

teat  ia 

diacontinued  after  r  = 

1  failure  occura,  and 

if  a  lot  with 

aean  life  0  ■  1 

ia  accepted  with  probability 

.95,  than 

a  lot 

with  naan  Ufa  0**  .074  ia  acceptad  with  probability  .50.  Hula  of  action 


is:  Accept  if  ‘S.  -  >  .052. 

l(n 

If  the  life  teat  ia  diacontinued  after  r  ■  5  failurea  occur,  and  if 
a  lot  with  aean  life  0  *  t  ia  acceptad  with  probability  .95.  than  a  lot 
with  aean  life  0  »  .422  ia  acceptad  with  probability  .50.  Rule  of  action 
ia:  Accept  if  >  *394. 


Number  of  failures 
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Table  2(c) 

Values  of  ©  accepted  with  probability  p,  when  8  =  1  is  accepted  with 
probability  „90»  Rule  of  action  is  to  accept  if  ©r  Q  >  X  g0(2r)/2.r. 


X’ 

-99 

>95 

,90 

.75 

-50 

.25 

.10 

.05 

— 

.01 

1 

10.550 

2.049 

1.0 

.367 

.152 

,0?6 

.  046 

-035 

.023 

2 

3.582 

1.496 

1.0 

.553 

.317 

.198 

.137 

,112 

.080 

3 

2.528 

1-348 

1.0 

.638 

.412 

.281 

.207 

.175 

.131 

4 

2  c  120 

1.277 

1.0 

.688 

.475 

.342 

.261 

,,225 

.174 

5 

1.902 

1-235 

1.0 

,722 

.521 

.388 

.304 

,266 

.210 

6 

1.765 

1.20  6 

1.0 

.747 

-556 

.  425 

.340 

„300 

.240 

7 

1.672 

1.186 

1.0 

,766 

.584 

.455 

-370 

-329 

.267 

8 

1.602 

1.170 

1.0 

.782 

.607 

.  43l 

-396 

.  354 

.291 

9 

1-549 

1.157 

1.0 

795 

.627 

-503 

.418 

-376 

.312 

10 

1.506 

1-14? 

1.0 

805 

.643 

.52  2 

.438 

.396 

.351 

15 

1.377 

1.114 

1.0 

,842 

.702 

-592 

.512 

,471 

.405 

20 

1.311 

1.096 

1.0 

.  863 

.739 

.637 

..561 

521 

.  456 

25 

1.269 

1.084 

1.0 

.878 

.764 

.669 

-  597 

.558 

•  495 

30 

1.239 

1.076 

1.0 

..888 

.783 

.694 

,624 

-587 

.526 

40 

1.201 

1.064 

1.0 

-903 

.310 

•  729 

.666 

.631 

-572 

50 

1.175 

1.057 

1.0 

.914 

.829 

.755 

-695 

.662 

.606 

75 

1.139 

1.046 

1.0 

930 

.359 

.  795 

.743 

.714 

.664 

100 

1.118 

1.039 

1.0 

939 

-377 

.820 

.774 

-747 

.701 

Example :  If  the  life  teal  is  discontinued  after  r  =  10  failures  occur,  and 

if  a  lot  with  mean  life  0=1  is  accepted  with  probability  ,90,  then 

* 

a  lot  with  mean  life  8  =  .643  is  accepted  with  probability  .50,  and 

•  * 

a  lot  with  neon  life  ©  =  .396  is  accepted  with  probability  .05.  •• 

Rule  of  action  is:  Accept  if  n  >  ,622c 
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VU.ua*  of  $  accepted  with  probability  p,  when  0-1  ia  accaptad  with 
probability  .75.  Rula  of  action  ia  to  accept  if  “9  „  >X  2-.(2r)/2 r 


5  .50 


28.750 

5.582 

2.725 

1.0 

.415 

.207 

.125 

.090 

.062 

6.475 

2.705 

1.807 

1.0 

.573 

.357 

.247 

.203 

.145 

3.902 

2.113 

1.568 

1.0 

,646 

.441 

.325 

.274 

.206 

3.C61  . 

1.855 

1.453 

.  1.0 

.690 

.496 

.380 

.327 

.252 

2.634 

1.710 

1.385 

1.0 

.721 

•537 

.421 

.308 

.290 

2.303 

I.0I5' 

1.339 

1.0 

.744 

.568 

.455 

.401 

.322 

2.181 

1.547 

1.305 

1.0 

.762 

.594 

.483 

-429 

.349 

2.050 

1.496 

1.279 

1.0 

•777 

.615 

.506 

.453 

.372 

1.949 

1.456 

1.259 

1.0 

.789 

.633 

•526 

.474 

•393 

1.871 

1.424 

1.242 

1.0 

•799 

.648 

.544 

.492 

.411 

1.637 

1.324 

1.188 

1.0 

.834 

.703 

.608 

•559 

.481 

1.519 

1.270 

1.159 

1.0 

.856 

.738 

.650 

.604 

.528 

1.446 

1.235 

1.139 

1.0 

.870 

.762 

.680 

.636 

.564 

1-395 

1.211 

1.126 

1,0 

.881 

.781 

.703 

.661 

.592 

1.329 

1.178 

1.107 

1.0 

.897 

.807 

•737 

.698 

.633 

1.266 

1.157 

1.094 

1.0 

•  907 

.826 

.761 

.725 

.664 

1.225 

1.125 

1.076 

1.0 

.924 

.855 

.800 

.768 

.714 

1.190 

1.106 

1.065 

1.0 

•934 

.874 

.824 

.796 

.746 

If  the  life  teat  ia  discontinued  after  r  -  10  failure*  occur,  and 
if  a  lot  with  naan  life  0-1  ia  accepted  with  probability  .75,  than 
a  lot  with  aean  life  0*  «  1.424  la  accepted  with  probability  .95  and  a 
lot  with  Man  life  0**  »  .544  ia  accepted  with  probability  .10. 


lot  with  aean  life  0  »  .544  ia  accepted  with  probability  .10. 

A 

Rule  of  action  la:  Accept  if  0^o  >  .773. 
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TiBLE  2(c) 


Values  of  0  accepted  with  probability  p,  whan  9  -  1  is  accepted  with  probability 
6 50 •  Rule  of  action  is  to  accept  if  9r  n  >  / 2r o 


30 

1.583 

1.374 

1.277 

1.135 

4o 

1.482 

1.314 

1.234 

1.115 

50 

1.418 

1.275 

1.206 

1.102 

75 

1.325 

1.217 

1.164 

1.082 

100 

1.274 

1.185 

l.lliO 

1.071 

»50 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1,0 

1.0 

1,0 

1.0 

1.0 

1.0 

1.0 

1,0 

1.0 

1.0 

1.0 

1.0 


.25 

.500 

.623 

.682 

.719 

,744 

•764 

.779 

.792 

.802 

.812 

<843 

.862 

.876 

.886 

.900 

.910 

.926 

.935 


.10 

o301 
cli32 
,502 
.550 
.58 a 
.611 
,633 

.652 

.667 

068I 

.729 

.759 

.781 

.798 

.821 

,838 

.865 

.882 


o0$  .01 

.231  ol50 

.354  o253 

.425  .318 

,474  .366 

<510  .403 

.539  ,433 

,563  .458 

,583  .479 

,601  .498 

616  .515 

,670  .576 

,705  .618 

,731  .648 

,750  .671 

<779  -706 

.799  .731 

.832  .773 

.852  .799 


Example!  If  the  life  test  is  discontinued  after  r  *  10  failures  occur,  and  if 
a  lot  with  mean  life  9-1  is  accepted  with  probability  ,$),  then  a 
lot  with  mean  lfJTe  9*  *  1.782  is  accepted  with  probability  .95  and  a 
lot  with  mean  life  9  »  .681  is  accepted  with  probability  ,10. 

A 

Rule  of  action  isi  Accept  if  n  >  .967  . 


TABLE  3(a) 

Values  of  E(X  )/0  ,  where  X  is  the  Smallest  Value  in  a 
r,n  r,n 

Randan  Sample  of  Site  n  Drawn  from  a  Distribution  whose  Probability 
Density  Function  is  g  e”*^9  ,  x  ,  0  >  0  . 

In  the  table  r  =  l(l)n  and  n  «  1(1)20(5)30(10)100  . 


1 


1.0000 


0.5000 

1.5000 


0.3333 
0.83  3  3 
1.8333 


n 


0.2500 

0.2000 

0.5833 

0.1i500 

1.0833 

0.7833 

2.0833 

1.2833 

2,2833 

0.1667 

0.3667 

0.6167 

0.9500 

lJt500 

2.h500 


H Hi 

8 

9 

10 

n 

0oll429 

0.125c 

0,1111 

0.1000 

0.0909 

0.3095 

0.2679 

0.2361 

0.2111 

0,1909 

0.5095 

U\ 

• 

O 

0.3790 

0.3361 

0,3020 

0.7595 

0,63li5 

0.51'56 

0.1790 

0.1i270 

1.0929 

0.881j' 

o.71j56 

0.6J/56 

0,5699 

1.5929 

1.2175 

0.9956 

0.8i<56 

0.7365 

2.5929 

1.7179 

1.3290 

1.0956 

0.9365 

2.7175 

I.O290 

2.8290 

1^ ^290 

1.9290 

2.9290 

1.1865 

1.5199 

2.0199 

3.0199 

0.0833 

0.17l»2 

0.27l»2 

0.3853 

0.5l0h 

0.6532 

0.8199 

1.0199 

1.2699 

1.6032 

2.1032 

3.1032 


•7  * 


23 

14 

15 

16 

17 

18 

0.0769 

O.O714 

0.0667 

0.0625 

0.0588 

0.0556 

0.1603 

0.1484 

0.1381 

0.1292 

0.1213 

0.1144 

0.251 2 

0.2317 

0.2150 

0.2006 

0.1880 

0.1769 

0.3512 

0.3226 

0.2984 

0.2775 

0.2594 

0.2435 

0.4623 

0.4226 

0.3893 

0.3609 

0.3363 

0.3150 

0.5873 

0.5337 

0.4893 

0.4518 

O.4197 

0.3919 

0.7301 

0.65*7 

O.6OO4 

0.5518 

O.51O6 

0.4752 

0.8968 

0.8016 

0.7254 

0.6629 

0,6106 

0.5661 

1.0Q68 

0„9682 

0.8682 

0.7879 

0.7217 

0.6661 

1.3468 

1.1682 

1.0349 

0.9307 

0.8467 

0.7773 

1.6801 

1.4182 

1.2349 

1.0974 

0.9896 

0.9023 

2.1801 

1.7516 

1.4S49 

1.2974 

1.1562 

1.C451 

3.1801 

2.2516 

1.8182 

1.5474 

1.356 2 

1.2118 

3.2516 

2.3182 

1.8807 

1.6062 

1.4118 

3.3182 

2.3807 

1.9396 

1.6618 

2.3807 

2.4396 

1.9951 

3.4396 

2.4951 

3.4951 

r 

n  * 

19 

1 

0.0526 

7 

0.4445 

14 

1.2644 

2 

0.1082 

8 

0.5279 

15 

1.4644 

3 

0.1670 

9 

0,6188 

16 

1.7144 

4 

0.2295 

10 

0.7188 

17 

2-0477 

5 

0.2962 

11 

0.8299 

18 

2.5477 

6 

0.3676 

12 

0.9549 

19 

3.5477 

13 

1.0977 

p 

r 

r 

r 

1 

0,0400 

8 

0.3764 

14 

0.7961 

20 

1.5326 

2 

5 

to 

O 

o 

9 

0.4352 

15 

0.8870 

21 

1.7326 

3 

0,1251 

10 

0.4977 

16 

0.9870 

22 

1.9826 

4 

0.1706 

11 

0.5644 

17 

1.0981 

23 

2.3160 

5 

0.2182 

12 

0.6358 

18 

1.2231 

24 

2,8160 

6 

0.2682 

33 

C.7027 

19 

1.3660 

25 

3.8160 

7 

0.3209 

r 

r 

1 

0.0333 

11 

0.4472 

21 

1,1660 

2 

0.0678 

12 

0.4999 

22 

1  2771 

3 

0.1035 

13 

0.5554 

23 

1.4021 

4 

0.1406 

14 

0.6143 

24 

1,5450 

5 

0.1790 

15 

0.6768 

25 

1,7117 

6 

0.2190 

16 

0.7434 

26 

1,9117 

7 

0.2607 

17 

0.8149 

27 

2.1617 

8 

0.3042 

18 

0.8918 

28 

2.4950 

9 

0.3496 

19 

0.9751 

29 

2.9950 

10 

0.3972 

20 

1.0660 

30 

3.9950 

i  (Coot  > ) 
40 


r 

r 

r 

V 

1 

0.0290 

11 

0.3169 

a 

0.7308 

31 

1.4496 

2 

0.0906 

12 

0.3514 

22 

0.7834 

32 

1.5607 

3 

0.0770 

13 

0.3871. 

23 

0.8390 

33 

1.6857 

4 

0.1040 

14 

0.42a 

24 

0.8978 

34 

1.8285 

5 

0.1318 

15 

0.4626 

25 

0.9603 

35 

1.9952 

6 

0.1603 

16 

0.9026 

26 

1.0270 

36 

2.1952 

7 

0o1897 

17 

0.5443 

27 

1.0984 

37 

2.4452 

8 

0.2200 

18 

0.5877 

28 

1.1753 

38 

2.7785 

9 

0.2913 

19 

0.6332 

29 

1.2587 

39 

3.2785 

10 

0.2836 

20 

0.6808 

30 

1.3496 

40 

4.2785 

n  *  90 

r 

r 

r 

r 

1 

0.0200 

16 

0.3810 

31 

0.9515 

46 

2.U59 

2 

0.0404 

17 

0.U04 

32 

1.00a 

47 

2.6659 

3 

0.0612 

18 

0.4407 

33 

1.0597 

48 

2.9992 

4 

0.0825 

19 

0.4720 

34 

1.1185 

49 

3.4992 

5 

0.1043 

20 

0.9042 

35 

1.1810 

50 

4.4992 

6 

0.1265 

21 

0.5376 

36 

1.2476 

7 

0.1492 

22 

0.5720 

37 

1.3191 

8 

0ol725 

23 

0.6C77 

38 

1.3960 

9 

0.1963 

24 

0.6448 

39 

1.4793 

10 

0.2207 

25 

0.6832 

40 

1.5702 

11 

0.2457 

26 

0.7232 

a 

1.6702 

12 

0.2713 

27 

0.7649 

42 

1.7813 

13 

0.2976 

28 

0.8084 

43 

1.9063 

14 

0.3246 

29 

0.8538 

44 

2.0492 

15 

0.3524 

30 

0.9015 

45 

2.2159 

2.93 

TABLE  3a  (Con't.) 
n  -  60 


r 

■  r 

r 

1 

0.0167 

21 

0.4263 

41 

1.1321 

2 

0.0336 

22 

0.4520 

42 

1.1040 

3 

0.0309 

23 

0.4703 

43 

1.2403 

4 

0.0604 

24 

0.5053 

44 

1.2991 

5 

0.0063 

25 

0.5331 

45 

1.3616 

6 

0.1044 

26 

0.5617 

46 

1.4283 

7 

0.1230 

27 

0.5911 

47 

1.4997 

0 

0.1410 

20 

0.6214 

40 

1.5767 

9 

0.1611 

29 

0.6526 

49 

1.6600 

10 

0.1007 

30 

0.6049 

50 

1.7509 

11 

0.2007 

31 

0.7102 

51 

1.8509 

12 

0.2211 

32 

0.7527 

52 

1,9620 

13 

0.2419 

33 

0.7804 

53 

2.0070 

14 

0.2632 

34 

0.0255 

54 

2.2299 

15 

0.2049 

35 

0.8639 

55 

2.3965 

16 

0.3071 

36 

0.9039 

56 

2.5965 

17 

0.3299 

37 

0.9456 

57 

2.8465 

10 

0.3531 

30 

0.9091 

50 

3.1799 

19 

0.3769 

39 

1.0345 

59 

3.6799 

20 

0.4013 

40 

1.0021 

60 

4.6799 

*• 


T 


1 

0.0143 

19 

0.3140 

2 

0.028* 

20 

0.3336 

3 

O.Q435 

21 

0.3536 

4 

0.0584 

22 

0.3740 

5 

0.0736 

23 

0.3949 

6 

0.0889 

24 

0.4161 

7 

0.1046 

25 

0.4379 

6 

0.1204 

26 

0.4601 

9 

0.1366 

27 

0.4828 

10 

0.1530 

28 

0.5061 

11 

0.1696 

29 

0.5229 

12 

0.1866 

30 

0.5543 

13 

0.2038 

31 

0.5793 

14 

0.2214 

32 

O.6O49 

15 

0.2392 

33 

0.6313 

14 

0.2574 

34 

0.6583 

17 

0.2759 

35 

0.6861  _ 

12 

0.2948 

36 

0.7146 

2< 


TABU  : 


r 

r 

1 

0.0125 

21 

0.3023 

2 

0.0253 

22 

0.3192 

3 

0.0380 

23 

0.3365 

4 

0.0530 

24 

0.35*0 

5 

0.06U 

25 

0.3719 

6 

0.0775 

26 

0,3900 

7 

0.0930 

27 

0.4086 

8 

0,1047 

28 

0,4274 

9 

0.U85 

29 

0.4467 

10 

0.1326 

30 

0,4663 

11 

0.1469 

31 

0.4863 

12 

0.1614 

32 

0.5067 

13 

0.1761 

33 

0.5275 

14 

0.1911 

34 

0.5488 

15 

0.2062 

35 

0.5705 

16 

0.2216 

36 

0.5928 

17 

0.2372 

37 

0.6155 

18 

0.2531 

38 

0,6387 

19 

0.2692 

39 

0,6625 

2D 

0.2856 

40 

0.6869 

I  t 


2.96 


TJUUt  34  (don't.) 
B  •  90 


p 

p 

p 

p 

1 

0.0111 

26 

0.3387 

51 

0.8290 

76 

1.8310 

2 

0.0223 

27 

0.3543 

52 

0.8547 

77 

1.9024 

3 

0.0337 

28 

0.3702 

S3 

0.8810 

78 

1.9794 

4 

0.0452 

29 

0.3863 

54 

0.9000 

79 

2.0627 

5 

0.0568 

30 

0.4027 

55 

0.9358 

80 

2.1536 

6 

0.0686 

31 

0.4194 

56 

0.9644 

81 

2.2536 

7 

0.0805 

32 

0.4363 

57 

0.9938 

82 

2.3647 

6 

0.0926 

33 

0.4536 

58 

1.02a 

83 

2.4897 

9 

0.1047 

34 

0.4711 

59 

1.0553 

84 

2.6326 

IP 

0.1171 

35 

0.4890 

60 

1.0876 

85 

2.7992 

11 

0.1296 

36 

0.5071 

61 

I0I2O9 

86 

2.9992 

12 

0.U22 

37 

0.5257 

62 

lol554 

87 

3.2492 

13 

0.1551 

38 

0.5445 

63 

1.19U 

88 

3.5825 

14 

0.1681 

39 

0.5637 

64 

1.2282 

89 

4.0826 

15 

0.1812 

40 

0.5834 

65 

1.2666 

90 

5.0826 

16 

0.1945 

u 

0.6034 

66 

1.3066 

17 

0.2081 

42 

0.6238 

67 

1.3483 

18 

0.2218 

63 

0.6446 

68 

1.3918 

19 

0.2356 

44 

0.6659 

69 

1.4372 

— 

2D 

0.2497 

45 

0.6876 

70 

1.4848 

21 

0.2640 

46 

0.7098 

71 

1.5348 

22 

0.2785 

47 

0.7326 

72 

1.5875 

23 

0.2932 

48 

0.7558 

73 

1.6430 

24 

0.3081 

49 

0.7796 

74 

1.7018 

25 

0.3233 

50 

0.0040 

75 

1.7«43 

2.97 

TABU  y ft  (mb'!*) 


n  ■  ICO 

► 


r 

t 

r 

r 

1 

0.0100 

26 

0.2994 

51 

0.7082 

76 

1.4114 

2 

0.0201 

27 

0.3129 

52 

0.7286 

77 

1.4531 

3 

0o0303 

24 

0.3266 

S3 

0.7494 

78 

1.4966 

4 

0.0406 

29 

0.3435 

54 

0.7707 

79 

1.5420 

5 

0.0510 

30 

0.3545 

55 

0.7924 

80 

1.5896 

6 

0.0616 

31 

0.3688 

56 

0.8147 

81 

1.6396 

7 

0.0722 

32 

0.3833 

57 

0.8374 

82 

1.6923 

8 

0.0829 

33 

0.3980 

58 

0.8606 

83 

1.7478 

9 

0.0738 

34 

0.4130 

59 

0.8844 

84 

1.0066 

10 

0.1048 

35 

0.4281 

60 

0.9088 

85 

1.8691 

n 

0.1159 

36 

0.4435 

61 

0.9334 

86 

1.9358 

12 

0.1272 

37 

0.4591 

62 

0.9595 

87 

2.0072 

33 

0.1385 

34 

0.4750 

63 

0.9858 

88 

2.0842 

14 

0.1500 

39 

0.4911 

64 

1.0128 

89 

2.1675 

15 

0.1616 

40 

0.5075 

65 

1.0406 

90 

2.2584 

16 

0.1734 

a 

0.5242 

66 

1.0692 

91 

2.3584 

17 

0.1453 

42 

0.5411 

67 

1.0986 

92 

2.4695 

14 

0.1974 

43 

0.5584 

68 

1.1289 

93 

2.5945 

19 

0.2076 

44 

0.5759 

69 

1.1601 

94 

2.7374 

20 

0.2219 

45 

0.5938 

70 

1.1924 

95 

2.9040 

21 

0,2344 

46 

0.6119 

71 

1.2257 

96 

3.3040 

22 

0.2471 

47 

0.6305 

72 

1.2602 

97 

3.3540 

23 

0.2599 

44 

0.6493 

73 

1.2959 

98 

3.6874 

24 

0.2729 

49 

0.6686 

74 

1.3330 

99 

4.1874 

25 

0.2860 

3° 

0.6882 

75 

1.3714 

100 

5.1874 

Mtjmalmt  Find  tha  sapactad  waiting  tint  fbr  tbs  10th  failura  in  a  aaapla  of  alia  20. 
&ntna  that  0  •  1000  hours. 

SoMism  n»»w»  E(X]n.2nJ  *  (.66«>  1000  -  668.8  hours. 
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Solution i  Using  r  *  15  ,  n  ■  1r  ,  wo  get  E(Xr  n)  »  („li000)(2000)  »  800  hours 


TABLE  3(o) 


Ratio  of  Expected  Wad  tine  Tine  to  Observe  the  r**1  Failure  in 
Samples  of  Size  n  and  r  ,  Respectively* 


'r*n 


Wr1 


\n 

r>\ 

l 

2 

3 

5 

10 

15 

20 

1 

l 

o5b 

.33 

o25 

.20 

.10 

$5 

O 

o 

.050 

2 

l 

o56 

«39 

.30 

.ill 

.092 

.068 

3 

1 

.5 9 

M 

.18 

.12 

.087 

h 

JL 

.62 

.23 

.1*1 

.10U 

5 

1 

.2e 

.18 

.125 

10 

l 

.35 

.23 

Example:  Compare  the  expected  waiting  time  to  observe  the  10th  failure 
in  a  sample  of  size  20  with  the  expected  waiting  time  to  observe 
the  10th  failure  in  a  sample  of  size  10.  The  answer  is  *.23 
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Thble  4(a) 


Values  of  Li  M  (2r)E(X„  )/2r  for  *•  •  »01  . 
A*  ®  r,U 


2r 

3p  ~ 

4r 

5 r 

6r 

7r 

8r 

9r 

lOr 

20r 

1 

.005 

.003 

,003 

,002 

.002 

.001 

.001 

.001 

.001 

.0005 

2 

.043 

.027 

.020 

.016 

.013 

.011. 

.010 

.008 

.008 

.004 

3 

.089 

.055 

,040 

.031 

.026 

,022 

.019 

.017 

.015 

,00? 

4 

.131 

.079 

.057 

.045 

.037 

.031 

,027 

.024 

.021 

<01] 

5 

.165 

100 

,,072 

.056 

.046 

.039 

.034 

o030 

.027 

.013 

6 

.195 

.11? 

.084 

«o65 

.054 

.045 

o039 

.035 

.031 

.015 

7 

.21 $ 

.131 

,094 

.073 

.060 

.051 

o044 

.039 

.035 

<017 

8 

.2l»l 

.ll»3 

.103 

.080 

,065 

.055 

,o46 

.042 

.038 

.019 

9 

.260 

.155 

110 

.086 

.070 

.060 

0O52 

Oo46 

.04l 

.020 

10 

.276 

.164 

117 

.091 

.075 

.063 

o055 

„o48 

.043 

.021 

15 

.337 

.1 99 

,142 

.no 

.090 

.076 

.066 

.056 

.05? 

;02C. 

20 

.377 

0222 

<158 

c!23 

.101 

.085 

.074 

,065 

.058 

,026 

25 

°h06 

.239 

.170 

.132 

.108 

.091 

079 

=070 

,06? 

,030 

30 

.428 

.252 

.179 

039 

.114 

.096 

.083 

.074 

.066 

,03? 

ItO 

0U60 

.270 

,192 

.149 

.122 

.103 

.089 

.079 

.070 

.034 

50 

.482 

,283 

201 

056 

028 

.108 

.094 

.062 

.074 

036 

75 

.518 

.304 

.216 

067 

.137 

.116 

OOO 

.068 

.079 

,039 

100 

.540 

.316 

.225 

074 

.142 

.120 

004 

.092 

,082 

o4o 

P. nark i  Truncated  non-replacement  tests  of  the  fora  accept  If 

X  _  >  T  «  8  Li  (2r)E(I  )/2r  have  virtually  the  sane  O.C«.  curve  as 

r9n  o  i*  ^  Pi^l 

tests  of  the  fora  accept  if  9*  „  >  C  •  oXf  (2r)/2r  ,  In  the  above 

t9t\  o  JLm  pt 

table,  <*  (the  type  I  error)  is  .01  ,  when  8  ■  0q  , 


Spaniel  Consider  a  non- replacement  situation,  where  we  start  with  n  *  18  items 
where  r  •  6  ,  #o  »  1000  hours  and  <x  •  .01  ,  The  test  procedure  is 
accept  if  *6,18  >  117  .  In  wordst  Aocept  if  the  6th  failure  has  not  yet 
occurred  by  117  hours  and  reject  if  the  6th  failure  occurs  before  117 
hours  have  elapsed.  Such  a  plan  will  accept,  a  lot  with  mean  life  1000 
hours  with  probability  <>99* 
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Table  4(b) 

▼aloe*  of  X?  w  (2r)S(X  )/2r  tor  oc  »  .05  • 

j>*m>  r ,  n 


IS 

2r 

3r 

4r 

5r 

6  r 

7r 

8r 

9r 

lOr 

20r 

1 

.026 

.017 

.013 

.010 

.009 

.007 

.006 

.006 

.005 

.003 

2 

.104 

.065 

.048 

.038 

.031 

.026 

.023 

.020 

.018 

.009 

3 

.168 

.103 

.075 

.058 

.048 

.041 

.036 

.031 

.028 

.014 

4 

.217 

.132 

.095 

.074 

.061 

.052 

.045 

.040 

.036 

.017 

5 

.254 

.153 

.110 

.086 

.071 

.060 

.052 

.046 

.041 

,020 

6 

.284 

.170 

.122 

.093 

.078 

.066 

.05? 

.051 

.045 

.022 

7 

.309 

.185 

.132 

.103 

.084 

.072 

*062 

.055 

.049 

.024 

8 

.330 

.197 

.141 

.110 

.090 

.076 

.066 

.058 

.052 

o025 

9 

.548 

.20  7 

.148 

.115 

.094 

.080 

.069 

.061 

.055 

,027 

10 

.363 

.216 

.154 

.120 

.098 

.083 

.072 

.064 

oQ57 

.028 

15 

.417 

.246 

.175 

.136 

.112 

.094 

.082 

.072 

.065 

.032 

20 

.451 

.266 

.189 

.147 

.120 

.102 

.088 

.078 

.070 

.034 

25 

.475 

.280 

.199 

.154 

.126 

.107 

.093 

.082 

.073 

.036 

30 

.493 

.290 

.206 

.160 

.131 

.111 

.096 

.085 

.076 

.037 

40 

.519 

•305 

.216 

.168 

.137 

.116 

.101 

.089 

.079 

-039 

50 

.536 

•  315 

.223 

.173 

.142 

.120 

.104 

.092 

.082 

»040 

75 

.564 

*331 

.233 

.182 

.149 

.126 

.109 

.096 

.086 

.042 

100 

.581 

.340 

.242 

.187 

.153 

.130 

.112 

.099 

.089 

.043 

Reaark:  Truncated  toots  of  the  form,  accept  if  X_  „>  T  *  ©_  X2,  (2r)E(X„  )/2r 

r%  n  o  A-®©  rtn 

have  virtually  the  saae  O.C.  curve  as  teste  of  the  fora,  accept  if 
©r^a  >  C  a  ©o  X2-aC  (2r)/2r  .  In  the  above  table,  oc  (the  Type  I  error) 
is  .0?,  shea  9  ■  ©Q  . 

Kxanple:  For  a  ■  18  ,  r  *  6  ,  ®0  «  1000  hours,  and  cc«  .05,  the  acceptance 


region  is:  Accept  if  X*  >  170  hours. 
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Table  4(e) 

7a lues  of  X?  (2r)E(X_  )/2r  for 


IB 

3r 

4r 

5r 

6r 

7r 

8r 

9r 

lOr 

20r 

.053 

.035 

.026 

.021 

.018 

.015 

.013 

.012 

.011 

.005 

.155 

.098 

.071 

.056 

.046 

.039 

.034 

.030 

.027 

.013 

.226 

.139 

.101 

.079 

.065 

.055 

.048 

.04 2 

•038 

.019 

.277 

.168 

.121 

.095 

.078 

.066 

.057 

.051 

.045 

.022 

.314 

.189 

.136 

.106 

.087 

.074 

.064 

.057 

.051 

.025 

.343 

.206 

.147 

.115 

.094 

.080 

.069 

.061 

.055 

.027 

.366 

.219 

.157 

.122 

.100 

.085 

.074 

.065 

.058 

.028 

.386 

.230 

.164 

.128 

.105 

.089 

.077 

.068 

.061 

.030 

.402 

.239 

.171 

.133 

.109 

.092 

.080 

.071 

.063 

.031 

.416 

.247 

.176 

.137, 

.112 

.095 

.082 

.073 

.065 

;032 

.465 

.275 

.196 

.152 

.124 

.105 

.091 

.081 

.072 

.035 

.494 

.291 

.207 

.161 

.132 

.112 

.097 

•« 

0 

a> 

vs 

.076 

.037 

.515 

.303 

.216 

.167 

.137 

.116 

.100 

.089 

.079 

.039 

.530 

.312 

.222 

.172 

.141 

.119 

.103 

.091 

.081 

.040 

.552 

.324 

-230 

.179 

.146 

.124 

.107 

.094 

.084 

.041 

.567 

•333 

.236 

.183 

.150 

.127 

.110 

.097 

.087 

.042 

.590 

.346 

.245 

.190 

.156 

.132 

.114 

.101 

.090 

.044 

.604 

.354 

.251 

.195 

.159 

.135 

.117 

.103 

.092 

.045 

Reaark:  Truncated  toots  of  tbs  fora,  accept  if 

X_  >  T  ■  9  X?  (2r)E(X  )/2r  hare  virtually  the  oaao  O.C. 

r#n  o  i*ot  r«a 

A  9 

curve  as  tests  of  ths  fora,  accept  if  >  C  »  90  (2r )/2r 

la  the  above  table,  oc  (the  Type  I  error)  is  .10,  shea  9  •  9  . 


Sxaaple:  For  B«l8,r*6,90«  1000  hoars,  and  oC  *  .10,  the  aoeeptaace 
region  is:  Aocept  if  Xg  >  206  hoars. 
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Table  4(d) 

Values  of  X?  ^  (2r)E(X„  )/2r  for  oc  *  .25 

OO  r  fS 


X 

2r 

3r 

4r 

5r 

6r 

7r 

8r 

9r 

lOr 

20r 

1 

.144 

.096 

.072 

.058 

.048 

.041 

.036 

.032 

.029 

.014 

2 

.281 

.176 

.129 

.102 

.084 

.071 

.062 

o055 

.049 

.024 

3 

.355 

.218 

.158 

.124 

.102 

.087 

.075 

.067 

.060 

.029 

4 

.402 

.244 

.176 

.138 

.113 

.096 

.083 

.074 

.066 

.032 

5 

.435 

.262 

.188 

.147 

.121 

.102 

.089 

.078 

-070 

.034 

6 

.459 

.276 

.197 

.154 

.126 

.107 

.093 

.082 

.073 

.036 

7 

.478 

.286 

.205 

.159 

.131 

.111 

.096 

.085 

.076 

<>037 

8 

.493 

.294 

.210 

.164 

.134 

.114 

.099 

.087 

.078 

.038 

9 

.506 

.301 

.215 

.168 

.137 

.116 

.101 

.089 

.080 

o039 

10 

.517 

.307 

.219 

.171 

.140 

.118 

.102 

.091 

.081 

o040 

15 

.552 

.326 

.232 

.181 

.148 

.125 

.108 

.096 

.086 

.042 

20 

.573 

.338 

.240 

.187 

.153 

.129 

.112 

.099 

.088 

.043 

25 

.587 

.345 

.246 

.191 

.156' 

.132 

.114 

.101 

.090 

.044 

30 

.597 

.351 

.250 

.194 

.159 

.134 

.116 

.103 

.092 

.045 

40 

.611 

.359 

.255 

.198 

.162 

>137 

.119 

.105 

.094 

.046 

50 

.620 

.364 

.258 

.201 

.164 

.139 

.120 

.106 

.095 

.046 

75 

.635 

.372 

.264 

.205 

.168 

.142 

.123 

.108 

o097 

.047 

100 

.643 

.377 

.267 

.208 

.170 

.143 

.124 

.110 

.098 

.048 

Remark:  Truaeated  testa  of  the  fora*  accept  if 

Xp#ft  >  T  ■  (2r)S(Xr^a)/2r  have  virtually  the  same  O.C. 

curve  as  teste  of  the  fora,  accept  If  $r  #  >  C  «  (2r)/2r  . 

la  the  above  table,  oc  (the  Type  I  error)  is  .23,  shea  0  ■  0  . 

o 

Example:  For  a  ■  18  ,  r  »  6  ,  ##  «  1000  hours,  and  oc  ■  .25,  the  acceptaase 
region  is:  Aeeept  if  Xg  >  276  hours. 


\n 

r\ 

2r 

3** 

4r 

5r 

6r 

7r 

8r 

9  r 

lOr 

20r 

l 

.346 

.231 

.173 

.139 

.116 

.099 

,087 

.077 

.069 

.035 

2 

.489 

.308 

.225 

.177 

.146 

.125 

.108 

.096 

.086 

.042 

3 

•549 

.338 

.244 

.192 

.158 

.134 

.116 

.103 

.092 

.045 

4 

.582 

.354 

«255 

.199 

.164 

.139 

.121 

.107 

.095 

.047 

5 

.603 

.364 

.261 

.204 

.167 

.142 

.123 

.109 

.097 

.048 

6 

.617 

.370 

.265 

,207 

.170 

.144 

.125 

.110 

.099 

.048 

7 

.628 

"375 

.269 

.209 

,172 

.145 

.126 

.111 

.100 

.049 

8 

.636 

.379 

.271 

.211 

.173 

.146 

.127 

.112 

.100 

.049 

9 

.641 

.382 

.273 

:,212 

.,174 

.147 

.128 

.113 

.101 

.049 

10 

.647 

.384 

.274 

.213 

,175 

.148 

.128 

.113 

.101 

.030 

15 

.662 

.391 

.279 

.217 

.177 

.150 

.130 

.115 

.103 

.050 

20 

.669 

.394 

.281 

.218 

.178 

.151 

.131 

.115 

.103 

.050 

25 

.674 

.397 

.282 

.219 

.179 

.152 

.131 

.116 

.104 

.051 

30 

.67? 

.398 

.283 

.220 

.180 

.152 

.132 

.116 

.104 

.051 

40 

.681 

.400 

.284 

.221 

.180 

.153 

.132 

.117 

.104 

.051 

50 

.683 

.401 

.285 

.221 

.181 

.153 

.132 

.117 

.104 

.051 

75 

.687 

.403 

.286 

.222 

.181 

.153 

.133 

.117 

.105 

.051 

100 

.689 

.403 

.286 

.222 

.182 

.154 

.133 

.11? 

.105 

.051 

Remark :  Truncated  tests  of  the  fora*  accept  if 

X  >  T  m  0  X?  _  (2r)E(X_  )/2r  hare  virtually  the  saae  O.C. 

r,n  o  i-oc  rfa 

curve  as  tests  of  the  fora,  accept  if  9  „  >  C  ■  #X,  „  (2r)/2r  . 

r  ,n  o  i*ot 

In  the  above  table,  <*. (the  Type  I  error)  is  *50,  when  ©  *  9  . 

i  0 

Example t  For  a  ■  18  ,  r  *  6  ,  9Q  a  1000  aours,  and  oc  *  >50  ,  the  acceptance 

region  is:  Aceept  if  Xg  >  37®  hours. 


t 
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Table  5 

Vain**  of  r  and  X?  ,  (2r)/2r  auot 

*"#6 

®_  _  >  0AXf  _  (2r)/2r  la  aueh  that 
r  |H  o  !•« 


«•  ■  .01  OC  a  .01 

/3»  .01  £  *  .o$ 


"A 

r 

*1- 

>#t<2  r)/2r 

r 

X2 

*1- 

(2r)/2r 

*  ®v 

5/2 

136 

.8114 

101 

.7831 

2 

46 

.6892 

35 

.6492 

3 

19 

.5445 

15 

.4985 

5 

9 

-3897 

8 

.3633 

10 

5 

.2558 

4 

.2058 

oc  a 

.05 

O t  m 

.05 

0* 

.01 

0* 

.05 

3/2 

95 

.8374 

67 

.8079 

2 

33 

.7319 

23 

.6834 

3 

13 

.5915 

10 

.5426 

5 

7 

.4694 

5 

.3940 

10 

4 

.3416 

3 

.2725 

OC  a 

.10 

OC  a 

.10 

£>m 

.01 

/«  = 

.05 

3/2 

2 

77 

26 

.8570 

.7583 

52 

18 

.8269 

.7123 

3 

11 

.6383 

8 

.5820 

5 

5 

.4865 

4 

.4363 

10 

3 

.3673 

2 

.2660 

i 

Ot  a 

.25 

OC  a 

.25 

0m 

.01 

0  0 

.05 

3/2 

52 

.9033 

32 

.8758 

2 

17 

.8275 

11 

.7836 

3 

7 

.7261 

5 

.6737 

5 

3 

.5758 

2 

.4808 

10 

2 

.4808 

2 

.4808 

that  the  aeeeptane*  region 
L(0o)  «  1  -06  and  1.(9^  <  fi 


oc  ■  <01  v  ■  .01 

£a  .10  /5  »  .25 


*  Xl-«(2r)/2r 

r  X2.xC2r(/2r 

83  .7625 

60  .7244 

30  .6247 

22  .5715 

13  .4692 

10  .4130 

7  .3329 

5  .2558 

4  .2058 

3  .1453 

.05 

o<  a  .05 

/S  =  .10 

/3  *  .25 

55  .7890 

35  .7391 

19  .6548 

13  .5915 

8  .4976 

'  6  .4355 

4  .3416 

3  .2725 

3  .2725 

2  .1778 

OC  a  .10 

oca  .10 

/9c  .10 

0  »  .25 

4l  .8058 

25  .7538 

'15  .6866 

9  .6036 

6  .5253 

4  .4363 

3  .3673 

3  .3673 

2  .2660 

2  .2660 

OC  a  .25 

OC  «  .25 

0  a  .10 

fi  •  .25 

23  .8526 

12  .7932 

8  .  7445 

5  .6737 

4  .6339 

2  .4808 

2  .4808 

1  .2875 

1  .2875 

1  .2875 

2.106 


Example:  find  *  lift  test  which  possesses  the  following  O.C.  curve: 

If  the  mean  lift  is  ■  900  hours.  it  is  accepted  with 

o 

probability  .95;  if  tha  aaan  Ufa  ia  0^  ■  300  hours,  it  ia 
aeoaptad  with  probability  <  .10  . 

Solution:  In  this  axaapla  ®0/®^  *  3  <  oc*  .05  .  and  /5  *  .10  <  therefore 
the  required  number  of  failures  la  r  a  8  .  The  region  of 
acceptance  ia  given  by  n  >  (900) (.4976)  ■  448  .  In  words: 
Stop  life  testing  after  8  failures  have  occurred.  If  the  mean 
life  based  on  the  8  failures  that  have  occurred  >  448,  accept; 
otherwise,  reject. 


Table  6 

Values  of  a  ,  the  sanple  else,  needed  In  truncated  replacement 
procedures. 


OC  a  *01  a  .01 


3  5  10  20 


331  551  1103  220? 

95  158  317  634 

31  51  103  206 

10  1?  35  70 

4*  6  12  25 


oC  a  .01  a  .10 

189  316  632  1265 

56  93  187  374 

18  30  60  121 

7  11  23  46 


Oi  a  .05  a  vQl 

238  397  795  1591 

72  ,120  241  483 

23  38  76  153 

9  16  32  65 


oC  a  .01  /$  a  .05 


10  20 


237  395  790  1581 

68  113  227  454 

22  3?  74  149 

8  14  29  58 

y  4  8  16 


oC  a 

.05 

£  -  .10 

s 

.05  /S 

a  .25 

3/2 

130 

216 

433 

867 

77 

129 

258 

517 

2 

37 

62 

124 

248 

23 

38 

76 

153 

3 

11 

19 

39 

79 

7 

13 

26 

52 

5 

4 

7* 

13 

27 

3* 

4 

8 

16 

10 

y 

4 

8 

16 

1 

2* 

3 

7 

O i  a 

.10 

£  «  .01 

OC  a 

.10  fi 

a  o05 

3/2 

197 

329 

659 

1319 

128 

214 

429 

859 

2 

59 

98 

197 

394 

38 

64 

128 

256 

3 

21 

35 

70 

140 

13 

23 

46 

93 

5 

7 

12 

24 

48 

5 

8 

17 

34 

10 

3 

5 

11 

22 

2* 

3* 

5 

10 
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Table  6  (coat'd) 


od  * 

.10 

.10 

06  a 

.10 

/$- 

.25 

vs*A 

3 

5 

10 

20 

3 

5 

10 

20 

3/2 

99 

165 

330 

660 

56 

94 

188 

376 

2 

30 

51 

102 

205 

16 

27 

54 

108 

3 

9 

15 

31 

63 

5 

8 

17 

34 

5 

4* 

6* 

11 

22 

3 

5 

11 

22 

10 

2* 

2 

5 

10 

1 

2 

5 

10 

od  * 

.25 

■ 

.01 

06  a 

.25 

.05 

3/2 

140 

234 

469 

939 

84 

140 

28o 

560 

2 

42 

70 

140 

281 

25 

43 

86 

172 

3 

15 

25 

50 

101 

,  10 

16 

33 

67 

5 

5 

8 

17 

34 

3* 

5* 

10* 

19 

10 

2 

4 

9 

19 

2 

4 

9 

19 

06  « 

.25 

/S  ■ 

.10 

OC  u 

.25 

£*• 

.25 

3/2 

58 

98 

196 

392 

28 

47 

95 

190 

2 

17 

29 

59 

119 

10 

16 

33 

67 

3 

7 

12 

25 

50 

2 

4 

9 

19 

5 

3* 

4 

9 

19 

1* 

2* 

3* 

6* 

10 

1* 

2* 

3* 

5 

1* 

1 

2 

5 

*Ro— t»V:  it  was  Indicated  that  if  one  uaea  the  BQ  in  table  (5)  and  aeta 


the  aanple  aise  n 


r® 

i 

7 

( 


i 


f 


then  it  may  happen  that  while 


L(9  )  >  1-*  ,  L(0, )  say  be  slightly  >  P  .  One  way  of  getting  around  this 

ia  to  use  n+l  itesa  (rather  than  n  itesa)  and  to  use  the  slightly  aaaller 

truncation  tiae  T"  «  Bjt?  (2r  )/2(n+l).  The  teat  baaed  on 

o  o  o 

Min(I  T")  will  hare  L(0>  -  1-*  and  L(0,  )<fi.  In  the  shows 

r  iAu|  o  o  *  * 

o 

table  suoh  an  adjustment  had  to  be  made  in  the  following  eases: 
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Table  6  (coat'd) 


0t 

fi 

o/t 

0  0 

V°i 

(n+1) 

t”/t 

0  0 

.01 

.01 

% 

10 

4 

.959 

.01 

.05 

3 

10 

3 

.023 

.01 

.25 

3 

5 

4 

.959 

.01 

.25 

5 

5 

7 

.914 

oOl 

.25 

10 

5 

13 

.984 

.01 

.25 

3 

10 

2 

.654 

.05 

.05 

3 

5 

6 

.985 

.05 

.05 

5 

5 

10 

.985 

.05 

.05 

3 

10 

3 

.818 

.05 

.10 

5 

5 

7 

.976 

.05 

.10 

3 

10 

* 

> 

.818 

.05 

.25 

3 

5 

J 

.818 

.05 

-25 

5 

10 

2 

,889 

.10 

.05 

3 

10 

2 

.798 

.10 

.05 

5 

10 

3 

,887 

.10 

.10 

3 

5 

4 

.827 

.10 

.10 

5 

5 

6 

.918 

.10 

.10 

3 

10 

2 

.798 

.25 

.05 

3 

5 

3 

.962 

.25 

.05 

5 

5 

C 

J 

.962 

.25 

.05 

10 

5 

10 

.962 

.25 

.10 

3 

5 

3 

.962 

.25 

.10 

3 

10 

1 

.863 

.25 

.10 

5 

10 

2 

.719 

.25 

.10 

10 

10 

7 

✓ 

.958 

.25 

.25 

3 

5 

l 

0863 

.25 

.25 

5 

5 

2 

.719 

.25 

.25 

10 

5 

> 

.958 

.25 

.25 

20 

5 

6 

.958 

.25 

.25 

3 

10 

1 

.863 

Example:  Find  a  truncated  replacement  plan  for  which  Tq  •  500  hours, 
which  will  accept  a  lot  with  mean  life  *  10,000  hours  at  least  90  percent 
of  the  tine  and  reject  a  lot  with  mean  life  e  2,000  hours  at  least  90 
percent  of  the  tine. 

Solution i  In  this  ease  ©0  ■  10,000  ,  0  *  2,000,  <*■£«  »10.  Sinoe 
0o/8j«5,  m  */•■•  10,  we  see  from  Table  5  that  the  rejection  number  is 

Corresponding  to  ©e/0.ja5,  0o/To«2O,  *  */S  *.10,  one  sees  from  Table  6 

that  the  sample  sise  is  n«22. 
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Table  6  (coat’d) 


Thus  the  derived  truncated  replaoeaent  plan  meeting  the  requirements 
ie  as  follows:  Start  the  life  test  vith  n *22  iteas.  As  soon  as  an 
itea  falls  replace  it  by  a  new  item.  Accept  the  lot  if:  mia(X-  -_t500)  * 
500  (i.e.,  if  3  failures  have  not  occurred  by  500  hours,  stop  life 
testing  and  accept).  Heject  the  lot  if:  mind.  ,-;500)  «  X_  __  (i.e., 
if  the  3rd  failure  runs  before  300  hours,  stop  at  the  third  failure 
and  reject). 

Example:  Find  a  truncated  replacement  plan  for  which  T  ®  500  hours, 
which  will  accept  a  lot  with  mean  life  *  1500  hours  at  least  95#  of 
the  tiae  and  reject  a  lot  vith  mean  life  »  150  hours  at  least  95%  of 
the  time. 

Solution:  In  this  case  0Q  *  1500,  ©1  *  150 ,  </  o  /s  =  „05<, 

Since  0^/0^  =  10,*«/»  o  „05,  we  see  from  Table  5  that  the  rejection 
number  is  rQ  ■  3»  Corresponding  to  9^0^  «  10,  0q/To  *  3, 

0*/Si  .05,  one  sees  froa  Table  6  that  the  appropriate  sample  sins 
to  use  is  3*  Since  this  number  has  an  asterisk  (*)  attached  to  it 
we  see  that  we  can  actually  use  the  smaller  truncation  time 
T”  *  .8l8T  •  (.818)500  ■  409*  Thus  the  desired  truncated  replacement 
plan  meeting  the  requireoenta  is  as  follows: 

Start  the  life  test  with  3  items.  As  soon  as  an  itea  fails 
replace  it  by  a  new  itea. 

Accept  the  lot  if:  Mia(X_  _;409)  *  409  (i.e.,  if  3  failures  have 

Pi/ 

not  occurred  by  409  hours,  stop  life  testing  and  accept). 

Bejeet  the  lot  if:  Min(X?^;409)  «  X?  (i.e.,  if  the  3rd  failure 
occurs  before  409  hours,  stop  at  the  3rd  failure  and  reject). 


7 


M«d*d  in  trmaoatcd  npUentat 
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Tati*  ?  (coat'd) 


el  *  . 

10 

Am 

-10 

ft  m 

*10 

A  • 

-25 

VVt. 

Vn 

3 

5 

10 

20 

3 

5 

10 

20. 

3/2 

121 

186 

351 

681 

69 

107 

201 

389 

2 

39 

59 

no 

213 

21 

31 

58 

n3 

3 

12 

18 

34 

66 

7 

10 

19 

36 

5 

5 

7 

12 

23 

5 

7 

12 

23 

10 

2 

3 

6 

n 

2 

3 

6 

n 

rf  «  • 

25 

.01 

et  * 

.25 

A* 

*05 

3/2 

168 

261 

% 

965 

101 

156 

296 

576 

2 

51 

79 

289 

31 

48 

91 

177 

3 

19 

29 

54 

105 

12 

19 

36 

69 

5 

6 

10 

18 

36 

3 

5 

10 

20 

10 

3 

5 

10 

20 

3 

5 

10 

20 

«i  *  • 

25 

.10 

•f  ■ 

.25 

A* 

.25 

3/2 

71 

no 

207 

403 

34 

53 

101 

196 

2 

22 

33 

63 

123 

12 

19 

36 

69 

3 

9 

It 

27 

52 

3 

5 

10 

20 

5 

3 

5 

10 

20 

1 

1 

3 

6 

10 

1 

1 

3 

6 

1 

1 

3 

6 

tjgggjle:  Find  a  truncated  non-replaceaent  life  teat  for  which  T#»500 
hours,  which  will  accept  a  lot  with  aaaa  life  «  10,000  hours  at  least  90% 
of  the  tlas  and  reject  a  lot  with  as an  life  «  2,000  boars  at  least  90% 
of  the  tiae. 

Solatioai  Za  this  case  ®o»10,000,  0^*2,000,  «(■/  ••10.  Since  0o/0^«9, 

«t  m f «.10,  we  fiad  froa  Table  9  that  the  rejection  aaaber  is  ro»3  • 
Corresponding  to  0o/Tq«20,  *  «/» -.10,  oae  sees  froa  Table  7  that 

the  saaple  sise  a *23  .  Thus  the  derived  truacated  non-replaceaent  plea 
seating  the  requireaents  is  as  follows:  Start  the  life  test  with  n*23 
iteas.  Don't  replace  iteas  as  they  fail.  Accept  the  lot  if  aiadj^i^OO)* 
500  (i.e.,  if  3  failures  have  not  oocurred  by  500  hours,  stop  life  testing 
and  aeeept)2teject  the  lot  if:  aind^^^OO'K}^  (i.e.,  if  the  3rd  failure 
occurs  before  500  hours,  stop  at  the  3rd  failure  and  reject). 


2°  113 


Table  8 

7«la es  of  r  (upper  numbers)  and  of  ,(2r  )/2  (lover  nuabere)  such 
that  the  teat  baaed  on  uaing  a  aaapllng  plan  with  sample  else  equal  to 
(X  ^--t(2r0)/2p0]  and  with  rejection  nnaber  rQ  will  have  an  OC  curve  such 
that  L(po)*l-«*.  and  I^p^)  L(p)  ia  the  probability  of  accepting  a 

lot  having  fraction  defective  p. 


Vp. 

at  »  .01 

*  *  .05 

«*  =  . 

10 

.01 

.05 

.10 

.01 

.05 

.10 

.01 

.05 

.10 

5/2 

136 

101 

83 

95 

67 

55 

77 

52 

4l 

110.4 

79-1 

63-3 

796 

54.1 

43.4 

66.0 

43.0 

33.0 

p 

46 

35 

30 

33 

23 

19 

26 

18 

15 

31-7 

22.7 

18.7 

24.2 

15.7 

12.4 

19-7 

12.8 

10.3 

5/2 

27 

21 

18 

19 

14 

11 

15 

11 

9 

16.4 

11.8 

9*62 

12.4 

8.46 

6.17 

10.3 

7.02 

5.43 

T 

19 

15 

13 

13 

10 

8 

11 

8 

6 

> 

10.3 

7.48 

6.10 

7.69 

5.43 

3.98 

7.02 

4.66 

3.15 

L 

12 

10 

9 

9 

7 

6 

7 

5 

4 

5.43 

4.13 

3-51 

4.70 

3.29 

2.61 

3-90 

2.43 

1.75 

K 

9 

8 

7 

7 

5 

4 

5 

4 

3 

J 

3*51 

2.91 

2.33 

3.29 

1.97 

1.37 

2.43 

1-75 

1.10 

10 

5 

1.28 

4 

-82^ 

4 

.823 

4 

1.37 

3 

.818 

3 

.818 

3 

1.10 

2 

.532 

2 

•532 

Example:  Find  a  life  test  having  the  following  properties;  I 
accept  at  least  90%  of  the  lots  for  which  the  probability  of  failing  before 
sows  tine  Tq  ia  <  .01  and  will  reject  at  least  95%  of  the  lots  for 
which  the  probability  of  falling  before  Tq  >  .10. 

Solution!  In  this  problem,  po».01v  p^w.lO,  a  *.10  and  fi  «„05.  Thus 
p./p^.  m  10  and  so  we  see  from  Table  8  that  r  «2  and  H»[.532/.01]»53» 

A  O  O 

Thus  the  Ufa  teat  is  as  follows:  Place  53  items  on  test.  If  2  or  sore 
failures  oeeur  before  tiae  Tq,  reject.  If  one  or  fewer  failures  occur 
before  tiae  Tq,  accept. 


/ 


I 

1 
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Table  9 

Volute  of  h0#h^,  and  s  for  various  values  of  *  t  p  ,  and  9^.  The  normalised 

value,  9  si,  is  used, 
o 


©1*2/3  Q^l/Z  0J-1/3 


•< 

hl 

h 

0 

hl 

ho 

hl 

.01 

.01 

9.1902 

9.1902 

4.5951 

4.5951 

2.2976 

2.2976 

.01 

.05 

5.9714 

9.1078 

2.9857 

4.5539 

1.4928 

2.2769 

.05 

.01 

9.1078 

5.9714 

4.5539 

2.9857 

2.2/69 

1.4928 

o05 

.05 

5.8889 

5.8889 

2.9444 

2.9444 

1.4722 

1.4722 

91 

8 

2/3 

.8109 

1/2 

.6931 

1/3 

.5493 

Example:  Find  a  sequential  life  test  for  the  case  when  *  «.05, 05, 
©  «300  hours  and  6. =100  hours. 


Solution:  For  this  case,  h_=h, =1.4722  (since  0,  »  ~  if  0^  is  normalised 
■  ■■— ■ o  x.  x  }  o 

as  1  ),  and  s«.5493* 

Therefore  the  region  (55)  is  given  by: 

J00 c-l. 4722 +.549 Jr]  <  V(t)  <  300[ 1.4722 +.5495^  .  1 

Simplifying  this  gives: 

-442+l65r  <  V(t)  <  442+l65r  .  j 

The  life  test  is  continued  so  long  as  V(t),  the  total  observed  life  up  to 
time  t,  satisfies  both  inequalities.  As  soon  os  the  inequalities  are 
vieVat^di  ,  one  accepts  Ho  (i.e„,  0^300)  if  V(t)  >  442  ♦  l65r  and 
one  rejects  (i.e.,  accepts  (0^«1OO)  if  V(t)  <  -442  ♦  I65r  . 

} 

1 


t 


2. 


^  Table  10 

Approximate  values  of  Eg(r)  for  sequential  tests  for  various  values  o£ 

ks9o/0lj{rf  and  p 


*-V®i 


Mil 


e  /3 


5 


5 


5 


5 


oo  any 


3/2 

2 

.01 

.05 

.01 

.05 

12 

8 

7 

5 

12 

8 

7 

5 

62.4 

40.3 

23*3 

15.1 

60.4 

36.7 

22.6 

13.7 

128 

82.7 

43.9 

28.3 

82.7 

52.7 

28.3 

18.0 

47.6 

44.2 

14.7 

13-6 

30.8 

28.0 

9.48 

8.64 

0 

0 

0 

0 

10.4 

10.1 

17.5 
11.3 

5.00 

3,23 


6.74 

6.14 

11.3 

7.18 

4.63 

2.94 


Example:  rind  EQ  (r)  if  cne  ie  testing  =  3  with  *  =  »05»  ® 

c 

yfi  *.05. 

Solution:  The  expected  number  of  items  failed  in  reaching  a  decision 
when  0a9  is  (r)  a  2.94. 
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